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Conversion  factors  for  U.S.  customary 
to  metric  (SI)  units  of  measurement. 


To  Convert  Fro* 

To 

Mul  t ipl v  Bv 

angstrom 

me  t  e  r  *  ( m ) 

1 . 000  000  X  E  -10 

atmosphere  (normal) 

kilo  pascal  (kPa) 

i  .oi  3  2$  X  E  ♦  ? 

bar 

kilo  pascal  (kPa) 

1.000  000  X  E  ♦? 

barn 

meter'  (m*) 

1  000  000  X  E  -?8 

British  thermal  unit  ( t hermochemica 1 ) 

Joule  (J) 

1.0S4  3$0  X  E  O 

calorie  ( t hernorhemlcal ) 

Joule  (.1) 

4.184  000 

cal  (thermochentcal)  'cm* 

mega  Joule/m*  (Hl/s‘  ) 

4 . 1 8-.  000  X  E  -2 

curie 

giga  be* querel  (CBq ) • 

3.  700  000  X  E  ♦  1 

degree  (angle) 

radian  (rad) 

1.74$  3.'*  X  E  -2 

degree  Fahrenheit 

degree  kelvln  (K) 

t  •  (f  f  ♦  4$9 . b  7 ) /  1 .8 

electron  volt 

Joule  (J) 

1 .b0:  19  X  F  -19 

•'8 

Joule  ( .1 ) 

1 . 000  000  X  E  -  7 

erg /second 

watt  (V) 

1 . 000  000  X  F.  -  7 

foot 

meter  (m) 

3.048  000  X  E  -1 

foot -pound- f ore e 

Joule  (J) 

1.3$$  818 

gallon  (U.S.  liquid) 

meter*  (m  *) 

3.78$  412  X  F  -3 

inch 

meter  (m) 

2. $40  000  X  E  -2 

Jerk 

Joule  (J) 

1 . 000  000  X  E  49 

Joule/k  1  logr.im  (.l/kg)  (radiation  dose 
absorbed) 

Cray  (Cv'** 

1 .000  000 

k  i  lot  ons 

t  era  Joule* 

4.18' 

kip  (1000  Ibf ) 

newton  (N) 

4.448  222  \  e  *3 

kip/ inch*  (ksi) 

kilo  pascal  (kTa) 

b.894  7$ 7  X  F  ♦ ' 

ktap 

newt  on-  sec ond  / *• 

(N-s/m;) 

i.ooo  ooo  x  e  ♦: 

•  ic  ron 

meter  (m) 

1.000  000  X  F  -b 

mil 

meter  (m) 

2. $40  000  X  E  -$ 

•lie  (international) 

meter  (m) 

l.b09  344  XE  ♦! 

ounce 

kilogram  (kg) 

2.834  9$2  X  E  -2 

pound-force  (Ibf  avoirdupois) 

newton  ( N ) 

4.448  222 

pound-force  inch 

newton-meter  (N-m) 

1.129  848  X  E  -l 

pound-force/ inch 

newton/meter  (N/m) 

1 . 7 $  1  2b8  X  F  *2 

pound- force/ foot } 

kilo  pascal  (kPa) 

4.788  02b  X  E  -2 

pound-force/ inch*  (psi) 

kilo  pascal  (kTa) 

b.894  7$7 

pound-mass  (lbm  avoirdupois) 

kl log tam  (kg) 

4.  $3$  924  X  F.  -1 

pound-mass-foot*  (sioment  of  inert  ia) 

ki  logi  its-met  er  * 

(kg*m- ) 

4.214  011  X  F.  -2 

pound-mass/ foot  * 

k 1 logram/met  er  * 

(kg/m*) 

l .*01  84b  X  E  M 

rad  (radiation  dose  absorbed) 

Gray  (Gv)*« 

1 .000  000  X  E  -2 

roentgen 

coulomb/kl logram  (C/kg) 

2.  $79  7b 0  X  F.  -4 

shake 

second  (s) 

1.000  000  X  F.  -8 

•  lug 

kilogram  (kg) 

1.4S9  390  I  E  *1 

torr  (ws  Hg.  0*  C) 

kilo  pascal  (kPa) 

1  .  3  33  22  X  F.  -1 

•The  becquerel  (Bq)  is  the  SI  unit  of  radioactivity;  1  Bq  ■  l  event/*. 

••The  Gray  (Gy)  the  SI  unit  of  Absorbed  radiation. 

K  mote  complete  listing  of  conversion*  mav  be  found  In  "Metric  Practice  Guide  E  '80-74, 
American  Society  for  Testing  and  Materials. 
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SECTION  1 
INTRODUCTION 

TRAP-ML  is  the  result  of  a  dual  effort  to  (1)  simplify  the 
TRAP  computer  code  (Reference  1)  for  studying  radiation  effects  on 
aircraft  components  and  (2)  provide  an  analytical  predictive  and  cor¬ 
relative  tool  for  use  in  the  Defense  Nuclear  Agency  (DNA)  Tri-Service 
Thermal  Radiation  Test  Facility  located  at  the  Air  Force  Materials 
Laboratory,  WPAFb,  Dayton,  Ohio.  The  appendage  "ML"  on  the  code  name 
indicates  the  special  applicability  of  the  code  to  the  Materials  Laboratory. 

The  DNA  thermal  test  facility,  described  in  detail  in 
References  2,  3  and  4,  includes  a  test  setup  consisting  of  a  quartz  lamp 
bank  for  thermal  radiation  and  a  wind  tunnel  providing  convective 
cooling  on  material  samples  approximately  four  inches  square.  The 
samples  are  mounted  flush  with  the  inside  surface  of  the  wind  tunnel, 
facing  a  window  in  the  opposite  tunnel  wall.  The  sample  is  irradiated 
through  this  window,  simulating  the  environment  of  an  aircraft  exposed 
to  nuclear  radiation  while  in  flight. 

The  TRAP  code  (Thermal  Response  of  Aircraft  to  Nuclear 
Radiation)  includes  the  material  modelling  and  heat  transfer  capabilities 
required  for  the  test  facility  application.  The  solution  method  is 
based  on  a  finite  representation  of  the  cross  section  by  thermally  thin 
"elements."  One  or  more  elements  make  up  a  "layer,"  characterized  by 
uniform  material  properties  and  thickness,  and  several  layers  are 
possible.  The  layers  make  up  a  "segment";  several  segments  are  allowed, 
introducing  the  possibility  of  two-dimensional  heat  flow.  The  heat 
transfer  mec  isms  in  TRAP  include  radiation,  convection,  and  conduction. 

The  problem  with  the  original  TRAP  is  that  it  is  too  general 
for  the  test  facility  application.  For  example,  it  also  includes  a 
nuclear  source  generator,  atmospheric  transmission  and  ground  reflection 
factors,  and  a  vulnerability  assessment  including  an  iteration  on  slant 


range.  Also,  It  is  designed  to  analyze  large  complex  aircraft  structures  - 
not  laboratory  samples.  Clearly,  these  capabilities  are  not  needed 
for  the  test  facility.  Nor  are  the  stress-strain  calculations  needed, 
as  the  mechanical  boundary  conditions  are  not  prescribed  in  the  laboratory 
setup . 

Certain  other  capabilities  were  required  in  TRAP-ML,  however, 
which  are  not  in  the  original  TRAP.  The  quartz  lamp  radiation  source 
was  modelled  as  a  function  of  wavelength  and  time  and  included  in  the 
program.  Provision  was  also  made  for  two  other  source  options,  including 
a  point-by-point  option. 

While  the  test  facility  setup  can  be  treated  as  a  one-dimensional 
thermal  problem,  the  capability  to  analyze  two-dimensional  heat  flow 
remains.  For  example,  the  airstream  does  tend  to  "heat  up"  downstream, 
introducing  a  nonuniform  axial  temperature  field.  So,  rather  than 
eliminate  the  2-D  capability,  it  was  retained  for  possible  future 
applications . 

Also,  while  convection  and  reradiation  are  still  important, 
these  heat  *-ransfer  processes  were  modified  slightly  to  better  simulate 
the  test  facility.  Input  and  output  were  also  modified  to  enable  the 
user  to  operate  the  code  with  less  than  a  60,000Q  core  requirement  and 

O 

receive  temperature-versus-time  plots  at  each  surface  and  layer  boundary 
of  the  material  sample. 

TRAP-ML,  then,  allows  the  user  to  simulate  a  multi-layer, 
multi-material  sample  subjected  to  a  general  radiation  source  (a  6000 
watt  tungsten  lamp  in  particular) ,  complete  with  two-dimensional  heat 
flow,  convective  cooling  (if  needed),  and  reradiation,  with  both  the 
source  and  material  absorptivity  specified,  in  general,  as  functions  of 
wavelength. 

With  this  capability,  the  pretest  calculations  can  be  per¬ 
formed  which,  in  many  cases,  may  reduce  the  number  of  test  exposures 
required.  And  where  material  properties  are  uncertain,  post-test 
correlation  with  tests  involving  these  materials  may  lead  to  evaluations 
of  the  unknown  properties. 
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Section  2  describes  in  more  detail  the  heat  transfer  processes 
in  TRAP-ML,  and  Section  3  documents  the  computer  program.  Section  4 
includes  a  sample  problem  which  is  compared  with  a  test  exposure  of  a 
.032  inch  aluminum  sample  coated  with  3  mils  of  white  polyurethane. 

This  limited  correlation  indicates  excellent  agreement.  Finally, 

Table  1  includes  a  conversion  table  of  units  (English  to  Metric), 
since  the  TRAP  computer  code  (and  hence,  TRAP-ML)  was  originally 
formulated  in  the  English  system  of  units,  including  input  and  output. 
Program  listings  make  up  Appendices  A  and  B. 
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Table  1.  Conversion  of  units  (English  to  metric). 


Value  Expressed 
in  English  Units 

Multiplied  By 

Yields  Value  Expressed 
In  Metric  Units 

Inch  (in) 

2.54 

Centimeters  (cm) 

British  Thermal  Unit  ( BTU ) 

252.0 

Calories  (cal) 

BTU/in2 

39.03 

2 

Cal /cm 

Temperature,  °R 

0.5555 

Temperature,  "K 

BTU/in2/°R 

70.  25 

Cal  /cm‘'/°K 

BTU/pound  mass  (BTU/lb-m) 

0. 5555 

Cal /gram 

BTU/lb-m/°R 

1.0 

Ca 1 /gram/ °K 

SECTION  2 

HEAT  TRANSFER  MECHANISMS 

2-1  RADIATION  SOURCE 

The  primary  radiation  source  included  in  the  program  is  the 
quartz  lamp  bank  consisting  of  6000  watt  tungsten  filament  bulbs.  The 
spectral  distribution  and  a  general  time  history  are  included  in  a  non- 
dimensional  manner,  such  that  the  user  need  only  specify  the  peak  flux 
level  to  be  achieved  and  the  time  of  exposure  before  the  lamps  are  shut 
off.  Provision  is  also  made  for  other  source  characteristics,  or  even 
very  simple  sources  such  as  a  step-on  and  step-off  time  history. 

The  nominal  spectral  distribution  is  based  on  spectrometer 
measurements  and  the  point-by-point  model  is  indicated  in  Figure  1. 

While  there  Is  some  absorption  of  the  radiation  by  the  atmosphere, 
it  is  small  since  the  lamps  are  normally  located  within  about  2  inches 
of  the  tiierm.il  sample,  and  is,  therefore,  neglected.  The  distribution 
is  considered  invariant  with  time,  although  it  is  thought  the  peak 
intensity  may  shift  to  a  higher  wavelength  when  cool.  While  this  is 
Ignored  in  the  nominal  source  option,  the  capability  exists  of  specifying 
the  source  as  a  function  of  both  wavelength  and  time. 

The  program  numerical] v  integrates  the  distribution  between 
0.2  and  5.5  microns  to  determine  the  total  flux,  then  normalizes  the 
entire  distribution  to  that  value  as  a  function  of  14  discrete  wave¬ 
lengths.  If  wavelengths  other  than  the  nominal  14  are  specified  on 
input,  the  program  linearly  interpolates  between  the  inputted  values 
to  determine  values  for  the  14  basic  wavelengths.  These  fourteen 
wavelengths  and  their  associated  band  widths  arc  indicated  in  Table  2. 
When  material  absorptivity  is  specified  as  a  function  of  wavelength, 
these  14  bands  are  used  to  numerically  integrate  the  total  absorbed 
radiation  between  0.2  and  5.5  microns.  Since  all  the  Incident  radiation 


Table  2.  Wavelengths  used  to  determine  absorbed  flux 


Wavelength 

Number 

Nominal  Wavelength 
(Microns) 

Wavelength  Band 
(Microns) 

1 

0.288 

0.2  -  0.375 

2 

0.5 

0.375  -  0.625 

3 

0.75 

0.625  -  0.875 

4 

1.0 

0.875  -  1.125 

5 

1.25 

1.125  -  1.375 

6 

1.5 

1.375  -  1.625 

7 

1.75 

1.625  -  1.875 

8 

2.0625 

1.875  -  2.25 

9 

2.5 

2.25  -  2.75 

10 

3.0 

2.75  -  3.25 

11 

3.5 

3.25  -  3.75 

12 

4.0 

3.75  -  4.25 

13 

4.5 

4.25  -  4.75 

14 

5.0 

4.75  -  5.25 

is  either  absorbed  or  reflected  back  into  space,  the  incident  absorbed 
relative  intensity  can  be  given  by 


Y  abs 


a  (  X )  1(A)  dA 


I(X)  dA 


(1) 


where  a  represents  the  fractional  absorptivity  (0  ^  or  _<  1)  and  X  is 
the  wavelength. 

Figure  2  indicates  the  nominal  time-wise  radiation  history, 
where  the  flux  has  been  normalized  to  an  arbitrary  value.  The  time 
history  is  divided  into  two  domains: 

1)  The  "rise"  domain  includes  a  very  sharp  rise  in  flux  and 
then  reaches  a  nearly  steady-state  value  for  large  times. 
After  3.87  seconds  the  program  assumes  a  constant  flux. 

2)  The  decay  phase  commences  when  the  lamps  are  turned  off  and 
has  a  rapid  decrease  in  flux  until  it  finally  is  assumed  to 
be  zero  at  2.85  seconds  after  the  exposure  has  ended. 

Regardless  of  when  the  lamps  are  turned  off,  the  entire  curve 
is  normalized  to  a  value  of  peak  flux  which  is  specified  on 

input.  And  the  decay  phase  begins  exactly  with  the  time  specified 
on  input.  The  incident  absorbed  flux  can  then  be  given  by 


abs  .  .  -r  abs  —  .  . 

q  (t)  =  q  •  I  •  F  (t) 

^  max 


(2) 


where  F  is  the  relative  flux  defined  in  Figure  2  and  t  represents  time. 

The  user  can  also  specify  a  simple  step-on  and  step-off  time 
history  or  a  general  point-by-point  description.  Linear  interpolation 
is  used  between  tabulated  values.  Even  these  distributions  are 
normalized  to  the  largest  tabulated  value  so  that  the  intensity  of 
the  radiation  can  be  varied  with  a  single  input  parameter.  The  entire 
specimen  is  assumed  to  be  irradiated  uniformly. 
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nominal  radiation 


2-2 


RERADIATION 


Reradiatlon  of  energy  from  each  surface  of  the  sample  back 
into  space  is  taken  into  account  by  the  following: 


q 


rad 

i 


-A^e^oT 


4 

i 


(3) 


where 

is  the  rate  of  heat  flow  (BTU/sec) 

e.  is  the  emissivity  of  surface  i  (dimensionless) 

1  2  4 

o  is  the  Stefan-Boltzmann  constant  (BTU/sec-in  -°R  ) 

T  is  the  temperature  of  the  surface  (°R) 

1  2 

A^  is  the  surface  area  of  element  i  (in  ) 

The  emissivity  is  considered  "hemispherical  total"  emissivity  and  is 
a  function  only  of  the  temperature  of  the  material.  As  will  be  noted 
in  Section  3-1,  if  the  absorptivity  (a)  is  independent  of  temperature, 
then  the  emissivity  can  be  taken  as  l-« .  But,  in  general,  the  emissivity 
must  either  be  specified  on  input  or  calculated  by  the  program  using 
the  absorptivity  as  a  function  of  wavelength.  The  heat  flow  from  the 
surface  is  considered  lost  forever. 


2-3  CONDUCTION 

Conduction  among  elements  in  the  model  is  described  in  detail 
on  pages  86-92  of  the  TRAP  report  (Reference  1).  It  is  sufficient 
here  to  point  out  that  both  lengthwise  and  depthwise  conduction  are 
taken  into  account  and  are  functions  of  the  material  conductivities 
and  the  exposed  areas  of  adjacent  elements.  Figure  3  indicates  an 
example  of  a  model  consisting  of  two  segments  and  a  total  of  nine 
elements.  The  element  "nodes"  are  always  taken  at  the  geometric 
center  of  each  element.  In  this  case  there  would  be  conduction  between 
elements  1  and  6,  2  and  1,  2  and  6,  2  and  7,  2  and  3,  etc.  In  general, 
the  conduction  between  two  elements  is  given  by 


cond 

qij 


V 


(A) 


16 


where 


is  a  conductivity  factor  between  elements  i  and 
is  temperature,  °R 


j. 


BTU/sec-°R 


As  indicated  earlier,  .  depends  on  the  relative  geometrical  configuration 
of  elements  i  and  j,  and  also  temperature.  It  is  assumed  that  there  are 
no  conductive  heat  losses  to  the  surrounding  environment. 


2-4 


CONVECTION 


Convective  cooling  is  provided  by  means  of  exhaust  fans  at 
the  outlet  end  of  the  wind  tunnel.  The  model  is  located  in  the  "throat" 
of  the  tunnel,  a  region  of  constant  cross  section,  1"  wide  and  4  1/2"  high. 

The  boundary-layer  heat  flux  to  the  surface  is  computed  using 
the  Newton  heat-transfer  equation 


conv 

q 


h  (T  -  T  )  A 
r  w 


(5) 


where 

2 

h  is  the  heat  transfer  coefficient,  BTU/in  /sec/°R 

T  is  the  recovery  temperature,  °R 

T  is  the  wall  temperature,  °R 
w  2 

A  is  the  surface  area,  in 

Before  discussing  h,  it  is  noted  that  the  recovery  temperature 
is  computed  with  reference  to  the  stagnation  temperature,  Tq  (which  in 
this  case  is  the  ambient,  room  temperature),  and  the  temperature  at 
the  outer  edge  of  the  body  layer,  T ^ , 

T  =  T.  +  r  (T  -  T,)  (6) 

r  <5  o  6 

The  recovery  factor  is  dimensionless  and  taken  to  be  0.88,  based  on  the 
discussion  in  Reference  1. 
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The  free-stream  temperature,  ,  can  also  be  computed  from  T  , 
using  perfect  gas  relations,  where  a^  represents  the  total  speed  of 
sound  in  the  room. 

at  =  ^1.4  (1716)  Tq  =  /2402.4  Tq  (7) 

=  To  {  1  -  0.2  (V/at)2|  (8) 

V  is  the  free-stream  velocity  in  ft/sec.  The  free-stream  pressure, 
which  will  be  needed  later,  can  also  be  computed  assuming  a  perfect  gas, 

3. 3 

P5  =  Pq  |  1  -  0.2  (V/at)2  }  (9) 

where  P  is  the  ambient  (room)  pressure, 
o 

This,  then,  leaves  the  heat  transfer  coefficient,  h,  to  be 
determined.  This  parameter  is  a  function  of  many  factors,  but  depends 
primarily  on  the  type  of  flow  in  the  region  of  interest,  namely  laminar 
vs.  turbulent.  This,  in  turn,  depends  on  the  tunnel  geometry,  roughness 
of  the  surface,  any  leaks  in  the  tunnel,  etc.  Ideal lv,  one  probably 
would  predict  laminar  flow  in  the  test  facilit\  setup,  but  tin.  discussion 
in  Reference  4  indicates  the  flow  is  actually  till  1  \  developed  turbulent. 

Due  to  the  factors  mentioned  above,  this  is  plausible  and  .uiepled. 

Based  on  turbulent  flow,  then,  the  Colburn  equat  ion  (Referenn  '' 

is  used  to  take  into  account  dependence  of  h  on  wall  temperature  and 

air-stream  velocity  (two  velocities  are  possible,  depending  on  whether 

one  or  two  blowers  is  used).  This  requires  either  knowledge  of  the 

effective  origin  of  the  turbulent  flow,  or  some  empirical  information 

regarding  the  value  of  h  for  a  particular  case.  Reference  ‘t  provided 

”  A  2 

the  latter,  where  an  average  value  of  1 .  t.  7  x  10  BTU/see/inV°R  was 
determined  for  an  air  flow  of  793  ft/sec.  The  effective  origin  (x^^) 
was  then  back  figured  to  be  1.15  inches. 
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w 


The  Colburn  equation  and  other  necessary  relations  are 
presented  below.  The  user  Is  referred  to  Reference  1  for  a  more  complete 
discussion. 

0.026  •  k  •  R  '8 

h  .  _ (10) 

Xef  f 

where 


k  • 

[  1.0  -  8.89  UO'6)  T"  +  3.58  (10~8)  T'2  1 

L 

r 

_J 

(ID 

•  [2  (10-8)  T'1,5/(T'  +  201.6) 

R 

R  •  V  •  x  ,, 
e  ef  f 

m  -  --■■■■  1 

(12) 

eeff 

12 

R 

e 

3.7  (106)  P{  (T'  +  198.6)/(T')2*5 

(13) 

T"  - 

0.5  (T„  +  T6)  +  0.22  (Tr  -  T6) 

(14) 

The  program  also  provides  the  user  with  the  optional  capability 
of  specifying  the  heat  transfer  coefficient,  h,  which  allows  him  more 
direct  control  over  that  parameter. 

2-5  TEMPERATURE  CHANCE 

The  computer  program  combines  the  results  of  all  the  heat 
transfer  mechanisms  at  each  instant  of  time  for  each  element.  The 
temperature  change  can  then  be  calculated  for  the  finite  time  interval 
according  to 


AT 


At 

PACp 


q 


where 

AT  1*  the  temperature  rise,  °R 
At  Is  the  time  interval,  sec 

3 

,  Is  the  element  density,  lbs/in 


(15) 
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Cp  is  the  specific  heat  of  the  material,  BTU/lb-°R 
A  is  the  cross-sectional  area  (volume)  of  the  element,  in 
q  is  the  net  heat  flow  into  the  element,  BTU/sec. 
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SECTION  3 

COMPUTER  PROGRAM  DESCRIPTION 


3-1  PROGRAM  OPERATION 

The  TRAP-ML  program  was  written  in  FORTRAN  IV  and  developed  on 
the  Control  Data  Corporation  (CDC)  6600  computer.  The  program  requires 
card  input,  which  is  described  in  Section  3-2,  and  produces  time  history 
response  output,  the  subject  of  Section  3-3.  The  program  also  has  the 
capability  to  produce  a  plot  file  of  temperatures  versus  time.  This 
data  file  then  serves  as  input  to  a  second  program,  APLOT,  which  produces 
paper  plots.  This  second  program  requires  no  card  input  of  its  own. 

To  execute  and  also  generate  plots,  then,  requires  the  execution 

of  both  TRAP-ML  and  APLOT  back  to  back.  The  primary  reason  for  separating 

the  two  codes  was  to  enable  the  user  to  run  with  less  than  a  60,000„ 

o 

cell  core  requirement.  This  makes  it  possible  to  operate  the  code  from 
remote  terminal  under  current  WPAFB  restrictions.  A  typical  job  control 
setup  is  shown  in  Table  3. 

The  routines  for  both  programs  are  listed  in  Table  A.  The 
basic  plot  routines  are  included  in  the  APLOT  package,  rather  than 
calling  on  system  routines,  so  that  the  program  is  less  dependent  on  one 
particular  computer. 

TRAP-ML  requires  approximately  53,000  to  load  and  execute. 

O 

Logical  files  TAPE1  and  TAPE6  are  used  for  input  and  output,  respectively. 
TAPE5  is  used  internally  to  store  the  input  data  subsequent  to  reading 
the  data  into  the  computer.  The  plot  information  is  written  to  file 
TAPE8  and  is  automatically  rewound  prior  to  and  upon  completion  of  job. 
More  than  one  case  can  be  run  during  a  job  execution,  in  which  case  the 
plot  output  for  each  case  represents  a  separate  file  on  TAPE8.  Approxi¬ 
mately  8  cp  seconds  are  required  to  compile  TRAP-ML  using  "0PT=1"  and 
"R»2"  options. 
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Table  3.  Typical  job  control  sequence. 


Job  Control  Cards  Comments 


JOB  CARD . 

MAP, PART 

ATTACH, A, TRAPML .  )  Attach  and  execute  TRAP-ML 

A.  j  binary  file 

ATTACH, B,APLOT .  {  Attach  and  execute  APLOT 

I  if_  plots  are  desired. 

ROUTE, PLOT, .  Direct  the  plot  file  to  remote 

terminal 

EXIT. 

7/8/9 


!  input  data  for 
TRAPML 

6/7/ 8/9 
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TKAl’-Ml.  tiiul  Al’I.OT  rout  ini's 


Program  APLOT  only  requires  the  file  TAPE8  for  input,  but  does 
generate  output  on  TAPE6.  The  program  also  ultilizes  file  TAPE40  for 
internal  purposes.  If  on-line  plots  are  requested,  then  the  file  PLOTS 
represents  the  resultant  plot  file,  ready  to  be  directed  to  a  plotter. 

APLOT  requires  approximately  50,000g  cells  of  memory  to  load 
and  execute.  Compilation  requires  only  about  four  cp  seconds. 

Execution  times  will  depend  on  the  size  of  the  model,  the  time 
interval  required,  and  the  duration  of  the  response  requested,  but 
generally  program  APLOT  will  execute  for  only  a  small  fraction  of  the 
time  used  by  TRAP-ML.  The  example  problem  discussed  in  Section  3-4 
required  about  40  cp  seconds  for  both  programs. 

3-2  PROGRAM  INPUT 

TRAP-ML  looks  for  BATCH  input  on  file  INPUT.  The  input  data 
are  specified  in  groups,  where  each  group  begins  on  a  separate  card. 

More  than  one  card  may  be  required  for  a  group,  however.  The  variable 
type  and  format  corresponding  to  each  data  group  is  given  in  parentheses 
in  the  input  instruction  and  is  always  in  fields  of  12.  For  convenience, 
floating  point  numbers  can  be  left  justified  in  the  field  as  long  as  the 
exponent  is  right  justified.  Also,  zero  values  can  be  replaced  by  a 
blank  field.  Columns  73  through  80  are  not  used  for  data  and  can  be 
used  for  card  identification  or  other  purposes. 

All  input  parameters,  where  appropriate,  should  be  compared 
with  the  maximum  dimensions  provided  for  in  the  program,  as  delineated 
in  Table  5.  This  is  very  important  since  the  program  does  not  attempt 
to  check  the  input  for  all  such  violations.  In  Table  5,  the  input  para¬ 
meters  are  delineated  with  an  asterisk  (*)  preceding  the  variable  name. 

The  specific  instructions  for  input  are  contained  in  Table  6. 

The  remainder  of  this  section  will  attempt  to  amplify  on  those  instructions 
and  provide  insight  into  the  applicability  and  limitations  of  the  program. 
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Table  5.  Description  of  variables  which  determine 
program  dimensions. 


VARIABLE 

LIMIT 

DESCRIPTION 

*NAT 

10 

Number  of  temperatures  at  which  absorptivity 
is  given. 

NCONDC 

122 

Number  of  conduction  connections  between  elements. 

*NCPT 

10 

Number  of  temperatures  at  which  specific  heat 
is  given. 

*NCT 

10 

Number  of  temperatures  at  which  conductivity 
is  given. 

NEl" 

75 

Total  number  of  elements. 

NFALL 

8 

Number  of  points  at  which  the  relative  flux 
is  specified  versus  time  in  the  decay  phase. 

*NFLUX 

15 

Number  of  points  describing  thermal  flux  time 
history . 

NLAY 

5 

Number  of  layers  in  segment. 

NMAT 

9 

Number  of  different  materials. 

NRISE 

8 

Number  of  points  at  which  the  relative  flux  is 
specified  versus  time  in  the  "rise"  phase. 

*NSEG 

3 

Number  of  segments. 

*NTSPEC 

3 

Number  of  times  at  which  spectral  distribution 
of  radiation  source  is  specified. 

NWL 

14 

Number  of  wavelength  bands  used  to  describe 
radiation  source  in  calculations. 

*NWLA 

10 

Number  of  wavelengths  at  which  absorptivity  is 
specified  for  new  material. 

*NWLI 

50 

Number  of  wavelengths  used  to  describe  spectral 
distribution. 

NOTES : 

*  -  Variable  is  specified  on  input. 

//  -  Dimension  also  limited  in  subroutine  SETUP  by  a  fixed  point 
number  before  statement  number  23. 
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Table  6.  TRAP-ML  input. 


Group  1:  (40A2)  (ID(I) ,  1-1,  40) 

Title  card  -  free  field.  (ID) 

Group  2:  (2112)  INOUT,  IPLOT 

Program  output-option  code:  (INOUT) 

0,  do  not  print  input  data 
1,  print  input  data 
Plot  code  (IPLOT) 

0,  do  not  generate  plots  of  temperature  vs.  time 
1,  do  generate  plots  of  temperature  vs.  time 

Group  3:  (112)  NSEG 

Number  of  segments.  (1  £  NSEG  <_  3) 

Group  4:  (2F12.1)  XSEG(IS),  YSEG(IS) 

x-coordinate  of  end  of  segment,  in  (XSEG) 
y-coordinate  of  end  of  segment,  in  (YSEG) 


Note  -  See  Figure  3  for  definition  of  coordinates. 


Group  5:  (112)  NLAY(IS) 


Number  of  layers  in  segment  (1  £  NLAY(IS)  <_  5). 


IL,  and  are  to  be  repeated 


Group  6:  (F12.1)  TLAY  (IS.IL) 


Thickness  of  layer,  in.  [TLAY (IS, IL) ] 
Group  7:  (2112)  KMAT  (IS,IL),  NELL(IS.IL) 


Code  defining  material  of  layer:  (KMAT  (IS,IL)) 

1,  Aluminum,  2024-T3 

2,  Aluminum,  7075-T6 

3,  Magnesium,  AZ31B-H24 

4,  Titanium,  Ti-8Mn 
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Table  6.  (Continued) 


If  materials  other  than  those  listed  above  are  used,  they 
should  be  assigned  numbers  which  are  consecutive  beginning 
with  5.  Material  properties  for  the  added  materials  are 
called  for  in  Groups  8-21 

Number  of  elements  in  layer.  (NELL(IS,IL)) 


Groups  8-22  provide  material  properties  data  for  each  added  material, 
NM,  and  should  be  repeated  until  all  added  materials  have  been 
completed.  If  there  are  no  added  mater ials,  skip  to  Group  23. 


Group  8: 


(112)  MCOR(NM) 


Number  of  material  to  which  new  material  partially 
corresponds.  This  number  may  be  any  number  lower  than 
the  number  of  the  new  material  for  which  data  are  now 
being  supplied.  Use  of  a  corresponding  material  allows 
the  analyst  to  skip  reading  in  material  properties 
which  are  the  same  as  those  of  the  corresponding  material. 
If  this  number  is  read  in  as  zero,  all  material  properties 
for  the  new  material  must  be  supplied,  starting  with 
Group  10  and  ending  with  Group  22.  (MCOR(NM)) 


If  MCOR  *  0,  skip  to  Group  10. 


Group  9: 


(112)  KCH 


Code  indicating  material  property  to  be  changed: 

1,  absorptivity 

2,  conductivity 

4,  specific  heat 

5,  melting  temperature 

6,  heat  of  fusion 
11,  density 


Under  this  mode  of  input,  a  change  card  is  read  in  for 
each  property  to  be  changed,  followed  by  data  for  the 
specified  material  property.  Group  22  terminates  this 
sequence.  (KCH) 


Groups  10-13  provide  absorptivity  data  for  the  new  material,  if  required . 
(MCOR  -  0  or  KCH  =  1) 


Group  10: 


(112)  NAT(NM) 


Number  of  temperatures  at  which  absorptivity  is  given 
NAT(NM) ) 
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Table  6.  (Continued) 


Groups  11-13  provide  data  for  one  temperature,  NT,  and  are  to  be  repeated 
for  NT  =  1,  NAT(NM) . 

Group  11:  (2F12.1)  TALF (NT, NM) ,  EMI(NT.NM) 

Temperature  at  which  absorptivity  is  given,  °R 
(TALF (NT.NM)) 


Emissivity  at  specified  temperature,  dimensionless 
(EMI(NT,NM))  (Omit  if  integration  of  absorptivity  by 
program  is  desired  in  order  to  obtain  emissivity) 

Group  12:  (112)  NWLA 

Number  of  wavelengths  at  which  absorptivity  is  specified 
(NWLA) 

Group  13:  (2F12.1)  WL(NW,NT),  ALF(NW.NT) 

Wavelength  at  which  absorptivity  is  given,  microns* 
(WL(NW.NT)) 

Absorptivity,  dimensionless  (ALF(NW,NT)) 


Repeat  Group  13  for  NW  =  1,  NWLA. 


Groups  14-16  provide  conductivity  data  for  the  new  material,  if  required 
(MCOR  =  0  or  KCH  =  2) 


Group  14: 


Group  15: 


Group  16: 


(112)  KLD(NM) 

Code  of  conductivities:  (KLD(KH)) 

1,  lengthwise  and  depthwise  conductivities  are  the  same 

2,  lengthwise  and  depthwise  conductivities  are  different 

(112)  NCT(NM) 

Number  of  temperatures  at  which  conductivity  is  given 
(NCT(NM)) 

(3F12.1)  TCOND (NT , NM) ,  CONDL(NT.NM) ,  CONDD(NT ,NM) 

Temperature  at  which  conductivity  is  given,  °R  (TCOND(NT,NM) 

Lengthwise  conductivity  of  material,  BTU/in-sec-°R 
(CONDL(NT,NM) ) 


*If  emissivity  has  been  specified  in  Group  11,  only  wavelengths  from  0.2 
to  5.5  microns  must  be  covered.  If  emissivity  is  to  be  found  by  inte¬ 
gration  of  absorptivity,  wavelengths  from  0.2  to  about  250  microns 
must  be  covered. 
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Table  b.  (Continued) 


Depthwise  conductivity  of  material  (Omit  If  KLD(NM)  =  l), 
B'l’U/  in-sec-°R  (CONDD (NT ,  NM)  ) 

Repeat  Croup  1 b  for  NT  *  1,  NCT(NM) 

G roups  17-18  provide  specific  heat  data  for  the  new  material,  if 
required  (HCOR  ■  0  or  KCH  "  4) 

Group  17:  (112)  NCPT(NM) 

Number  of  temperatures  at  which  specific  heat  is  given 
(NCPT(NM) ) 

Group  18:  (2F12.1)  TCP (NT , NM) ,  CPM(NT.NM) 

Temperature  at  which  specific  heat  is  given,  °R 
(TCP (NT,NM) ) 

Specific  heat  of  material,  BTU/lb-°R(CPM(NT,NM)) 

Repeat  Group  18  for  NT  =  l,  NCPT(NM) 

Group  19  provides  melting  tempo  rat  ure  data  for  t  lie  new  material  if 
required  (MCOR  =_  0  or  KGH  m  5) 

Group  19:  (FI  2.1) 

Melting  temperature,  °R  (TMELT(NM)) 

If  the  melting  temperature  is  changed,  usually  all  of  the 
material  properties  will  have  to  be  changed.  The  reason 
for  this  is  as  follows:  the  me  ting  temperature  for  each 
material  for  which  data  are  pro  ided  in  the  program  is 
the  same  as  the  maximum  temperature  for  the  material 
property  data.  Suppose  that  tb  melting  temperature  is 
increased.  Unless  the  other  m.  .trial  properties  are  also 
changed,  the  melting  temperatu  will  be  above  the  maximum 
temperature  at  which  material  i  >perties  are  available. 

The  program  will  then  produce  ..  error  stop. 

Group  20  provides  heat  of  fusion  data  for  the  w  mat  or ial ,  i f  requi  red 
(MCOR  -  0  or  KCH  -  6) 

Group  20:  (FI  2.1)  HOF (NM) 

Heat  of  fusion  of  material,  BT1  b  (HOF(NM)) 


U) 
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Table  b.  (Continued) 


Groujy  21  provides  density  data  tor  the  new  mater  in lj  it  required 
(MCOR  -  0  or  KCH  -  11) 

Group  21:  (F12.1)  RHOM(NM) 

Density  of  material,  lbs/in'  (RHOM(NM)) 

Group  22:  Blank  card.  Is  required  for  each  new  material  NM,  either 

at  the  completion  of  Groups  10-21  for  MCOR(NM)  *  0, 
or  at  the  completion  of  the  sets  of  data  beginning  with 
Group  9  for  NC0R(NM)70. 

Group  23:  (4F12.1)  VF.L,  PO,  TAM,  HCONO 

Velocity  of  tree-stream  flow,  ft/sec  (VE1.) 


Ambient  (room)  pressure,  psi  (Pit) 

Ambient  (room)  temperature,  degrees  F  (TAM) 

Convective  heat  transfer  coefficient,  RTU/in*'-sec-0R  (HCONO) 


Note  -  Program  computes  coefficients  for  each  segment 
it  !U'0_NO=0.0. 

Group  24:  (112)  NOPT 

Code  designating  thermal  flux  time  history  (NOPT) 

1,  bOOO  watt  tungsten  source  (Figure  2) 

2,  Step  load. 

3,  Arbitrary  point-by-point  description. 

Group  25:  (FI  2.1)  QDMAX 

a 

Maximum  flux  achieved,  BTU/in^-sec  (QDMAX) 


If  NOPT=3,  skip  Group  2b 
Group  26:  (FI  2.1)  TCUT 

Time  at  which  thermal  source  is  cut  off,  sec.  (TCUT) 


If  NOPT-- 3,  skip  Groups  2  7-28 . 


Group  27:  (112)  NFLUX 

Number  of  points  describing  thermal  flux  time  history. 
First  point  must  include  time*0.  (2yNFl.UX^15) 
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Table  6.  (Continued) 


Croup  28:  (2F12.1)  TFLUX(l),  FLUXT(l) 

Time,  sec  (TFLUX(I)) 

2 

Relative  flux  corresponding  to  TFI.UX,  BTU/ln  /sec. 

Magnitude  of  flux  is  adjusted  to  QDMAX.  (FLUXT(I)) 

Repeat  Croup  28  for  1  =  1 ,  NFLUX 

Group  29:  (2112)  NWLI,  NTSPEC 

Code  designating  spectral  distribution  of  flux  (NWLI) 

0,  6000  watt  tungsten  characteristics  for  all 
times  (Figure  1). 

>0,  number  of  wavelengths  used  to  describe  distribution 
( 2 -NWLI < 50) . 

Number  of  times  at  which  spectral  distribution  is  specified. 
Assumed  to  be  1  if  NWLI-0.  (1<NTSPEC<3) . 

If  NWL 1*0  or  NTSPEC* 1,  skip  Croup  30. 

Group  30:  (bF12.1)  (TSPEC(I),  1-1,  NTSPEC) 

Times  at  which  spectral  distribution  is  specified. 

First  time  must  be  0.0.  (TPSEC) 

If  NWLI-0 ,  skip  Croup  31. 

Croup  31:  (2F12.1)  (WAVE(J,I),  QR(J,I)) 

Wavelength  at  which  flux  is  specified.  First  wavelength 
must  be  0.2  microns;  the  last  must  be  5.5.  (WAVE) 

2 

Relative  flux  corresponding  to  WAVE,  BTU/in  /sec.  Only 
relative  magnitudes  are  important.  (QR) 

Repeat  Croup  31  for  J«l,  NWI.I ,  then  repeat  that  sequence  for  1*1,  NTSPEC . 

Croup  32:  (3F12.1)  DELT1M,  TSTOP,  PRINT 

Integration  time  interval  (if  read  in  as  0.0,  program  will 
determine  time  interval  required  for  stability),  sec 
(DELTIM) 

Time  at  which  computations  are  to  stop,  sec  (TSTOP) 

Number  of  time  intervals  between  printouts  (0.0  will 
give  no  printout),  dimensionless  (PRINT) 
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Table  6.  (Concluded) 


To.  execute  another  problem,  repeat  Groups  1-32. 

Group  33:  (A3)  I END 

To  terminate  job,  place  "END"  in  columns  1- 


(IEND) 
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TRAP-ML  has  the  capability  to  assemble  a  two-dimensional 
heat  transfer  model  of  a  particular  material  specimen.  This  implies 
that  the  specimen  and  the  thermal  loads  are  uniform  in  one  coordinate 
direction  -  defined  here  to  be  the  z  direction.  The  x-coordinate 
direction  is  defined  as  in  the  direction  of  the  air  flow,  and  y  is 
defined  to  be  from  the  unexposed  surface  through  the  cross  section 
toward  the  exposed  surface  (see  Figure  3) . 

The  model  consists  of  segments,  layers  and  elements,  where 
elements  are  the  basic  "building  blocks"  and  must  be  thermally  "thin"; 
i.e.,  the  temperature  gradient  is  accurate  only  to  the  extent  that 
enough  elements  are  included  in  the  model.  Geometrically,  the  elements 
are  horizontal  slices  of  the  cross  section  of  constant  thickness. 

Layers  consist  of  a  set  of  elements,  all  of  which  exhibit  identical 
material  properties.  The  model  can  contain  up  to  five  layers  in  the 
cross  section  per  segment.  The  segment  is  defined  as  a  depthwise  set  of 
layers  representing  a  rectangular  portion  of  the  cross  section.  Elements, 
layers,  and  segments  all  have  constant  thickness.  If  only  one  segment 
is  included  in  the  model,  the  program  is,  by  definition,  reduced  to  a 
one-dimensional  analysis.  It  is  expected  that  only  one  segment  will  be 
used  to  model  the  test  specimens. 

In  the  general  model  indicated  in  Figure  3,  segment  1  has  only 
one  layer  and  five  elements.  Segment  2  consists  of  two  layers,  each 
having  two  elements. 

Now  for  the  input  instructions  in  Table  6.  Group  1  represents 
a  free-field  case  descriptor  which  appears  as  a  title  in  the  output. 

The  first  40  columns  will  also  appear  on  plots  if  requested.  It  should 
be  pointed  out  that  data  groups  1  through  32  are  repeated  for  as  many 
cases  as  the  user  desires. 

Group  2  defines  the  output  options.  By  selecting  1N0UT=1,  the 
output  is  reduced  somewhat  and  the  user  still  gets  a  raw  data  listing 
at  the  beginning  of  the  output.  The  plots  which  are  generated  if  IPL0T«1 
give  temperature  at  both  surfaces  and  at  each  layer  boundary  versus  time. 
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Croup  4  defines  the  location  and  size  of  each  segment  by 
specifying  the  lower,  left-hand  and  right-hand  coordinates  of  each  segment 
as  defined  in  Figure  3.  The  coordinate  origin  can  be  located  anywhere  - 
one  might  just  as  well  start  with  (0.,  0.)  for  the  first  segment  as  not. 
Groups  5  and  6  define  the  layers. 

Group  7  introduces  the  material  properties.  If  materials  or 
material  properties  other  than  those  associated  with  the  four  materials 
indicated  are  desired,  the  user  can  accomplish  this  by  calling  for  one 
or  more  additional  materials.  The  first  new  material  is  number  5,  the 
second  number  b,  etc.  The  properties  corresponding  to  the  four  built-in 
metal  alloys  are  indicated  in  Table  7  and  Figures  4  and  5.  In  this  regard, 
it  will  be  noted  that  the  selection  of  these  four  materials  was  arbitrary 
and  the  user  is  referred  to  Reference  1  for  a  complete  discussion  of 
property  specification.  It  should  be  noted,  for  example,  that  the 
absorptivities  specified  for  materials  1-4  are  all  0.5  for  an  unpainted 
condition.  If  the  user  wants  to  specify  a  different  absorptivity,  he 
would  need  to  call  for  a  new  material  and  proceed  with  the  input. 

Also  in  Group  7  the  user  must  decide  on  the  number  of  elements 
per  layer.  The  program  automically  divides  each  layer  into  elements  of 
equal  thickness,  depending  on  the  number  (NELL)  specified.  In  general, 
thicker  layers  with  lower  conductivities  will  require  more  elements. 

One  guideline,  other  than  checking  the  resulting  output  to  see  that  the 
temperature  gradient  is  not  too  large,  is  to  check  the  At  calculation 
(see  the  discussion  of  At  in  Group  32) .  Here  it  makes  some  sense  to 
end  up  with  a  At  requirement  more  or  less  the  same  for  all  elements. 

At  least  two  elements  per  layer  are  recommended  so  that  a  good  estimate 
of  temperature  can  be  made  at  the  boundary  between  layers  and  at  the 
surface. 

Group  32  requires  the  selection  of  a  At  to  be  used  in  the 
numerical  integrations.  This  number  must  be  small  enough  to  give  a 
stable  solution.  If  At=0.0  is  specified,  the  program  will  compute  a 
At  associated  with  each  element  and  then  select  the  smallest  of  these. 
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Table  7,  Material  properties  in  TRAP-ML 


Material 

Aluminum 

Aluminum 

Magnesium 

Titanium 

Property''^-^^^ 

2024-T3 

7075-T6 

AZ31B-H24 

Ti-8Mn 

Absorptivity 

0.5 

0.5 

0.5 

0.5 

Emissivity 

0.5 

0.5 

0.5 

0.5 

Heat  of  Fusion 
(BTU/lb) 

170. 

170. 

158. 

188. 

Density 

(lb/in^) 

0.100 

0.101 

0.064 

0.171 

Melt  Temp. 

(°R) 

_ 

1680. 

1680. 

1662. 

3495. 

Thermal  conduct 
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The  user  is  referred  to  Reference  1  for  a  complete  discussion  of  the 
formulation,  but  the  resultant  formula  used  is  the  following: 


At . 


1 


0.8  p  Aa 


(16) 


where  the  summation  on  j  extends  over  all  elements  which  touch  element  i. 
These  numbers  are  then  printed  out  for  the  user's  information. 

Groups  8  through  22  permit  the  user  to  specify  new  material 
properties.  Note  that  if  the  new  material  (called  for  in  Group  7)  does 
not  have  totally  different  properties  from  one  of  the  four  basic  material 
options,  then  oniy  those  properties  that  change  need  be  specified.  For 
example,  if  one  were  to  specify  the  absorptivities  for  7075-T6  Al,  then 
MCOR  (Group  8)  would  be  2  and  KCH  (Group  9)  would  be  1.  The  absorptivities 
would  be  specified  in  Groups  10-13. 

With  regard  to  melt  temperature  (Group  19) ,  the  program  treats 
the  melting  of  an  element  by  eliminating  it  from  the  solution.  The 
next  element  then  becomes  exposed  to  the  radiation  source,  with  the 
implicit  assumption  that  the  melting  always  proceeds  from  the  exposed 
surface  toward  the  opposite  surface,  element  by  element .  The  output 
reflects  this  event  by  setting  the  temperature  of  a  melted  layer  to 
0.0  in  the  printed  output,  and  a  small  vertical  arrow  indicates  the 
event  on  the  plotted  output. 

3-3  PROGRAM  OUTPUT 

TRAP-ML  output  consists  of  three  logically  distinct  sections. 

The  first  output  is  a  card  image  listing  of  the  entire  input  deck, 
exactly  as  it  is  read  in,  except  that  the  cards  are  numbered.  The 
second  section  describes  the  input  data  (if  IN0UT=1)  in  terms  of  its 
meaning  in  the  program. 

The  third  section  is  the  time-history  response  output.  At 
each  printout  time  the  following  values  are  printed:  time,  incident 
flux,  incident  fluence,  and  a  table  of  (1)  temperatures  and  (2)  a 
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breakdown  of  heat  flow  factors  for  each  element.  Figure  6  represents 
typical  time-history  output  for  three  consecutive  printouts. 

If  plots  are  requested  (IPLOT-1),  the  program  also  stores 
between  100  and  200  sets  of  temperatures  (evenly  spaced  in  time)  for 
both  surfaces  of  the  model  and  for  each  layer  division  for  segment  1,  or 
for  segment  3  if  there  are  3  segments  in  the  model.  Linear  interpolation/ 
extrapolation  is  used  to  estimate  these  temperatures  based  on  the  avail¬ 
able  values  at  the  center  of  each  element.  Each  curve  plotted  is 
identified  by  a  unique  symbol.  These  symbols  are  defined  in  Table  8. 

The  final  two  lines  of  printed  output  indicate  the  number  of  data  points 
per  curve  and  the  number  of  curves  written  to  disk  file. 

Error  messages  are  largely  self-explanatory,  although  the 
linear  interpolution  routine  (INT1Z)  only  indicates  a  code  number 
indicating  the  source  of  a  problem  if  a  variable  exceeds  the  limit  of 
the  table.  Table  9  identifies  the  routine  corresponding  to  the  error 
code. 

Program  APLOT  generates  a  list  of  the  data  plotted  and  indicates 
the  completion  of  each  graph  as  it  is  actually  generated. 
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Figure  6.  Time-history  output  for  sample  problem 


Table  8.  Plot  symbols  used  to  identify  plot  curves. 


Curve  Number 
(front  face  first- 
back  face  last) 

Plot 

Symbol 

1  (front 

face) 

□ 

2 

A 

3 

X 

4 

O 

5 

+ 

6  ^ 

t 

o 

Note  -  A  vertical  arrow  ( f  )  indicates  the 
time  at  which  the  element  melted. 
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Table  9.  Error  codes  from  INT1Z. 


Error 

Code 

Subroutine  which 
called  INT1Z 

Location 

Code  It 

1 

SPECT 

300'2 

2 

FLUX 

300 

3 

SPECT 

700_1 

4 

SETUP 

-3 

228  J 

5 

TSTEP 

12 

6 

TSTEP 

13 

7 

DTEMP 

7 

8 

DTEMP 

10 

9 

XHEAT 

62 

10 

XHEAT 

83 

11 

DTEMP 

25 

12 

SETUP 

117'1 

13 

TSTEP 

12-4 

14 

XHEAT 
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//The  location  code  is  read  as  follows:  a  code  S— ° 
is  interpreted  as  FORTRAN  statement  number  S  plus 
or  minus  n  cards. 


SECTION  A 

SAMPLE  PROBLEM  AND  COMPARISON  WITH  EXPERIMENT 

The  sample  problem  described  in  this  section  is  designed  to 
(a)  provide  the  user  with  a  check  case  for  exercising  the  code,  (b)  provide 
an  example  of  input  data  preparation,  and  (c)  make  a  limited  comparison 
with  experimental  results.  The  example  case  selected  is  based  on  an 
experiment  performed  at  the  DNA  test  facility  and  documented  in 
Reference  '4.  A  four  inch  sample  of  2024-T3  aluminum  substrate  .032 
inches  thick  was  coated  with  .003  inch  of  white  polyurethane  and  sub¬ 
jected  to  a  793  ft/sec  wind  tunnel  flow  and  a  two  second  pulse  from  the 
quartz  lamp  bank. 

The  input  data  for  TRAP -ML  is  indicated  in  Figure  7.  Although 
the  data  corresponds  to  the  test  data  as  closely  as  possible,  the 
documentation  was  not  totally  complete.  Figure  8  represents  the  plot 
generated  by  APLOT  of  three  temperature-time  histories:  (1)  the  exposed 
surface,  which  obtained  a  peak  temperature  of  1019<>R,  (2)  t he  division 
between  coating  and  substrate,  and  (3)  the  unexposed  surface  (backface), 
which  obtained  a  peak  temperature  of  900°R.  The  last  two  curves  fall 
virtually  one  on  top  of  the  other  due  to  the  relatively  high  conductivity 
of  the  aluminum.  Remember,  the  plot  symbol  codes  are  tabulated  in 
Table  8.  The  printed  output  in  the  vicinity  of  where  peak  front-face 
temperatures  were  computed  can  be  found  in  Figure  b. 

Figure  9  compares  the  backface  temperature  with  that  presented 
in  Figure  26  of  Reference  4.  Also  shown  are  the  results  of  the  ASTHMA 
code  from  Reference  4.  in  general,  TRAP-MI,  seems  to  predict  the  temperature 
response  well  -  its  only  shortcoming  is  a  failure  in  this  case  to 
predict  as  large  a  peak  temperature.  This  minor  deviance  could  easily 
be  explained  bv  any  number  of  uncertainties  in  the  data,  o.g.,  the  exact 
shape  of  the  radiation  pulse  time  history. 
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Figure  7.  Sample  problem  input  data. 
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TEMP.  C  DEG  R: 


•  HITE  P0LTURETHANE  0  VER  .132  IN.  2I2M-T3 
PEAK  FLUX,  BTU/1N2/SEC  *  .821 


TIME  l SEC! 


Figure  8.  Computer-generated  sample  plot. 
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TRAP-ML  predictions  to  experiment 
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appendix  a 

PROGRAM  LISTING  OF 
TRAP-ML 
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wK  K?l  I- 


P3OSPAM  T9flPM(_  <  IN°1JT  .OUTPUT  ,T  APF1=INPUT,  T  APEF  =  OUTPUT, 

1  TAPE5=129,TA3E«=5l4l 

tp43-v(l  IS  4  S 3E CT  A  L  VERSION  OF  TRAP  OE VELOPEO  p 0 »  THE  SIR  F09CE 
ISTroiflLS  LA30RAI0PY,  WPAFB,  3Y  KAMAN  AVIOYNE,  9UPLINGT ON , MA . 
Vc  R  r  T  ON  1.0,  RE3  T .  ,  1978. 

TE  M  3  FR A  T  URE  VS.  TT  ME  PLOT  INFORMATION  IS  WRITTEN  TO  TAPES. 

3  R  OGR A  M  APLOT  IS  3rOT'3EO  TO  GENERATE  THC  °LOT3. 


COMMON  A<75),ALr(10,lC),B(12?) ,CKO ( 75 > , CKO M ( 9 ) , 0 KL < 75 > , CKLM(°) , 

1  C0M0AG(l?2),C0NO0<lC,R>,C0Nni(l?2» , CONO J ( 1 22) , Cp ( 75 1 , CPMA TI9) , 

2  •)-LTIH,r)TIM(75>,DTIM£,OTMIN,EL<7R),  FLUXT  C 1 5 ) , MCONO, 19  <  2  0  )  , 

3  I :  L  <  3 )  ,IFLAY(5,3)  ,  I  FNO  ,  I M  A  T  (  7*5 )  ,  I  MELT  (75)  , I NOUT , I  PLOT , I SE G< 75 ) , 

4  I  STOP, JFLC?> , <CH ,KCONOIfl2?» .KCONOJ (122) , KCPI( 122 > ,KC3 J  <1 22) , 

5  <hat(3»5)»KjJR»LAvN(75),MC0R,NC0NDC,NEL,N£LL(3*5)  ,NEL513>  , 

5  Nr  LUX , NL  AY (  T ) , NMAT, NOOT,NPRTNT,NSEG,NTSPEC,NM.A,NWLI,  PO , 

7  3?INT,0<75) ,0A9( 3) ,0C0N0(75) ,QC0NV<3>  .QDMAX.JplU, QP (50 , 3) , 

8  ORAO,QRER(2,3),OPP(14)  ,SEGL  ,  T  (75)  ,T  AM, TOUT  ,TEr<-3)  ,  TEMP  (  751 
COMMON  TFLUX  < 1 3 >  , TFS,TIME,TLAY (3,5) , TP(75) ,TR,TSPEC <  3) , TSTOP, 

1  Vi  L.WAVE  (50 ,3)  ,WAVFL(14)  ,VIL(  10,10)  ,WT  (75)  « X  (75)  ,XC0N<75), 

2  X  S  EG ( 4 ) ,  XTR , Y ( 75) ,  YSEG (4) 

COMMON  /BLOCK/  A  LrAT(  io,  14,9)  ,CONOL  ( 1  0 , 9 )  ,  CP  M  (1  r> ,  9 )  ,  FMI  ( 1  0  ,  9)  , 

1  MOF  (9)  , KLO ( 9 ) , N AT (9) ,NOPT ( °)  ,NCT(9) ,NMATT, NWL, PI,RH0M(9) , 

2  MLF(10,9),T:ONO(10  ,9),TC°(10,9),TMFLT(9),WAVE8(15) 


CALL  tin 
5  CALL  SETUP 

IF  (  T'-TOP.EQ.l)  GO  to  100 
IF  ( OFLTIM.EO.O. 0)  CALL  TSTEP 
CAL  .  YH c  A  T  M0) 

OO  10  IE  =  1 , NE  L 
TFM3(IF)  =  TAM 
IME. T  (IE  )  =  0 

io  o < i - )  =  o.o 
ttm-  =  o. 
nprtnt  =  o 
OFLJ  =  0. 

CALL  XHEAT  (2) 

call  DTEHP(O) 

if  (IPLOT.GT.C)  CALL  RPL  OT ( 0 ) 

70  IF  (TPLOT.GT.0)  CALL  RPLOT(l) 
IF  ( MPPINT .GT. 0)  GO  TO  75 
IF  ( PRINT ,FO. 0.0 )  GO  TO  77 
CALL  PpRIMT 
NPRT  M T=P RINT-P . 5 
GO  TO  77 

75  NPRT NT=NPRINT-1 
77  TIMZ=TIME*DELTIM 

IF  ( T IMF, GT. TSTOP)  GO  TO  ao 
CALL  XHEAT  (2) 

CALL  OTEMP(l) 

OFLJ  =  :)FLU  ♦  3RA0*DELTIM 
GO  TO  70 

90  IF  (TPLOT.GT.0)  CALL  poL  OT  ( 2 ) 
GO  TO  5 

100  ENTITLE  s 
REWIND  S 
120  CT03 
EN9 


TPAPML 
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'M  )  :  K  PAT  A 


('  ’mm  pm  /PI  (1P.K/  M  t  AM  10,  14  ,  M  >  ,  cnNDL  (  1  C,  H*  .  ('“'mo  .  o>  ,rMT  (  1  ?f  <A|  , 
)  If  ('M,Kl  l(‘n,NAT(^|  .NC'flOl  ,NPT<H)  ,  NMA  T  T  ,  NHL  ,  t>  I  ,  $H(1*1  (  <5  >  , 

*  "  l  I  <  !  0, 9)  ,  TMIN'l  <  10  ,M>  ,  TCfM  1  0,91  ,  TMFl  T<  U  ,  WAVrl  <161 

A  Klfl'l  I  V  !  H, 

'MT  A  NAT  /<.  *  V  .  Al  r  A  T  0.  r-/  .  T  Al.  f  <  1  ,  1»  /  *0!)  .  0  /  , 

i  mi.  *  ,  1 1  /  i.,,ic. a  /  .  r  alf  <  1  ,?>/ .u  o.o/,  r  Air  <?,?»/  ihio. o/, 

’ML  (I.D/MO.O/.TAIF  (7 ,  Xl  /  1  1,1,?.  o/,  T  ALT  U  ,  41  /  90  0  .  0/  , 

>  T  A |  ‘  17,41/  94  9  4.  p  / 

P1'1  ”  t  f  i  i  mi  AT,  1TU/1  H-C1FT.  <?. 

MATA  MC 11 T  79,9,9, 9/ 

MAT1  lM  A I  1,  7074- T9. 

H  A  T  ('PM  /  0  .  t  74,0.191,0.  7  1 4  ,  0 . 7  7  1 ,  d  .  7  4  4  , 

I  '1 .  ’  4  0  ,0  .  79  7,  0  .  i  i  a ,  0  .  14  7  ,  n  .  9  9  7  , 

MAT  PM  A  l  7  C  7  4  -  T  6  . 

’  J  .  t  /  4 , 0  .  1  n  0 ,  0  .  7  1  3  ,  0  .  7  90 , 0 . 2r-  0  , 

fP.  ?M  ,0.  10  7,  0  .  M7,  0 . 99  9, 0.  9?n, 

MAT  PtAL  A ritn. 

4  0  .  M  r  ,  0  .  7  4  6  ,  0  .  7  6  4 , 0  .  7  4  1  ,  c  .  7  9  9  , 

4  p  .  1 !  7  ,  0  .  9  7 4 ,  0  .  116,0.9  Mi,  ,  (1 .  .1  Ml- , 

MAT  F--  I  A  l  4  ,  TT-4MN. 

■  0.  t  4  ,  .  l  ?/•, ,  (1.  141,,  0  .  I'M  ,  0.  l'-P  , 

ro.  i  >4,n. ?n4,o.  mi  ,  o .  <  ro ,  o.  a  r  r/ 

**  A  T  FM.U  1 

PAT  A  TC'*/  '60  .  0, 4  no  .  C,  h6  0  .P.S6C  .  0,  1060.0  , 

I I  .• ■> ' .  a,  t  i'« 6 o .  moo.  o,  i *' s j .  o , 
ff a r  i.'iai  >,  ru/‘,-r>,. 

7  l>,0.  fir, 460. 00, 6*0.  0  I',  060. 00, 1060.0, 

M  .0.1  At.  0.  0,  1  4  1-0  .  0,  1  6  4  n.  0  .  1  >'0  0 . 0, 

MAT  P  I A  l  A  ,  A  7  .<  1  n  . 

4  i>- ii .  p  i; , u t. p  .10,  i>!<  o .  o  i  ,  ” o .  n o ,  i  o 6  j  .  o , 

MM  .  0,  1  41.0.  0,  li.fM.  p,  1  bf,  7.  0  ,  1 1,  67.  0, 

MAT  p  T  A  l  4,  Tt-MN, 

<>  MO.  o  o, t>t.o  .no, 960 .  or ,  1 1 1, o.  o ,  1  760  . o, 

r\  A  Ml  .  0,  1  9  60.  o,  7  4  6 1,  .  o,  94  <14.  0  ,  .14  96, 0/ 
ft  .’mai  rnNiuir ti  vt  i  v,  ‘ATii/TN-'iF'i-nrr.  p. 

n  A  r  1  NCT/1C, 10,9,7/ 

MAT  P  T A l.  1  ,  70  74  -  T  .1  . 

MAT  A  ro  Mill  /I) .  001  ?4  7, 0.001. “7  4, 0.00  7767, 0.007471,0 .00  7440, 

1  P.  ll.  741  1 , 0.007  17.1,  0. 0  07  1  n<',  0  .p  07070 , 0 . 001494, 

“AT  F’lAl  7,  rO/4-16, 

*0,00  1  U  7, P.001MO, 0.0070  so, 0.007?  94,0.007177, 

M.  n  ? .14 n  ,  7 .00  7  74  C ,  0. 007P91 , U . 00700  0, 0 . 001  <96  , 

MAT  'MAI  1  ,  A  /  *  1  9  . 

4  0.1  r  •»  01  ,  0 . 0  0  0  M  0  ,  0 . 0  0  1  1  6 1  ,  0  .  0  0  U  ?6  ,  o  .  0  0 1 4  7 0  , 

■  '.  J  .  1  '•,14, 0  .00  1  6,9,  0 .00  IF,  9P,  0  .  0  01700,  0  .  0  01  7  00  , 

UA  T  PM  A I  4,  TT-imm. 

4.1. OP  0177,  '1*1.000411  4/ 

MAT  PM  A l  1  ,7074-T  ». 

hat  .  Tcnwn/1 6n.fi, t> 6n.o,96n. 0,960.0,1060.0, 

I  t  1  ,“.0,1  7 1.0. 0,1  9  60. 0,1460. 0,1 6  00.0, 

»AT  lM  A l  7,7074-16  . 

7  .160. Of, 6  6  0.00, 740. OP, 460. 00, 96  0.00, 

I I  P, ..  .  0,  1  740. 0,  1  4  6  7 . 0,  14;,r  .  0  ,  1690 . 0, 

MAT  PM  A  l  1 , A  7  1 1  1  . 

4  16 J.  P 6,460.00.6  4  0.  0  0,  96 p. 00 , 1060 . 0, 

41  4,1  .  O,  1  46C.  0, 1460 . 0,  1667.  0  •  1667.  0, 

MAT  PM  A l  4.TT-S4N. 

6.160.  P.o*  94*4,0/ 

I  «T ss  TV  I  IV. 

DAM  I  Ml  /4  0*0. 4/ 


M 


”1  ncKP 


HFJI  np  FUSION. 

DAT  *  HOT  /  5  70  .  0,  1  7C  .  0,  15°  .0,  16“  .  0/ 


r*r 'i ' !  ty  . 

OAT  A  3HOM/  C.  IPG,  0  «  1  r:  1  »  0  •  3fc  4  ,  C .  1717 

•1r'.T  tj  m  pr  da  T'J3  7  S  FOR  GTVFN,  fKGRFFS  ®. 

par  A  Tmfit/16*0. 0,1^  0.0, 166?. 0,346  5.0/ 

DATA  <Lrl/4*t/,  '1»&TT/U/ 

para  nw  /\  <♦/ 

DATA  WAtfF4/Q,2,  0.376,  9.625,  0.«76,  1.125,  1.375,  1.625, 

1  >.?6,  7.76,  3.76,  7.7 5,  4.25,  4.76,  5.75/ 

DATA  14169265  3F0v»«/ 

CMD 


1.675, 


Dl  OCKO 


ri.n-n'iTI  fir  COMV- :  <  I  FI?  ST,  t»0  ,TH,  VEl  ,  Tfmp,  EL  .  X'0\ ,  H1'  0N3 ,  T  r 

T  :  ,  QCOM> 


ruIc  ^OJTTN"  C4.C.M  AT  E  S 
'*03,  t  nc,  . 


Tur  CONVENTION  LOOSES  T^  to  ip.'-TPri!M 


ITI’rMSION  Tr^nj^t  ,  F  L  <  7  5  *  ,  XI"  0  N  (  7  5 ) 

IF  (  TfI-j?T  ,GT.  Ot  GO  T  n  40 
TA4  =  TIM  +  45  9. 6 0 
vr  =  AMCYI (VrL, ic  .  * 
rc  =  ri  < 

AT  -  SQ^T ( ?402. 4* T5 4 

Tr5  =  T<?»(  1.0-3.  2*  (W0/ftTt**?j 

0A-«  =  P  0 *  ( 1 ,  0  -  C  .  ?*  <  VO  /  AT  )  **  2 )  *•  3.  5 

T:  r  TFG  ♦  .?a*(T<?  -  TF  G ) 

X^Ff  -  3.15 
'\'r'lT  r  (4, ic*  T? 

1A  F  0  ? u  A  T  (ipco^ov-’V  T''1P-PiTl!P',  1-0.  R 
P-‘T  jp»j 


c15.  5  > 


oalgulatf  po.v/ErrTvr  moat  flow. 

4  0  Hr‘1M  =  hFOVO 

r F  f  wro*J0.  GT.0.0  )  GO  TO  6? 

tp  ?  L M  =  . 5* ( TF4° ( T  t  *  ♦  TF'.J  ♦  ,??*(TP-TFG* 

°'I1  =  3 . 7rF*3M‘  (  TOC  T  ‘14.1  03  )  /TCOJM*  *  ?  .  5 

4<=<1.Q-P,  •P9E-6«TP'>TMv-«.^»F_^.>TFPIM**2)*2.0E-8*TPPI'-l**i.t;/ 
1  (T59  T  ‘*+ ?  0 1  .6) 

P';3(in*«r*XFFF/i  2. 0 

4"  Vj  =  0.  j?4*a<*RF  *«.  0  ,  :  /y  r  Fr 

5]  ^CON  =  ipCN‘  (Tp’F^p  ,  IE)  >*EL(IFI 
53  T  J  PM 

Fi>n 


CPNVEO 


53 


G'|T  T!IT  I  N1'  0TF'i:>(<Cl 


OT--F  C  Al  -~IJL  A  Tf  ^  “FAT  FL01  ()Y  CJNotlOTTON  ANO  TNT  F  RNAl.  ’4^1  ATT  ON 

vd  opoa  h  s  rriJF5.Hu->"';  co”  ‘nrH  rt 

rOM^O”  3<7^l,£Lr<l''.i:»,T«l??»»OKn«Tci>,C<V4»0)t:Kl<7^),'KLNCf'», 

«  : H'p/j ;  << ,;)N"n  <1  r, 51 ,coNor  < i2? i  ,  ro no  m  i zz\ , cc  ( 7ci , ;f  *a t (o>  , 

?  o  m 'I  .O’-! '>«  7' ,nTMTN,  FL  <  7^  )  .FLLXT  (!<;),  HrONO,n<  ?  0»  , 
x  I  l  (?»  ,  »  FLAY<  *.  7  )  ,  IfN0,T>*fl  T  <  75)  f  TifiT  <73  I  ,  I NOUT , I  PLOT,  I  SF  Gf  75  )  , 

U  T  '  TpF,  JF  L<  3)  .  <N,  vr  O\o  j<  12?)  ,<rONOJ(l  22)  ,<CPT<  122)  ,K35J  (1  771  , 

x  <vr  t ( T ,f  > a VN (75) ,MCO*,NCONDC,NFL.NflI(3 , 5)  ,Nfi_5(  3)  , 
b  i:  i  'iv ,  ►'(.  av  (  »» ,  n'*at  ,  von ,  n*>c’t  nt ,  NF^r, ,n  Ts°Fr;  ,  nw_  A  , nhli,  =n , 

7  *  -  75)  ,  JA"  (  3)  ,  OOON.M  75)  ,  OCOMV  (  T)  tcnMAXtOrL!)*GP(50«  5)  « 

c  Pf  1,  )°f  <5(7,  31  ,  ,'~Fr,L  ,  T1  7>H  , TAM, TOUT, Tff< , TFUP(  *5) 

'•OUPN  rrinx(l-?»,Tr«:,riHF.TLflV(^,';»,TF(7S),T»fT;‘»tci,T»,TST0o, 

1  J  l  ,  Wrt  Vr  (5  1,  31  ,  .IAW>_(  1  l*»  ,WL  (  12,  1?)  , WT(75),X(75)  ,  XrON(  7  *  )  , 

■>  X  ~  FG((«  )  ,  XTP,  V  (  75  )  ,V  SFO  (<*) 

n0M  ,n»  1Lfhhj,14,',|,  OONOl  <  1  0  ,  *  I  ,0PM  (10 ,S  )  ,  F«T  (  1  0  ,  O'  , 

1  -n<-  (0)  ,  <10  ( "J)  ,  N  A  T  C°)  ,  MCr,T  (  Q)  ,  Nr*  T  (  9  )  ,  NMA  T  T  ,  N  Wl  .  P  T  ,  PMOM(  O  )  , 

>  T’lHlO.ot.T'VCMH  ,'H  ,70f  ( IT  ,9)  ,  TMTt  T(  1)  ,  WAV- 0  ( 15) 

inr\T'S  CrlTTOAL  TF"'PF(?AT  l»F  '~>Atio  IF  KD.oF.P 

SFf  I'"’  vairFIAL  CONDUCT TVTTj £5  ANO  SPECIFIC  H F fl T s 


TllsNFL 

oi  'r  ir=i,ru 

IF  (1VFLTI  TEl.FQ.l)  00  TO  70 
N“=T «AT (  TF) 

TFL=T"MP(  IF) 

Tr  (TrL,r>c,TMFi_r  <*“*))  go  to  ?7 

7  ra.'.L  INTI7  (  7,  t:l  ,  NOT  (n  n  ,  tfoNO  (1  ,NM)  ,nONiL  (1  ,NM)  ,C<L(IE)  ) 
TF  ( KID ( nm ) , fo. 7 )  go  to  i" 

0<  1  (  T  v  )  r  r  <1  (  I  T) 

GO  '0  <?: 

1C  F  1 1_ H'TI?  <  <?.  ,T  ~L  ,  NC.T  (\|1)  ,  Tr-C-in  (1  ,NM)  ,00Niin  <1  ,‘H)  ,r<0(T?  )  ) 
?:  IF  (  tf(_  .  L  T  ,TMF|_t  (*'*>- to.  0)  go  TO  ?5 

( r  r )  = ; .  i  *-nr  (  im  i 

GO  TOT 

25  C-'L'.  IN  i  1  7  (  11,  TFL  ,  KCPT  (NT)  ,T.'.o(l  ,MM)  ,t*PM  (1  ,NM)  ,DP(TF)» 

r  '■>  T  o  3  - 

2  7  J:\ljT  (N'*) 

OK. ( TF ) rOONDL ( I, NT) 

IF  (  »LD(N,').rO.?)  GO  TO  2- 
CKO ( Ti ) -OKI (  IF) 

GO  TO  2- 

?»  one  tf)  =  roNon  ( I  *  v* ) 

’-J  I-N'OT  CjT) 

GO  70  2  5 
rn  contt»(if 

OIL  3 1 L  A  T  F  CONn'JCTTVc  MFftT  C|_nn 
10 o  i  =  f'O  oNor 

oo  ‘  r  rcc=i, T“iu 

T  =  < '  n»'OI  (toe) 

J=<'.  O'/D  »(TC''i 

TF  <  T  Tf  L  r  I  T)»TTFLT(  J).Gt,1)  GO  TO  JJ 
T0  (F|'C1(  H".'  .FO,?)  GO  TO  U  0 

rVT  : r  KL  (  T ) 
go  to  5  : 

.■«o  ckt  -  Fvr'  (  n 

53  IF  (KCPvM  ICO.fO.?)  GO  TO  <-0 
CK  |  =  C<L  (  J) 

GO  TO  7  ’ 

F,  T  rT  J:FT5(  j) 

7-  OKI  )-rONOT  (  ICC)  /  OKT»CO’;OJ(  TOO /CXJtCONOAG(  TCO) 
r:(T'M.’(  n-rf(u(  j))/oktj 

O^O-JOtTl  -  O^-TNOIT)  -  OG 

o  c  o  n  17  ( j  >  =  icoNHj)  ♦  ,r 

0(  T ) =0  (  T)--K 
0)  J  )  -  F  (  J)FO~ 

■n  FONT  T  mjf 


OTFfiP 


5A 


'J3TiTf  rr-'Pr^Ty^r'S  (I'O  L  T  0  0  wr»T  T  I  0  NO 


)  £  L  T  T  M  /  (  WT  ( 


TC  ( Kr.r/ .  C)  GO  Ti  It," 

0"  1  1 ..  1^=1,  IJ 

if  '■’•'tin  rii.ra.ii  go  to  in 

M*=I  fi  ft T  (  T F  > 

T’*=  T  NrLr  <NM»  -13.0 
r"i  r  rr-io  (rri 

T-^l  T  r) :TF[I o (  Ii *  *  D  £  L  T  X 

TF  (  T.-'L.  GT.T**  -,0  to  10  u 
IT  IF  I  .  L  "  .  ™  )  GO  Tr  1 

Tl:'  I  TC  ).  (T'*-  TF  L  **WT  (  ID 

TF*0(  =  TI'1/(  HT( 

I r  (T'“-l  TE  »  .  L  T  .  tMFLt  C  'M  1 

T 0  =  T  G *- r,  (  TC| 

if  ( tf.mf . tec m ) >  go  to  i o o 

T«rLT  (  Tr>  =1 

JFLHS)  r  IF  L  <  T  5  )  *•  1 

IS*  -  NELSITG)  -  1 
It  =  UVH(T-) 

"l-LL  •  IS,  TL  I  =  'ni<Tr,IU  -  1 
T-1’<  TC)  =  0«  0 
0<  TF)=:.0 


TE*  *GP(IF») 


•:F(IF)  /" 

Ir)  *0  .  l*»  JF(  \IN»  ) 
GO  TO  10G 


IF  <MrLS( 
T  r  U  <  T  G )  = 
j-'_  f  r  g)  = 

Or  \|  r  T  ►'Ll  ■" 
r0NT  T  HU r 


'  )  . GT . j  I 


o  -  r  |  o  v 
r  MO 


OrrMD 


Oil.)  00.1  .KK)  i)  0CHio;)i>.>ilC).)U.K)00 


sm  rout  inf  cL(l<  (  t first  ,nopt,mflux,  TFiux,cL'JXT,iOMav,T!;  jr  ,t,pfao» 

THIS  ROUT  I N-  rA.CULATFS  the-  TNSTANTANFOUS  FLUX  AS  A  FUNCTION  OF 
ri^E.  THF»F  THPFF  OPTIONS. 


I  FI  >  <"  r  -  -3,  F I  ?5  T  PASF 

1,  ML  OTHFS  04SSES. 

NO3  T  -  l,  6003  <J  A  T  T  TUNGSTEN  LANP  CHA  9ACTF  RI STT  *5  . 

2,  STfp  r  U  NO  T I  ON . 

3,  AROITRAPY,  POINT  iy  POINT  (TCUT  NOT  JSEO)  . 
NFLJX,  tflUX,  FLUXT  -  OF  FINE  FLUX  T4«MF  VS.  THE. 
m-My  -  ppa<  cl J x 

TO  J  T  -  TTVF  ftT  XMTOH  TH  f  pMfl  l  FXPOSUPF  IS  OUT  OFF, 

T  -  THE 

OPfl 1  -  ImctAnt4VF0US  V4  LUr  OF  FLUX, 


01  If  MSI  ON  TFLUXdSI.  FLUXT(151,  TPI SE  (3  )  ,  0  OR  ISE  <  8  )  ,  TF4L  .  <  ft  1  , 
1  nF&LL<»> 

OflTa  NRISF/3/,  NFfl|.L/«/ 

oan  tr isf/o.,o.i 4, a. 22,0.32,0.47,0.72, 3. « 7, ioo./» 

1  30PTSF/Q.,4.9,5.5,5.‘»9,ft.l2,6.2  2,ft.4,6.4/ 

OATl  TFALl/Q.*.lft,«?3». 35, .54, 1.25, 2.85, 103.7 
CAM  00 FALL7 0.3, 2. 4,1. 5,. 33,. 47, .18, 2*0.0/ 

IF  ( IFIRST.GT.31  GO  TO  1 1C  0 
IF  <NOPT-?>  lC0»ftl0,700 


OPTION  1. 


IF  C  Tf?IFF(  II  .or.  TOUT!  GO  TO  30  0 
200  CONY  I  mijf 

TOO  CALI  INTt7{2,T:jT,NPISE,TPIFE,anRlSF,00‘i» 
OORAT  =  OOM4X/OOM 

onniitn  =  qon 
IT  =  II  -1 
OO  400  1=1,11 
TF.JXTIl  =  TRISECT) 

400  Ff.lX  Tfi)  =  ')ORISFCI)*OD:56T 
III  =  IT 

DO  5  P  o  I=l,NFftL'- 
III  e  in  *  1 

T^L JX (III)  =  TFftLLCI)  ♦  HUT 
500  fl  J v  T  (HI)  =  03FALL  CIHOORAT 
NFLJV  r  II  ♦  MF4LL 
GO  TO  100C 

OPTION  2. 

600  TfL  (XU)  =  0. 

TCL  I  X ( 2 )  =  TCUT 
TFLJXCII  =  TCUT  ♦  l.E-5 
TFLJVC4)  =  100. 

CL JXT (1 )  =  QOMDX 

FL JXT (2)  =  OOMjx 

FLJXTC3)  =  o. 

FL  J<  T ( 4 1  =  0. 

NFL  J X  =  4 
GO  TO  ICO” 

OPT  TON  3. 


71Q  304  =  0. 

no  *  0  0  I  =  1 ,  NF L 'J x 

IF  CFLUXTCII  .GT.OOMI  Q On  =  FLUXTTII 
8C0  CONTINUE 

(HRAT  =  QONAX/OOM 
00  -A  C  O  I=1,NFL'JX 
930  FLJXT(I)  =  FLUXMI  »  »OPRftT 


1000  JL  F  2 
Rct jom 


FLUX 


56 


FLUX 


57 


<ui>*oijTii>'r  inti?  ( Troop , x,nx ,x t,yt,d 

^MPOUTIMf  TMT1  IS  a  LINE*!?  INTFPPOtaTTON  O'nJTTXE". 

r,  T  7 :  M  a  VHJF  X,  TND  pcTU<?M<;  TM£  n ®RF SOONO T NS  VALUE  Y. 

I'TTf  -  rooo  r?DM  mil ng  pfoopa.m  fo°  tpasf3AC<  in  casr  of  cr',o=’.. 
nx  -  oiyenfidn  df  y-taplf  bno  y-t arl f  in  calling  prof,°ay. 

yr  —  tA^LF  or  X-VALUfS  Hi  CALLING  PROGP  A*. 

V T  r  TAOLr  OF  Y  -  V  a  L 1 11  *5  IN  GALL  INS  PROGPAN. 

P  =  Pt-TUPN  VALUE  *  . 


1 

? 

2 

4 


7 


4 

9 


OT  v  MSisn  *t ( 1) , vT  ( 1) 
tc  iy-xTiim,3,? 

:HlF(6,in  X,  (<T(  I>,T=1,NX) 

SO  TO  1J 

IF  (  X-xf  (NX)  I  4, Q,  3 
j~ITr(F,,it )  x,  (*TI  I) ,  1=1, NX) 

GO  Ml  1C 
00  3  I=?,HX 
IF  (V-VMIM  7,3,3 
CONr  t MUr 
P  =  Y  T  CD 

PFr JON 

p=YMT-1)*(X-XY(T-i))»(YT(I)-vt(I-1))/(XT(I»-XT(I-1)) 

pry J  o fi 

°=  xT » 1) 

B-T JON 
o-Yr  (NX) 


pi-r  i on 

wo  nr  (s,i2)  (xr  (i)  ,i=i,  nxi 
wttf  (3,i4>  icoor 

STJP 


FO  2  “  A  T  STA  7fiMr'lr  1  _ 


11  ft?-xt  iqm  x-vaijp  Fio.5,ix,i4mi)T5inE  taolc./ 

1  IX,  OPX - T A OL E  =. / (3 (2X,El?,o) ) ) 

12  FOB  4  a T  (PH  VTA1LF  =  /  (1X,3£14.S) ) 

14  f oR-a:  ( s^hc to \OFiarx  fops?  coor  in  tnti7  =  14) 

r»n 


TNT17 


S'I33PUT  INF  pplOT  (<P> 

TUI'.  -  OUT  T  NE  SETS  fJF  TMc-  °LOT  F7LFS  FO?  Tr^p  v3  TIME. 
r  PM'in‘J  A(76),ALr<lG,lC>,'m?2>,CKO(75>,ri<:}M(9),:KL<76),;KLM(c», 

1  s  omoag  (i??),:^M[n(ic>3».:oN?m??)tnoNTj(i?2),cc<75)f':'=MflT(p, 

2  1~  L  TI^  ,DTI'«<  75  )  ,  OTTMtr,  r,* -.1IN,  pL  (  7q,  ,  FLUXT  (IF)  ,  -(CONO,  I1(  2  3)  , 

3  I :  L (31  ,TFLA* (5, 3 ) ,IFN  O.TNAT (75) ,IMELT (75 ) , TNOJT,IOLOT, ISEGt75> , 

4  r- T"F, jfl( 3>  ,  <pM  ,  KCONOT  ( 122)  ,<COMOJ(l?2)  ,<CFI(122),tO»J(l?7), 

5  <M/>  7  (3  ,  5  ) ,  KS'J>,L  A*N  (  75  )  ,  i^O0  ,  NCONOC  ,  NFL ,  NF  LL( 3  <  51  .  MEL  5  (  3)  , 

=•  NrLL'X,  ML  AY  (  »)  ,  NMAt.N0Pt,NPPTNT.NSEG  ,NTSPFO , NH LA  ,NNLT,PO  , 

7  >?  I  NT,  0(  75  )  ,  1A  5(  3)  ,  CJCOMi  (75)  .OCONV (3)  ,Q3MA  X  ,  !)r  LtJ ,  OP  (5  C  ,  3 )  , 

"•  V  An,  )pro(  2.  3)  ,TPf  (  14  >  .SFGL,  T  (  75)  ,  U'l  .TOUT  ,  TF?  (  3)  .TEN3  (  75) 

~  0-M  ON  TFLUX(  15)  ,Trs,TIME.TLAY(3,5)  , T o ( 75  > , T p , TSPFC ( 3 > , TS TOP, 

1  /  r  l,W  A',r  (=50,  3)  ,  wAvrn  i  4)  ,Ul_  (  ic,  10  )  ,  WT  its  )  ,  x  (75  >  ,  XCON(  75  )  , 

2  X5FGI4)  ,xro,V(75)  ,  V  SE  ^  <^* » 

TT  -IT  ►,3IOK!  TCM|(,| 

N”1TMT  -  '‘AyiMjM  N'J“nFR  OF  POINTS  TO  PE  PLOTTFt. 

OAT  i  »'PO  INT/20C'/ 

Tr  (k-P-n  1!)C,40C,«CC 

100  NF'_  0  T  r  C 

M  r  (TSTOP  «•  OELTI  M)/0FLT1M 

MHO'  OT  r  3*N/N30TNT 

IF  (  M'JPLOT.EO.")  MIJPLOT  =  1 

srr  IIP  °LOT  FOP  FftCH  SURFACE  A  NO  LAYFQ  OO’JMOAPY. 

DO  3T0  IS=!,MSES 
ILJ  -  ML  AY  (IS) 

I-LAV(i,lS)  =  I  r  L  (  I  c ) 

00  ?  C  C  T  L=l,  IL'J 

?'0  tpl5  v  ( I  L  ♦  1 , 1  S »  =  IFLAYflulS!  ♦  NELL  C  T  *5 , 1 U  > 

300  IFL'-  v  (TLIJ+1,  IS)  -  IFL  AT  (  itUTl ,  I3>  -  1 

TO  =  j 

IF  (NSEG.GT.l)  IS  =  2 
TL  I  =  NLA* (TS) 

ILJ1  -  ILO  ♦  1 

Nop  t S  the  NUM3ER  Oc  SFGMrNTS  PLOTTFO. 

“'S3  ~  1 

NPirc  (■M  (10(1)  ,T  =  1, 20)  ,-'0MftX, TOUT, ILUl.NOPT, MSP 

NOT-  -  0 

PFT J  °  N 

400  Ic  (  m°l  n  t  .  GT  .  0 )  r,n  Tr  7C0 
NPL  0 T  =  NUPLOT 
N^T'  =  NPTS  ♦  1 
our;p  sjppfiOF. 

IF  =  IFL AY  (1,T3» 

T-3(l)  =  TFN°  ( T E ) 

TF  (tfm<1>  .EQ.0.0)  O’O  TO  440 
T r  (  mf.LL  (IS,  1)  .IF.  n  GO  TO  44C 

Tr3(l)  =  TFMP(I-:»  ♦  0 .5*  (TEMPI  IE)  -  TEMP(TE*1)  ) 

INN- c  SJPFAOF. 

440  I-  =  TELAY (ILU1, IS) 

TFN( TLU1 )  =  TEM° (IF' 

I p  (  TFM(  tluH  .£3.0  .  C)  SO  TO  433 
TF  ( NFLL ( IS, ILJ) .LF.l)  GO  TO  433 

TFN(IL'Jl)  =  TFN3  (IE)  ♦  .5*(TrMP(IE)  -  TFMO(IE-t)) 

L-AY'p  o011MPA°IE3  . 

4°0  TF  ( TLU . cO. 1 I  SO  to  650 
OO  6  To  I L  =  2, TL  J 
IE  =  IFL Ax  <IL,I5> 

TA  =  TfNr(IE-i) 

Ic  ( Tfl.EQ.O. 0)  SO  TO  520 
T1  =  TA 

IF  ( nfll (IS, IL-1 ) .LF. 1)  GO  TO  520 
TA  r  TA  ♦  .5* ( T  A  -  temd(TF-2>) 

520  TP  =  tf^pcie) 

TF  (Tq.ro, 0.01  SO  TO  530 
T?  =  xq 

IF  (MFLL(IS,IL).LE.1)  GO  tO  560 
TO  =  T«  ♦  . 5 *  ( T  3  -  tEmF ( IE ♦ 1 ) ) 

560  IF  (Tfl»TO.E0.0.3)  GO  TO  5»0 
TF  =  C.5*(rfl  ♦  toj 

IF  (T2.GT.T1)  03  To  565 
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PPLOT 


T  -  :M*!(r‘,T’| 

T-  z  ^I*'i  (r“,rj  > 

no  f  r  <;  7 

e~r  Tr  r  J'UV<  (t.-jk 
T"  :  I'll  M  irr,T'| 

CJ7-,  t-  4(  j  =  Tf 

vO  TO  F,  ..  r 

t-«ckj  r  £i4yi  (  Ta  ,  ▼<>  i 

e.  •!?  ~-mi 

«,c;  -  I'n'i'  T‘JUC 

Ml  TI«r,  JT-m<TI  ,  T=l,RUl» 
r j  =  molot  -  i 

o .-  T  |OM 

«v  f  \  i c  T  i  r  ? 

M'ITr  (e>»lCi«C»  I L  U 1 «  NP  T  s 

Tr.”  r  rnis 

1^.'?  ru<jT  (T^iiiMimrs  of  nUP>'P5  nFMFPATro  r  lu  f 
\  x  W  Ti')“  ’.fo  OC  rrjjfjTS  nu^vt  =  14) 

orr  r*] 

f'lO 


PPLOT 


60 


V.  [  N"  POO^IT 

0  01  '  05,  4(75)tALr(lG,lC),cMl??),!'Kr)(76),P<'),nq»,;;KL(75),rK:L*('’>, 

i  odmdag  a??> , :o\no du ,9 >  ,con->i  a  ??>  ,(  ono  m?2> , cp(  75) ,  o^s  T»n  . 

’  )  -  LH  1,0TT>M75>  ,  OM'-'F.,  0T4IU,  FL(  75)  ,PLUXT  115>,  HC0N0,  IK  ?  C>  , 

3  rim  , iflav  e 9 , x> ,  i  pui^at  ( 75) .  tyflt  (7-> > , inojt , ipi or, tgfg»  75> , 

4  I  o  TOP,  jr(. '  3)  ,  ,  KOOUin  1??)  ,<COfnj  ( 12?>  ,<CPI  (  1  ??)  ,<m  (l?’)  , 

5  <■*  t  T  (  3  ,  5)  ,  KG  J?  ,  l  4  VM  (75  )  ,  ''T'V’,  NCOMOT  ,  MEL  ,  NELL*  7,  5)  ,  \'Ft_  5  (  T>  , 

5  '|F  LUX  ,  f.L  AY  (  7)  ,  N^A  T,  igQP  I,  NDPI  MT,  N?PG  .UT'nFC  ,  MW  4,'WLTi3ni 

7  .  TN'T,  0  (  75  >  ,  -ja  rs  (  ■?»  ,000*0  (  76)  ,  iro*'V  (3» ,00*4X , 1plU,0P(5', 3» , 

*  '»:>Aq,-)OFP(2,3>,0‘?m4),SFGL.T(75>,TA‘.',T~!JT,T=;n7),TPYP(75> 

CO'YYflN  TFLUX(15),TF5,  TiM6,U_AY<3,5),TD<75),TP,TSot0(3),TFT7>P, 
t  -  L  ,  Wfl  Vf  (5  0,3)  ,w4VFL(14),WL(lC,10>,'IT(7  3  ),X(751,Xr:PM(7c;>, 

?  <  6  p  G  (  4  )  , XT  R , f (  75)  ,YSfG(U) 


H”  I  T 

'•■n  I  r 

O'* 

T3  - 
in 

T?  ; 

Tf"  ( 
IF  ( 
WF>Y  r 

Gn  r 
TP  ( 
\r  t  r 
l  l 
GO 

63  .(~Ir 
n  OOMT 
“"»■  J 


m,0FLu 


*3 


F  <5,10CC)  Y7k.f, 
r  (5,1100 
0  Ic=l,MFL 

TGFG(  t£) 

=  I  <T  L  (IS  I 
-  JrL ( IS ) 

IF.FO.TF!)  3P  TO  50 

It".EQ.  IF?)  GO  to  50 

P  (5.12P0)  IS, Tr,TFMP(T<n  ,QCONP(IE)  ,0(tf) 

n  mQ 

TE.f0.I5?)  }PCP(1,I~)  =  OP  FP ( 1, I S )  ♦  1PFP (?, IS ) 
c  (5,1300)  I S , IF, Tfmcjjf ) , nan<i5 ) ,npFP (1 , I~) , 

OMV  (  IS1  ,  OPD'IO  (  I F)  .0  (  or) 

'0  90 

P  (6,1400)  IF, IF, TPMnfxr ) , nPFPr ?, IS) , OOINO(TE) , OIT'I 
r  MIJC 
0,1 


1C00  Fii^JT  (13H0TX9F,  ~rr~  _  r1(,,F/ 

1  34H  P'CIDEMT  FLUX,  n  f  l  *  / 1  ►)*■  *  7/SFC  =  E14.6/ 

?  34H  T,irT0F9T  p  L  Uc  NIC r  ,  3TU/IM**?  =  614.6) 

no:  PP’YA  T  (  3 1  HO  SEGMENT  FLt  MEMT  TF WOF 5  A T Ua p , 9X , 

1  37hm  F4T  CL09,  otu/sfO/ 

?  2 7 P  NUH3FP  UUrfRF»  oegofEG  R,4X, 

3  4'.«4PS0pptI09  PFP4nIJTT0N  OONWECTICN  C  ONOUC  T  I  ON,  5  X  ,  X  K'CT  ) 

1230  ppp^AT  (1H  ,15, 6X , I ?, 4X ,Fll«4,41X,Ell.4,2X»E 11.4) 

1  300  FORMAT  (1H  , I6,5X,  I?,4X  ,E11.4, 5(?X,P11,4)  I 
14::  FORMAT  (1H  , 16, 5X. I?,4X,F11.4. 15X.F11 .4,13X, 2(?X,P11.4) ) 
fm3 


pppI NT 


S'Jl  ?  OUT  I  N,c  SFTJ3 

rOMNPN  A  (75)  ,  ALr  ( 1  0  ,  10)  ,  o<122)  ,GK0<?5 )  ,  CK1MC  S  ),CKL<75) ,  G<LNl°>  , 

i  gvioag  ( i??> , go n^d  a :,'-*> .conn t <i??)  ,oonoj(122> ,  cp«  7*  >  ,  gpma  t  (°» , 

>  V  L  Tl«  ,  "IT  I  Mi  ,  ntHlM,  PL  (  75)  .FUJXT  115) .HCONO,  T1  (  2  0  )  . 

T  i r  l  <  n  .  may  <=> ,  x » ,  1 3v-i,  inf  r  <  7«;»  ,im£lt  < ^ ,  indut  ,  i°lot  , TSEGcrs > , 

4  T~  TOP,  JR*  7)  ,  <fH  ,  KCONOI  ( 1 2  -*  *  ,  <C0N1 J  (1221  ,  <  GPI  (  1  22  1  ,  <3  3  J  ( 1  2?  *  , 

*  <■*  AT  (  »,  SI  ,  KSU’,  l.  a  VN  (7*  )  ,  M^OF'  ,  NCONOE  ,  NEL,  NrLL<  X.  5)  ,NEL*  (  X)  , 

F  N'  LUV,  NL  AY  (  XI  ,  N"e  T,  njPT.NPEINT-nSEG.N  T5PEG  ,N<lA  ,NWLI.  3  % 

’  3’  TUT  .o  (73  )  ,  3  a  3  f  x»  , 0CONO(  75)  ,  TC07'V  C  x»  ,00M4X  ,1- L'JtO®  (5  , 

?  V.A0,-jp«-Q(2,?)  ,1  Off  (14*  ,SEGL,T(75)  ,TA«,T'uT,Tr:-<,n  ,tenp(  751 

SOM  NOV  'Ftuxus*  ,Trs,Tl'«r,TLaY(3,5»,T0(75),TP,T3PrC(T),  TSTOO, 

1  4  L  .WAVE  (5  0,  T>  ,  Wi  VEL<  14)  ,  WL  <  13, 1P>  ,  WT  (75)  ,  X  (?R  ,XCON<  75  >  , 

7  X5rG(4) ,XTP,Y( 75) ,YSFG(4) 

erHHPM  7<u  05  </  ALC  AR  1C  ,  1-,,%)  .COUnL  <  10,  <5)  ,  CPM  (10,9)  ,EMI  ( 1  0,9!  , 

1  -OE  (Q)  ,<L3(0),  VAT  (O)  ,  vert  (91  , NCT(T>  ,  N  MA  TT ,  NWL  ,  P  I,  =>H  CM  (  3  )  , 

2  T\L‘r(lP,Q*,TG.)N,i(lC,5),TCP(lC,t)),T«ELT(9),Ma7E!5(l5l 
on- f’SJon  prco?(S)  ,Fi}k-r-t  (??| 

oara  ffff c n*'/icuaLiJMT  nu m ,  , io**fil  umi  num,  , iommasnesium  , , ioht it antum, 

1,  13H?C24-T3  ,10M7C75-7f>  ,10**  a7  3 1 R- H?4 , 1  0  HT I  -  *  MM  / 

nail  r’ff<rH/iOMa353ffPTI7I,13HrONOUC'CIVI,10HTHER>iaL  EX,13MSPECTCIC  H 
1,13  ( MEL  T  IMG  TF,  lC-HHrTar  nr  F!J  ,  1  OH  MOOUL  US  OE  ,  1  C HT=  N  S I L r  T  I , 

?1"-(T  F'iSILE  JL  ,  10  H’JL  TIMaTF  S.ICHOENSITY  ,  2HT  Y ,  2  H  TY , 7W=>  a  NS  I O*’ , 

J3H-)MCH,'Pr74TJ3'  .uhSION,  1GH  EL  1ST  ICTY ,  IOHElI  STRESS, 

4  1  3  H  f  Y  **A  TE  GTff,PHTffaTN),iH  / 

)AT'*  IF  M 1 / 4MfNl  /,  TE  M2/4H  / 

1  POR'*flT  ( FT  12) 

7  eq?«AT  (FF12.1) 
i  ft? ^ a t  (2ca4) 

RTOPrp 

?F43  (5,7)  (  in(I  )  ,  1  =  1,20 

T  c  (  T  o  (  )  )  ,NE.IEN1.0R.I0(2)«NE,TFN2)  GO  T  0  4 

TST3PR 

GO  T  n  6  4  0 

4  ffPl'ME,!)  IN1UT  ,  T  DLOT 

Ic  (INOUT.GT.O)  wc-PE  (6,351)  (TO(I) .  1=1,20) ,IPLOT 

osi  rn?4flr  ( l«l,?ca<»7/.TG“  rot  EDPE  (0  c0«  NO  pl6ts>  =  I?) 
oo  3  t=1,mwl 

5  wa7rL(T)  =  .F*(WC7EP)T)+wavri)(T4R) 

X  MI ''  =  10.  Or  10 

XMlXr-1 3 , 0  E 1 0 
TK=  1 
MEL  =  C 

NM a  I =  NMA TT 
<5)3=1 
NS  =  1 

RraKc,n  nseg 

I SU  =  NSEG ♦  1 

On  1 3  T5=l, ISJ 

?cao(5,?>  XSEG ( I S ) ,vSEG(IE) 

XMIMrflHIMl  (XMI'I,  V5EG(  15)  ) 

13  xmay=amaxi txMax.xsEGiis) ) 

T5J  =  1 5 IJ  -  1 

10  ? F>  IS=1,ISU 

<?rn(5,D  nlay<is> 

IlJ: MLfly (IS) 

no  14  IL  =  1,IL’J 

?raj  (5.3)  Ttayds.TL) 

RrA  3 (5,1)  KM4T(I5,iL) ,NtLL(I5,TL) 
n>«a  r  =  «axo (NMar,<MaT (is, il> ) 

14  GONTTNUE 

GALOULATE  SEGMENT  LFNGTHS 

1X=Y5EG( IS»1) -X5EG ( IS) 

Ov=Y  5FG( IS  +  1) -YSEG ( IS) 

5rGL=SlPT(PX*»2»nv**2) 

OXrjy/s-GL 
DV  =  3Y/SGGL 
ixL=nx 
iyl=pv 

SET  UP  ELEMENT  NUMlEffING  ST S TEM, T AL CUL ATE  ELEMENT  THIGXNFSS, 
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SETUP 


4  '£  ft 


10  HT,  Li  POSITION 


I r L  <  T  ?)  =  MFIM 
09  <"=  TL  =  1,tlj 

rMi'^=TLftvCIS,IL)'FLOflT(NrLL(IS,TL»» 
a»-  4  =  THiCK*"cr,i 
TF  <tl.3t.d  r,o  t"  ?? 

Irr - i fl (TO) 

V<I-T)=i.-»tyf_:'(j:)+Y-cr,,I^  +  1)  +  THrC)<#r,YLt 
T<I-T)=j.S*(v^rr,  (TSH-¥S''0»IS«-i)-T4Tr;i<*CUL> 

2?  IrL.  =  ^1  +  3 

Tr_  )-vFL  +MFLL (13, <L' 

T-  (I  FL'J.Lr.TF)  30  Tn  ?3 
N  ~  T  T  r  <F,c4ll  M  3 
c> T  1 3 

’3  Do'  ?u  Tr=I'LL»TrLU 
LIY-'MIF)  =  IL 
T'MT  /  IF)  =  K«AT  <  T3  ,  TL  » 

T o -  3 <  T  £  I  =  TC 

T  (  T  )  =T  u  TC  i 

EL  <  I r»=0TGL 
\  <  I  '  )  =  4  '•  F  4 

tp<  r  f  ,E0  .TFT)  33  T9  ?u 
<  Cl  :  )  =  <<  Tr-1  )  *0. 5-  <  T(  IE-1>  +  T  <  I!  E>  ,  *9YL 
Yd  -  )=v<  iF-ii-a.o*  (T(tf~d+t(if)  i*txl 
z\  p'TMTTNUr 
?3  N-  L  =  T  *L ! J 

MFi.  0fTFI=TFLI|4-l-IrL<TF) 

JCL  (  IS)  =  Mil 

f'PNlTTN'Ui 

I o  <  I  MO <JT  ,  El  •  0  •  00  TO  23 
^nTMndUl  M3rG 
roj  =  nseo  ♦  i 

U’rr(f.,or.M  <  XSEM  IS  > ,  10=1,  ISU) 

W~ITF<b.F:5)  < YSEG  < IS  I , IS  =  1, IOU) 
rj  ~  TGU  -  1 
00  d?  TCrl.IdJ 
IL  J  =  H  L4  Y  (IS) 

W’I't  (6,3!t)  T  S  » I L  U 

w»I  T  F  <  <3  ,  J  1  J  (IL,TL4Y(lS,iL),<YflT»IS,IL>,NELL<IS,IL)  ,TL=1,TU)> 

912  0ONr  tk||J£ 

0AL3ULATE  HP  AT  iONO'IOTION  CONSTANTS  F0°  ELEMENTS  WITHIN  99ANOH 

23  r)o  +0  T  0=1,  ISJ 
103  =  1 

T‘rL_I  =  IiL(I0> 

IEl JI=IELLT-1+NELS  (IS) 

IF  (IS. IT. ism  30  TO  31 

IF<i(SFG»ll  .NE.X3EG (IS+l)  .0» . YO EG( 1 > . NE . YSEGl IS*1 »)  GOTO  33 
TEL. J=IEL ( 1* 

TrL <J=I-LL J-1+NEL0 < 1» 

GO  TO  32 

T1  TrL'_J=IEL  (  IS  ♦  1 ) 

IEL  J  J  =  IELl  J-l+NELS  (ISHI 
12  031=0.0 
OT Jz c.0 
ooj  =  T  <IrLLJ) 

IEJ= IFLLJ 
GO  in  34 
33  103=0 

3-4  so  ->o  IEI=IELLl »  TELIJI 
I c  (TSC.FO.0)  33  TO  3« 

OTI  =  ooi 

br»i  =  rTT*T(iEi» 

rLE ye nts  of  adjacent  segmcnts 

’5  0TJ  =  A«IM1  (D3I,33J»  -  AMflXl  <oTT,r)Tj) 

KF0NrHT'O=IEI 

<0  3'l  o  j  (  TK)  =IFJ 
9<I<  )=*I J 

OOM3T  »r<)=FLITEI)/<?.0*9IJ) 


OfTUP 


Uv) 


OONU<I<)=EHIrJW(?.G*3TJ> 

CON  tag  (TO  =  0.0 
r’r  (  ki  =  t 
<”  j<  to  =1 

r<=T<*i 

!<•  (  n'>j.GT.T'I)  Gn  Tn  33 
I f  < If J.NC.IELUJ)  GP  Tn  35 
TS1=  0 
GP  TP  33 
T5  lEJ^TFJH 

PT jrpoj 

mj=  OTJ»T(  IEJ) 

Tc  <  ptj.lt  . “>q i »  go  TP  35 

ELEMENTS  of  SANE  SFGMENT 

35  IF  ( IFI. EO.IEL JT >  GP  to  40 
<rONPI(TK)=IEI 
KP3ND  »(IKf=IEI*’l 
3TJ=FL(TrI> 

3(I<  >  -91 J 

CONTI (I<I=T(tcii /(?,0*9TJ) 

C  p  N 1 J <I<)=T(iri»l)/»2.a*9TJ> 

coNOftG(io=o.o 
KC3I (  IKI  =3 
KC3  J  (  T<) =2 
I<=T  V  *1 
T3  CONTINUE 
40  CPN  T  I  *JUC 

-IPONHC  =  I<-1 

SO  UP  MflTFRIftL  PROPERTIES  TftOLES  FOP  600^0  MATERIALS 


NML  =  N“ATT»1 

IF  ( NML • OT . NMA  T |  GO  TO  2  5C 

00  350  NM=NML,NMAT 

RP40  (5,1)  MCOR 

IF  (MPOp.pp.O)  GO  TO  115 

TP  <  TNOUT.EQ. 0)  GO  to  113 

TF(HC09.GT  .  N  m  fl  T  T  )  GO  TP  112 

W p I  T  F  ( 6  *  q  35)  NN , PPPOD( MOOR)  ,  c RCOR (  MC  OR  f-4 ) 

GO  TO  11 3 

112  W3I  T  F  <5  ,Q351)  NM.MPOP 

113  RE  4  0  (5,1)  <CH 

TF  (KPH.FO.l)  GO  TO  115 
NOT ( MM) =NAT (WCO3) 

M  T  J  -  M  A  T  (  N  M  I 

on  114  NT=1,NTJ 

Tfll-(NT,UH)=TfiLr(NT,MCOR) 

FMT(fiT,NM>=FWr(Wr,«PPR) 
on  114  nw=i,nwl 

114  ALFA T (NT ,NW,N<) =ALFAT  (NT ,MW,MCOR) 

GO  TO  1?0 

115  3CA0  (5,1)  N  A  T  (  N  M ) 

NT  J  =  MflT ( NO 

IF  (INOt)T.El.l)  W=ITf  (  5 , 0  35  )  PPKCH(l)  ,PRKCH(12) 

IF  (INOUT.EO.il  WRITE  (5,018)  NAT(NW) 

00  11*  NT=i,ntj 

RPftT  (5,P)  TALF(NT,NM),FMT(NT,NM> 

RE  ft  0  (5,11  NWLft 

OO  115  VW  =  1,NW'_  ft 

115  PraO  (5.21  WL ( NW ,MT) , ftLF (NW , NT  ) 

OP  1 17  NW=1 ,NW_ 

CftLL  INT17  (12,  WflVEL(NW) , MWL ft  ,  WL ( 1 , NT > , ftLF (1 , NT)  , Al FA T ( N T , NW , NN )) 
117  CONTTN'JE 

IF  (  TNOUT.FO.O)  GO  Tp  ns 

W3I  TF  (5 ,3m  TftLF  (NT  ,MN)  ,ENI  (  NT,  NM)  ,NWLft  ,  (WL  (NW,NT  )  ,  ftLc  (  NW,MT)  ,  MM 
1=1, N  WL  A I 
11?  continue 

120  IF  (NCOp.fq.o)  GO  TO  125 


Ir  (vpw.LT. 2)  REAP  (5,1)  <CH 
IF  (<rH.ro. 2)  GO  TO  125 
NCf  (  NO  =  Mr  t  (  MO 0 3  ) 


..J 


64 


SETUP 


<1  )(  NK'I  -  k-  L 
MTJ-  f'f'T  (  MM  ) 

00  1  ?2  'IT  =1,  MT  J 

T  r  )  \l  r><  M  T  ,  t  J  M  1  -  T  3  -)  f.  ■,  (  f,  T  t  k.  ^  ~P  I 

'P'l  'L  <UT  ,NM)  =C'1N n  <  MT 

it  c mp  (  h")  .  r c.  i  >  r.p  m  \?? 

con  iprMT,NM)rr'iMrn(vr,«T'p-'* 

1  J?  C'VM  T  T  M'JC 

Go  r  n  1 1  o 

i ?5  7"m  (r>,i)  < n ( 'i » 

Tf  M  <c  ,  It  ‘<10 T  (  M 11  > 
ur  jr  mpt  (f'x) 

TP  !  ?  '  MT=1*MTJ 

,r  5  1  IG,?)  TOONO  ( ‘IT  ,NM)  ,  r-'N'It  (  NT,  NMl  ,  CONO  1  (NT  ,UM  ) 

IT  (  M.n  ( NJM)  .n. !  I  CP(  DO  I  MT  ,NM>  =  Crpnc  (NT.NM) 

t  ;>■»  CT>M  I  1  Mil" 

I r  (  TNOIJT.  rp.  Cl  G  C  TO  1  TC 

M'-’ITf  (n.jTft)  o?<r  H  (?  I  .  PP^CH  <  1  T) 

«>I’f  PP?'!  <L  0  (  MM)  ,NCT<  NM  ;  .  (  TCOMT  I  .  t  t  C  CNO  L  (NT.NM)  ,  ',CNpn(NT, 

1 NMI  . MT:l  ,UTJ) 
i  (3  tm  mco&.  rp.a)  go  to  ms 
t c  (KrM.ro.'i*  go  To  mi 

IF  (  TpH  .  L  T  .<*  )  O'AO  (G.l)  •'CM 
T c  (  KTH.ro. U)  GO  TO  ms 
141  r><“ •> r  cm- )  =  ncptmcom 

NT.)  =  N  0  r  T  INII 
no  «  up  \,T  =  1,NT'J 

r"->i  -  r 3 ° ( “i t  )Hr)5' 

m7  ?p«*r'T,NMi  -  tpm ( v t ,*c.o  ' ) 

on  rn  tsc 

14S  o-  1  !  IS  .  «  )  NC.PT  (  \IM  ) 

N T J - M CP r  cum) 

00  i  /.  '  UT  =  1  ,  »|T  ( 

14  a  ■’"51  IS.’)  Trn(  NT, MU*  .CDM(NTfNM> 

TT  I  TNOUT.TT.0t  GO  TO  ISO 

■•'"‘I’F  mp3  'll  o?KC  M  14  )  ,  nr<c-  ( 1  G) 

»o['r  (n,a-n|  \|prT(N|Mj,  (Tfo  ( ‘j  -  ,  \|M  »  ,  p  dm  (  NT  ,  NM  )  ,  NT  -  i  ,  kj  T'J) 

mi  if  t  urn-' ,mt.  ot  go  re  i-m 
t r  iKTn.rn.ci  go  to  mi 
it  it,:h.lt.m  ?-m  <r'.i)  Kpw 
t  r  (KTH.fO.S)  GO  TO  1  r-S 
T.t  tm-  T  f M  -  )  r  TM"l  r  (  m  '  .ok'  ) 

GO  '0  1  nO 

1  T’S  .ITJ  1  (T,,  ?)  fM'LT  IMM< 

tt  ( T nou T . t^ . on  to  imp 

'•1^ I r  r  (-  ,otm  s?<r  -t's),r,p<CM(if,t 

M'irr  (r,t( ,?3)  I  itlt  (►,•») 
mO  TO  ("TO'.Tn.f'i  on  T r  m3 
tp  KrH.i".i'l  GO  tt*  m» 

I"  (K  PH.  IT.  St  pap  n;  ,;i 
t t  ( y 'h . r  0 . G I  Go  TO  1 SO 
1  r,  1  H  TP  I  ' '  M )  HOF  I  V  -  0  3  \ 

r,r  t  n  <  'C 

m.G  o  -  4  0  (S.?>  H  0  p  I  N  M  I 

IT  (TMOJT.FG.ct  00  ’P  1 70 
’•’■'I'T  »?  KG  M  (S  )  ,  m  (  «  71 

rp  ls,:’4|  U0T(*p) 

1  t.  TT  I  ■••'OF.  ro.  .3  |  GO  TO  ?lr. 

T  T  I  Tim  ,FO  ,0  j  GO  TO  cn 
T  T  I Tfu , |  T  ,ij )  » T  « r>  (0,1)  <  *  u 

tt  (ttm.io.iu  or  n  21 r 

211  ''M  T  -  (  *|»i )  -  JHOM(M'O-) 

G't  *o  ??(. 

’t*  n  m.p)  f  ho  m  u"  t 

’t  »rt)pn,n.ji  g p  t(.  >i  / 

ii  pp  is,  17s t  sii(Cu(,ii,o:'<r,M i??* 

’I’T!"  (>.  ,-70)  p  ( 0  'M  ■  • ,  | 

2V  pp  at  (r,,j)  I'm 

T _  mrM.FG.G)  G  l  tp  ’2G 
W  T  •  ~  ( <-  ,  ?  U  t  ‘M 

op.  rn,'.’.T  (  f  nnoi  r  p  •••  t  i  c  c.'oprj.y  OHeijp.Tc  -mi  tooftl*  gofctfipd  m  mat 
i  r  t  i  ’  l  :  ’  t 


•;rTUP 


65 


orT  °“TNF  FHS5I  VI  T  Tr$  nr  AonrD  NATFRIALS  IF  PF-)ulprn 

2’G  IF  < F“I ( 1 ,MM) .GT . 0  .  ci  $0  *0  ?C0 

MT  J-  ‘'AT  (  T"M 

no  "ur  ntm.nt,j 

M‘">-  ’ 

FMI'>f  =  3.  o 

’111  =  ?r'9  0  0  0 .0/TftL  f  (  NT  ,MM) 
no  *  t.  n:’, l  13 
F  Ms  M 

7FT*.  =  0.  l*rN 

■1 L  ?  -  W 1. 1  /  F  ‘I 

T  c  (  wi  ?  .  L  T  .wi  I  1  ,  MT  )  .03.  WLO.GT  .  ML  I  NHL  A  ,NTII  r,n  TO  ?2<K 
CAL'  TNT17(4,NL;,,MWLft.NL<l,NT>,ALF<l,NT),AL) 
f'‘T?^-r-lT<',c»AL*7rTA**'»/(rvn(7‘-Ta»-l.C» 
r.n  Tp  <m  o 

??«  T r  4MF.GT.CI  00  Tp 
Nra=  1 

N7!T  r  (6,??0)  NN.X'T 

FORMAT  <  14H3r0?  MAT^KIAL  .  I  2*  1  6H  IT  T  F*i  PI  R  A  TO  RF  ,T  2, 

17-H,  TN  t  F  F,PA  T  x  on  TOP  r M  T  Or  T  V  ITY  RFO'JIRFS  A  noOR’T  I VI  T  If"  OUTSTOE  OF 

?  p4\ir.r  ppovt nr ■>, /4UH  ao^o^pt ivttie«:  taken  as  o  *THfPr  unavailable) 

7T1  Ofvjr  T  Myr 

?43  mNTTM'J- 
?r.T  '•p‘1'  tmup 
?(  :  TFJ=MrL 

00  07C  TF-l.IFJ 

T'4r  T  MAT  (  TF  ) 

?7Z  WT  (  :  r»  =ft  ( If)  *pmn  (  NM) 

Of  T  r  =1 »  TF  J 

Ta0  T‘iELT<Ir>  ~  C 

OF  T  UF  0  I S  T  A  N0r  S  FROM  L  F  APING  FOGF  FOP  CONVECTION 
TO J: NKEG 

no  ■« f  f  10=1,  toj 
T  F  - 1  r  1  (  t  c  1 
VOO'I  (  TF)  -  v  (  TF) 

I~1  =  IF  ♦  1 
JrJ  -  IF  ♦  NFL  0 ( I n  )  -  1 

TF  (  TFU.LT.IFH  GO  to  4  « 0 
nn  '.<-0  i-Tri,T-:j 
463  X~0N  (  T>  =  XC  ON  (  Tr  ) 

4H3  '■PN'TWllr 

Pi'Al  <5,?)  V  F  l  1  5Pi  TAH.Hn  nvo 
T  F  |  1  NOUT  .TO.  G  >  no  TO  60(1 
c  T  T  F  (h.oF?)  l/zL.PO.TAN 
Tp<  (rONO.GT.  0.  C)  WRTTF<6 ,053)  FTONO 
6nj  Pr  A  ■>  f  6,  1  I  ►  OPT 

IF<  T  MCUT  .F.T.G  )  rt  P T  Tf  (  6 »  J  54 )  NOPT 
prAi«^,?)  nn«flx 

TF(  I  mdiit  .r.T.O)  <*PITF  (F,  nor,»  nn*AX 
IFINOTT.rQ.3)  GO  T0  6  6n 

PrAT»S,2>  TOUT 

IFJT  MOUT.GT.  0)  rfRITF <6,056)  TO'JT 
no  »  n  f.  j  r 

r,r.g  NFLUX 

I F  T  T  M  OUT  • GT . 0 )  4RTTf(F,,307) 
no  n?t  t-i.nflux 

Rr  A  1  *  r-  «  ?  I  TFLUX  <  I  (  ,  FLUX  T  (  T) 

IF(  t  HOIIT  .  O.T.  0*  Tr  |  I ,  TPLUX  «  I  *  ,  FL  JXT  «  T» 

6°0  PON7  T  MU  f 
60  3  ‘it  n a  rc;  =  i 

PrA:H6,l)  MH L  T,NTSpF0 

lr< vpl! . ro. j)  .o  to  top 

TFf  1  MOUT  ,GT.  P)  iT^TTF  (6,0501  NWLT.NTSPFC 

Tr<  NT^o-r  .{  T,  Y)  r,0  T0  6K  0 

“'  A  A  IF,?  )  (Tc.n;3  <  T  )  ,  T  -l  .MTSPrr  ) 

T  F  (  J  ‘i  OUT  ,  G  T  .  3)  /4RI  TF  ( 6, 06f )  <TS0£C<T)»I  =  1*NT*1PES) 

663  IFCTMfmr.GT.Ot  4PTTf«6,RFi» 
no  ,kc  i  =  ‘  ,ntsp-:o 
nn  j«g  j-i.nnlt 

SETUP 

6b 


'7ri!<5.7)  WA  7E(J,I  )  ,0'»  J,T> 

TF'TMnijJ  ,GT.  1)  JPITr(fc,  ?=,-»>  I,J,wsvf(J,I)  ,  03  (  J*  I ) 

«•*}  r«-.jr  Tj-'jr 

M'i  ’-aifc,  -.*)  nrl.  TT^,  TST0-'i  r'"'T‘,'T 

TF<  1  sail  .  r.T.  0)  oelTTN,  TSTOO,  point 

->’r,r- 1.  c 

ffc,crT) 

rr  i  =  ” rL 

n-(  -  4.  T  "=1,150 
-  I'-EGIIF) 

Tr  <  tf,ne. tel  its)  )  r,o  to  «■  jo 

TL  =  ! 

r.o  ’o  i<?: 

•Ml  IF  f  T-.U  jrp)  .F1.lMnT(T£-n  .ANO.T(IE)  .EO  .TCTE-U  *  GO  TO  3?0 
1 1  *  TL  +  1 

N-'I'r  (h,'02»  IE  *IS,TL,tMAT(TE)  ,  T  ( T  F  ) 

<<*:•  com'  ’  .vij" 

FP  >rrjj\ 

r  5  FORMAT  (IV, 14,  11X ,T4, 11X,IU,12Y,I4, 11X,F7.4) 

orj  cnv<M  (?FM1  fjhharv  r,F  T  h  £  r  m  a  l  MOOEL /8  MOELEMEMI ,  7V  ,  7H3  EG  NFNT  ,  l  X  , 

1  jv,  *  mi  AY  E  ^  ,1  av,  3  HV  iTf  ..-x  4L,  TX  ,^HTHrr<MFSS/ 

?7H  JII-SFF,  A  V  ,f,HNUM3E°  ,QX  ,  6HM'JMPFp,H  X  .4HC0DE,  11  X, 

74M  t  ,i,//  ) 

Ofi,  FOMA-  <1HH  X  COORDINATE  , IN/ (IX, 5tE15. 7M 
or,  rn'J'if  t  ne-H  Y  OOOPOTA  ATE  ,TN/  <1X,5E15,  7)  ) 

-313  F  0  i  A  t  (  2GH0N'mFP  OF  SEGNENTS  ,151 

(’ll  F-)J'M’  <  1  r’R1ScGNFN  I  NUOOF-3,  xxypp.h  MUMOfR  OF  LAYERS  IN  S  E  G  ME  NT  / 

1  (  7  C  T  S  )  ) 

PIT  FORMAT  (l*H0L4Y-5  NIIMOCP  ,5X,29H  THICKNESS  OF  LAYFR  (  TL  4  Y  )  ,  IM ,  5  X  , 
1>’^  o  00r  PFFINING  MATERIAL  OF  L  A  Y FR  ,  5  X ,  7“  M  NUM3ER  Oc  ELEMENTS  IN  L 
PAY"'  /  t>3X,T?,i?x,E15,7,?'ix.li5,  ?0X,T15M 
PIT  '■nr.'  V  {U~  M  pp  t  r  mf  rc  A  T  ij  RE  S  F0»  A  0  S  OR  P  TI  V  I  T  Y  ,  t  FI 

pip  ^ r a t  (2c,h  tcmpf'p 4T(Jpp  (tulfi  ,nrG.R,ri5.  7,<;x,ithfmi5SHItty  iemd, 
1c15.7/2EW  M'J-PiIR  ?r  WOVE  LENGTHS, I5/25H  Ml  VF  LCN  GTH  f  WL)  ,  NINONS  , 

7  1 1  -  M  ra  S  0  R  °  T I  YTTY  (ALC)/(5X,E15.7,17X,E15.7>> 
pro  f.OJlf.T  (5H  E0Or,T5,44M  Vp  MUMPER  OF  Tc“DCC>A  TURES  FC°  Z  ON  0UCT IV I T  Y 
1,TS/?M  TFMPfPATu°F  (TOOMPI  ,PFG.R,f<,  TGH  LENGTHWISE  ICONOL)  ,PT'J/TN 
2“SFT“r'EG.  -  ,5X,I3H  orpTMUTSE  (  CONDO ),  MTU/ TN-SEC-OEG.  ? ,  5X/ 

7  fr'  <  .  r  1  f  .  7 , 2 IX  ,  El  5.  7*  2  6X  ,  El  5 . 71  I 
P?P  for  '  A  ’  f  u  ih  NU*7p->  OF  T r  wr E  A  T  URE  S  P0R  SPECTEI3  HEAT, 15/ 

7  p.,H  TI  HOfP/iT'jij-:  (  T  cpl  ,0rs.  ?,  1  ?X  ,  ?TH  SPECIFIC  H- AT  {  C.P  M)  ,  3  TU  /L"<- OF  G 

<  ,o/ (  r<v,F1ci7t2x)(,rts.  n  I 

47p  FORMAT  (  34H  MELTING  T  F**d  Er  A  T  Ua  F  (  THE  L  T  >  ,  D  CG .  D  ,  ?X  ,  F 1 5 . 71 
7  >4  POM"  I  (?*H  HFir  of  FUST  O'1  (  HOF>  ,  0TU/L0,2X  ,  E  1  5.  7  ) 

473  pom.U  OFNS I T V  (FMOM)  ,LOP-IN**3,  ?X,F15. 7//) 

P7S  t'07'JAT  (70H3NFW  MATr;TAL  N')W3-  R,  I7/ 

1  sew  MATEPTll  to  which  vjEW  MATFPIAL  CORRESPONDS  MCO’!, 

r  ?uiai 

p  T  f  i  rn?iJT  (  ?  0  H  0  N  E  M  M  A  T  r=  TA  L  NIMppp,  T3/ 

1  r!H  Wfli-Diai  TO  Hit  Tow  NEW  MATERIAL  COR’rSPONOS  (NrORI,I3» 

p<c,  FOR-M.r  (  p-MawATEPrAL  F4>(0DF5>tY  IQ  7F  CHANGED  (<CH»,2A10) 

•04  j  TnriJ;  (l7HflTTMF  tncp“'H1Tt0N/ 

4  4  ?w  ‘  INTEGRATION  TT  ME  INTERVAL  C  0  EL  T  I'M  ,  SEC  ,  2X,  E  15  .  f  f 

GR'TH  STOP  tthx  MStOP)  ,  E  F  F  ,  7  X  ,  E  15 . 7  /  7  GH  NJWR1ER  Op  TINE  TNTE°V 

1  A  L  S  »‘»'Tf'T),2V,Fl4.7) 

phi  forma  i  (yhiS'-ctio"  t?,35h  contains  more  elfyents  than  allowed! 

007  FORMAT  (  '4M0  (/."LOCI  TV  OF  AT  7  STREAM,  FT/SrC  =  FIS.  ft/ 

PH  A  1  o  I  E  N  T  (POr«H>  pf  E  3r  J 5  E  ,  PST  =  F15.5/ 

7  T  17  A-MTFHT  (ROTP  TFMF-RATURF,  peg.  F  =  F 1 5.  p ) 

-ft  FOR"  -  (F  v,H0CO"VrC  T  I  »*F  It  ■  A  r  TP."NSrFP  COcFFICTENT  ( D  TU /I  N*  *  ? /Ffc/ OR 
H  -  ‘:F.-I 

or',  F-pvt.n-jpjunrHFRHAI  FLUX  OcTTON  13) 

3'IS  FOR'  V  (  irut  Ml  XT  v  u  M  FLIIX  (  P  T  IJ/ 1  N»  *  ? /S  eC  )  =  E15.E) 

s  r  o  ?  “  A  T  (  7  o  h  -  i_  U  v  CUTOFF  AT  T  (SFp)  =  F15.5) 

OF  7  ^OR  M  ’  ( 1  w:  ,1  4X  ,  4WT  I  V|-  ,  i  3X,  UHFL'JX) 

'•CH  Fn-  *  A  M  t  7.  ;>c  is.F  ) 

P-3  F03"Ai  f<,  rup<;  or  C  T  ?  A  L  r!  IS  T  7T  OU  T  I  ON  OF  THFSMAL  SDJRCc/ 

1  7*h  UU'TPR  -)r  NAME  L  -  NGthc  =  I  ?/25H  NU'MER  OF  T  IN*  s  =  I?> 
rir-jrijf.iiu  TTxr-  -  3017,.*,) 

pt  1  F'-ov  T  r  1  w  •  ,3  Y,  IH  J  ,  3  V,  in  J  *  '  X,  lr  HNAVFL  ENGTH,  *x,  4-<rLUX» 

p.-.o  FOP'tA  '  (  >  14 ,4X  ,  ?F  1  f  ,  f  ) 

-■p 
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PrTUD 


SUfOUTIME  S»nr  <  IEIRST.VTSPEC, TSPFC.NWLT, WAVE, *3R, NWL, 

i  wa vfl , waver, t, opp* 

OFVPhINE  SPECIAL  01 STR  T^'J  T  ION  OF  SOURCE  AS  A  -UNCTION  OF 
14  WAVELENGTHS  OETWFEN  c.<?  mICR1NS  and  s.s  microns. 

IFI^r  -  0,  FT5ST  PASS. 

1,  ALL  OTHER  FAS^EF. 

NTS5FC  -  NUMBER  0^  TTMFS  THE  SPECTRAL  0  TSTRI  TUT  TON  IS  3  5  E  C  T  ETrp 
TFJ-f  -  TIMES  AT  W  H  T  0  H  T Hr  SPATIAL  0ISTRI3UTI0N  IS  SPECIFIFO. 
MHLT  -  MUw,ITrR  OF  INPUT  TABLE  ENTRIES  (JFfiming  0 1 S TR  I bi jt  1 0 N.  TF 
0,  USE  -3JTLT-IN  OTST^IOUTION  FOp  6000  WATT  TUNGSTEN  LAMP. 
WAV  -  INPUT  WAVE!  ENGTHS  . 

-17  .  input  INTESITIFS. 

NWL  -  NJMSFR  Or  STANPARO  WAVELENGTHS  -  PROBABLY  14. 

W  A  V"  L  -  14  STANDARD  RAVELr  NGTHS. 

WAV-0  -  IS  HAVrLFNGTHS  OFFJMTNG  RAND  WIDTHS  OF  HAVEL. 

T  -  Tint  . 

drp  -  relative  intfnftty  fop  14  wavelengths. 

01 Hc  nf I  ON  OR (50, 3) , WA VE (SO, 3) . MAVEX( IS. 3) , Q»X (IS, *1 , 

i  VPY<T,i4>,rm(i4»,TSPErm,c<5> 

BIN- N'T  ON  HAVEL  1 14) ,W AVERT  IS) 

0  A  T  A  NWLX/1S/,  WAVEVG.  2, 0.375, 0.625, 0.875, 1.125, 1.  375,  1.625, 

1  1. a75,2. 25, 7S, 3.25, 3. 75,4.25.4.75,5.25,  30*0.0/ 

OAT  3  QRX /0 . 1, 0. C 85, C. 63 ,0. 99  5,  0.6  9,  0. 59 ,3.44,0.  31,0. 205,0 .06, 

1  0.  055, 0.  037,0.  02,  0.  r'15,''.  0,  30*0.0/ 

OATS  C/l.C,4.0,2.0,4.0,l.C/ 

IF  (  TFIRST.r.T.O)  GO  rn  63C 
TF  (NWLI.GT.O)  SO  TO  2C0 
USE  4  KW  TUNGST-N  LAMP. 

NTS PrC  =  1 

NWL  I  =  NWL  y 

00  100  J=l,NTFOEC 

00  1  0C  I  =  1 , N  WL I 

WAV  (I.J)  =  WAVEX(I.J) 

ioo  qp(  i ,  j)  =  ooxa,  j) 

INTEGRATE,  US  TNG  STMPSON'S  PULE,  OVER  EACH  -1ANOWIDTH. 

200  00  50C  J  =  1 »N  TS3E  C 
OTOT  =  0.0 
00  '4  Oil  1=1, NWL 
XI  =  WAVEB(I) 

X5  :  WAVF'MT*1> 
ncLX  =  (X5-XH/4. 

FP  =  0.0 
OP  <0C  <=1,5 

CALL  TNT  17 (1,  XI,  NWL I, WAVE(1, J) ,DR (1, J) ,F) 

XT  =  <1  ♦  PELX 

3C0  FF  =  FE  ♦  E*CC<) 

0°P< I )  r  FF* (X5-X1) /12. 

400  OTOT  =  OTOT  ♦  ORR(I) 

OP  ^OC  1=1, NWL 
SCO  1D7  X  ( J,  T  )  =  ORRm/OTOT 
00  ’  5  C  T  -  1  ,NWL 
553  QPR(T)  =  OPRX(l,T) 

PrT  JPN 

Spr;TPAL  DISTRIBUTION  VARIES  WITH  TIME. 

600  TF  < NTSPrr.EQ.l)  RETURN 

Tc  (  T.GT.TSPEC(')TSPFC))  GO  TO  800 
00  700  1=1, NWL 

CALL  inti  7  (3,T,mtfpeC  ,  TS  °E  3  .  OR  RX  1 1 ,  I  )  .H^III  ) 

700  CONTINUE 
PET  JPN 

8oo  no  goo  t=i,nwl 

POO  ORR(M  -  OPRX  <  NT  SPEC  ,  I ) 

p<-r  jpn 

FNO 


SPFCT 
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GUI? OUT  INF  TIM 
?'Ji  IN  TM°UT  OUa. 
dtn-msion  M(40) 

M- IT f  (5,1330) 

or^TMp  5 

J  =  1 

r3  a D  (i,i)  (  ami  ) ,  i=i  ,40  ) 

Tr  ( F  OF (  1 )  )  400,000 
200  WJI  r  r  (5,2030)  J,  (  4A(T)  ,  1  =  1,40) 
W*ITF  (5,1)  (AM I)  ,1=1,43) 

J  =  J  ♦  1 
GO  TO  13 C 

4^3  M?I T F  (5,3003) 

3CKTM3  5 

j°n 

i  ( 4 n a? ) 

looo  p c ? u a t  < ihi, /  /, icx ,  15ht  5>  a  p  _  M 

1  15M0***  INPJT  04Ta/t 

2000  F03-1AT  (  IV  ,I4,?<  ,U0A2) 

3  0  3  0  TO?iaT  (2-»H0***  r  NO  OF  TNPJT  04Ta) 
Fhn 


sm’n'jTT*T  tste3 

CONDON  A  ( 75)  ,ALr(i&,lC)  ,0(122)  ,CKD<75)  ,  CKDM < 0 ) , CKL ( 75) , CKLM (a)  , 

1  3i«cr,  (122  ) , CON3  ](li,3),  COND  Ul  2?),  C0N0J(122),CP(75),:r>iAT(a), 

2  •»-  LTIM,pTTM(75)  ,  DTT-E,  r>rMIM,EL(75l  , FLUXT < 1 5 ) , HCONO , ID < ? 0 ) . 

7  I  L  (3)  ,TEL  AV(=>,  7)  ,  TFnD,TMAT(  75)  ,IMf  LMX5)  ,  I  NO JT , I  PLOT ,  I SEG ( 75 )  , 

4  I  5  TCP, JEL(  3) ,<CH, KCON  01 (122) ,<CONDJ (122) ,KCPI ( 122) ,<C®J (12?) , 

5  <7  A T  n ,5) ,<sjp, L AYN< 75 ) , MCOP  »  NCONDC  ,  N  FL ,  NF  LL  ( 3  «  5) , MELS ( 7) , 

5  M  7  L IJX  ,  NL  A  v  (  7)  ,  N *  A  T ,  no PT ,  NPR  T  NT ,  NSE  G  ,  N  TS^C  ,  N*'_  A  ,  N WLI .  PO  , 

7  35  TNT  ,0(75)  » QA  D  (  7  )  ,OCONn(75)  ,QC0MV  (  ?)  ,00  M  A  X  ,  D7 LU  ,  O7  (5  C  ,  7)  , 

rt  D?A0,0PFcm2,3)  ,0RP(14  >,~FGL,T(75> ,TAM,TCUT,TE7(3)  ,TE7P(75) 

0  07  7  ON  T  FLUX  (15)  ,TFF,TTS1C«TLA7  (  T  ,  5 )  ,  T  o  (  7  5  )  ,  T<?,T3PFC  (  7)  ,  TST0°, 

1  V7L,KAVf(5C,7)  , W AVFL( 14) ,WL  < 10, IP) , WT (75) , X(75>  ,XCON(75 ) , 

?  X'FG(4)  ,XT3,r(75)  ,V$EG(4) 

0  017  0-  /9LPCK/  ALE  AM  io  ,14,9)  ,  CONDL  (10,9)  ,CPM(10,9>  ,EMI(10,9) , 

1  7  D  c  (9>  ,KLD(Q),NAM9),NCpT(9>  ,  NC  T ( 9 )  ,NMATT,NWL»PI,DHCM(9)» 

2  TUF(1C,9)  ,TCDN3  (10 ,9)  .TC®  (1C,9),  T-FLT(9> ,  WAVE8(15) 


T  2=1  .  C 
TUsNFl 

00  5  IF  =  1  ,  I'J 
N7=I MAT< irj 

Tl  =  ilAXl  (T1,TCDMD  (  1,NM)  ,TrP(  1,  N**)  ) 

5  T z- a  MAX  1 ( TJ, tiel t < N-)  ) 

no  ir  ip=i,iu 
10  0TT7(  Ir)=lOO.C 
OO  ' 0  NT =1,10 

TrL  =  Tl  +  (F|.OAT  (NT-1)  79 .3)  *(T?-T1) 

00  ^0  N1  =  1,N**AT 

TF  ( tel . lt  .TDONO ()  , N- ) .op.TFL .LT . TCP ( 1 , NNI )  GO  TO  20 
IF  (  TrL  .  r,T  .TMELT  (-7)  )  r,o  =0  15 

CALI  TNT  17  < 13,r FL ,NCT(MM) .TCOND ( 1, NM) ,CONOL ( 1, MM) ,CKLM( NM) t 
IF  (<L0(NM).FQ.2)  GO  TO  12 
rxo- ( NM) =  CKLM  <N7 ) 

GO  TO  17 

12  CALL  TNT  1 7  (5,TTL,NCT (NM),T00ND(1,NM),C0NDD(1»N7)  ,CKDM(NM)) 

13  CALL  INT  1  7  (6,TEL.  NCPT(NM)  ,  TCP  (1,NM)  ,roN(l,NN)  ,C0f7AT  (NN)  ) 

GO  TO  23 

15  T=N'T(NM) 

CKLM'NM) =OONCL(T,NM) 

IF  (ELO(NM).FQ.?)  GO  TO  1? 

CKO  7  (-7)  =0  KL  7  (  N7  ) 

GO  TO  1  a 

17  CKD7 ( N  M ) =C ON 30(1 ,NM) 

1°  T=N~PT(NM) 

C3M\T (NM ) =  OPM ( I , N«) 

30  ''ONTJA'IJF 

DO  25  Ir=l,TU 
’5  Q(  TF  )  =0 . 0 
TPC  .)  =  NC  ONDC 
DO  50  icc=i,iccj 
T  =  <~  ONDT  (TCC) 
j=<co-nj (tcc* 

MM= I  MAT ( I ) 

IF  (KCPI(  ICDl.FQ.?)  GO  TO  70 
0  X I  =  C  KL  7  ( N  M ) 

GO  TO  32 
1C  CKT  =  C  <PM ( NM) 

72  N«=TMAT(J) 

TF  <KCPJ<  ICC). ED. 2)  GO  TO  4  0 
CK  J  =  CKLM (MM) 

GO  to  4 3 
40  0<J= CKOM (MM) 

4?  OKTJ  =  COwni(  ICC) /CK  T*C OND J (  ICC ) 7CK J* COND AG<  ICC) 

00  =  1 . Q/OKT  J 
0(  I)=0<  X ) »  DC 

go  3(  j)=o(  j)*or 

DO  50  IF=l,IU 

IF  (0(  IF).E0.C.0»  GO  TO  50 

N7  = I  MAT (  TE) 

DTT7(  IF (  =  AMI N 1 ( 0TI M(  IF),WM  IE)*C',M4T(N7)'*0.<,73(  IF)) 

50  CONTINUF 

73  C0NrINUE 


Tr  TEP 
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■'-'I  Tf  <0,110* 

11’  **>->•*  t,r  (  ?<SH0  T  TM:  "'Tu-M/fL  I  '.""O3’*/*  TI  ON  * 

T’*  =  'T1. 

'•!’T'r  < ^ , i  ?- >  nTi«<  Tr* ,  t r = 1 1 t t j » 

1?1  =-n=*  v*  1»  r  <  INTERVAL  RcOtJIREO  8Y  EACH 

1  /  <  1  *  1  U  r  1  ?  .  >  *  * 

■?T^T  *i  =  TTImCD 
1  =  1 

T r  (tu.co.i)  r-n  to  it«> 

n  :*c  l?  =  2,l'i 

Tr  (HTHI  TF  *  .  .  "1  MI  N*  r,c  to  13? 

V‘‘1T  V  r  r'TTl(T“) 

T  =  I  ' 

tv  r  t  \i(j- 
i7-  ,j  =  rTf  ct.iuo* 

1 1* :>  C-.5HT  <??W]V(T'JT  MU’*  TI<F  INTERVAL  IS  FIT. 5, 

1  •> •  h  ^-roNOS  FO?  FlrMF'lT  IT* 

1-0  ~  M T  I  ►  U  ~ 

-i-i  T  J  V-  njM  TM 

T!1**’  =3FL  txm 

v'T  ,I?M 
r*’1 


ELE**ENT  / 


TFTE  o 
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SDRVPUTTNF  X  Hc  4  T  (  KF) 

CALCULATES  EXtERNAl  RAOIATTON  ANP  CONVECTION  HEAT  FLOWS 
FP?  ^aCH  ELEMENT , 

<R  :  C,  SfTJP  ONLY. 

1,  CONVECTION  A NO  R  EC- A  P I  A  T I  O  N  A °E  CONSIDERED. 

2,  CONVECTION,  EERAOTATION,  AMO  EXTERNAL  RADIATION  ARE  USED. 

COHN  on  A  (75)  ,  A  L  r  *10,10)  .0(122)  ,  CKO  (  75  )  ,  COM*  9)  .  OKU  79 )  ,  CKLM  H). 

1  SO  ‘'DAS  (122)  ,  CONDO <10  »al  ,  CONDI  <12?  )  ,CONO  J  ( 1  22)  ,  CP(  751  .  CPMA  HR)  , 
3rLTIM,DTIM(79)  ,  OTTmF-DTMIN.ELI  7<?)  .FLUXT  (19).-*C0N0,  10*20)  , 
Irini,TFLAY(=.,T>,IEN0,IMAT<75>  , I  MELT  ( 7^  )  ,  I  NO  J  T  ,  I  PLOT  ,  I  SEC,  I  75  >  , 
TsTOP,JEL(3),<Cn.KC0NDT(12?) , KCONO J ( 1 22 ) , KCPT ( i 22 > , KC ° J ( 1 2? I , 
<nat ( ^ ,5), KSJR.  ,LAyM79 ) ,  MOOR , NCONOC »NEL»NELL(3.5) .NELSC3I, 

N- IIIX.NLAY  (3»  ,  MM  AT  , U0° T , MFE I  NT , NSF G ♦ N  TSPEC , NHL A , NWLI , PO , 

’  ?  I  NT  ,0*79>,"}AR(  3)  ,  DCOND  (  7  5  >  ,9C0NV  <  3 )  ,00H  A  x  ,  J- LU  ,  QP  (5  0 , 3>  , 

AC,  3fER(2,3  »  ,oqo(  14)  ,?EGL,  YUS)  ,T  A* , T CUT , Tpc  (j>  ,  Tf*<o<  75) 
C0H"OM  TEL'JX  <151  ,TEc,TIHF,TLAY  (3,5)  ,to(  751  , TR  , T5  PEC  <  3)  ,  TS  TOP, 
t  V“L,HAVr(50.3»  , HAVEL <141  . WL < 1 0 , IP ) , WT ( 79  ) , X  (75 > , XCON* 79 ) t 

2  X  CF',14  )  ,XTR,r  (  75  I  .YCf-;  « t, ) 

''ON*  ON  /  91  0C</  AlpAt(10,14,9),  CONDL  *10.9),  C°M  *13,9)  ,  EMI  *10,9)  , 

1  MIC *9)  ,  ^LO  *  o )  , N  A  T  (° ) , NC^T  *  9 )  ,NCT<9>  , NM A T T , NHL , PI , RHOM *  9 ) , 

2  TA  Lf*10,9)  ,t:dnD  *1C  ,9)  ,TC.o  UC  ,9)  ,T*FLT*9)  .HAVER  (19) 

TE’J  =  MEL 

DO  =i  IE  =  1,I*U 

Dr  3N  o  ( If; )  =  c. 

9  0*Tf )  =  0. 

no  7  ts=i,nseg 

D  r  D  N  V  IIS)  =  3. 

0S3ITEI  -  0. 

QP"  3  *  1, 1 S)  =  C. 

3  DRE?»?,TS)  =  0. 

IE  (K-o-1)  10, 33,  ?G 

10  CAL.  SP-Ct  *0, NTSppC.TSOEC.NHLI, HAVE, DR  .NHL, HAVEL, HAVE3, TIME, ORR) 
CALL  FLUX  (0, NOPT, NELNX , TFLUX, PLUXT, OOMAX.TCUT.TIME.ORAD) 

CALL  CONVFC  *3,30,TAM,VcL,tEMP,El,yC0N,HC0N0,TE,TP, DC  ON ) 

RF  T  J  p  N 

?r  r  ALL  SPECT  fl,NTSPEC,TCPEC,NHL  I, wav»E, OR  , NHL,  HAVEL, HAVER, TIME, QRR) 
CALL  PLJX  *l,NO°T .NELUX .TPLUX, FLU  XT .ODMAX, TCUT , T I  ME , ORA  j) 

33  ISJ  =  NSEC 

oo  on  IS=1,ISJ 
Tr  =  JEL  *  I  9) 

IE  (TE.FO.O)  GO  TO  90 
INN: p  SUCP  ACE  . 

N«  =  IMAT(IE» 

TFl  =  TEMp  *  T  F) 

IE  * TEL.LT .TMELT (MM) )  GO  TO  42 
I  =  MflT(NM) 

EMI  SC  =  EMT* I » N  N  ) 

GP  IP  43 

4?  CALL  tNT 17  *14,TEL,NAT<NM» ,TALP*1,NM) ,EMI*1,NM),EMISS) 

43  QRE?*?,IS)  =  -* 1 . 0E-1 4/3  .C ) *rL  * IE)*£MlSS*TEL**4 
3*  I E  )  =  O  *  IE )  ♦  0PEP*?,IS) 

90  IE  =  IEL  *TS) 

IE  IIP. ED. 0)  GO  TO  oo 

calculate  oomvestivf  wrar  floh 

CALL  CONVFC  ll,PO,TAM,VPL,TcMP,EL,XCON,HCONO,IE. TR.DCON) 

3C0NVIIS)  ^  QC 5v 
0  €  T  ~  )  =  DUE)  »  DCOM 

CALCULATE  PAOIATION  EPOM  SURFACE  ELEMENTS 
NM  =  IMA  THE) 

TE.rTTMP (If) 

IF  (TEL. LT. TMELT (NM))  GO  TO  62 
T  =  NA  T  (NM ) 

FMT?S=EmI(T,NM) 

GO  TO  63 
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j 


1 1 

*.  3 

oa 


0*^11  1NT17  (  9  ,  Tf  L  ,  NAT  (Mm  *  ,t  ft(.c  »  ,  EMI  ( 1  ,N“* ,  =  MTS'-, » 

r\3’  *?, is*  =  -(  i  .  o  t-i  u/-> .:  i  *  ft.  ( in* -Kr<;<;»TCL  **4' 

KT;l  =  FfTF*  ♦  Tc-t"c(i.I^> 

rr  (K'-.if'.i)  ri  rn  -n 


ronnr  °."Net.s  a m 0  sfrTn»i^t 


►n=i 

r,n  . 

tr  < 
T  =  N" 
ai"i 

m  * 

v?L_ 

Tfl  3  * 
0(T- 
rnMX 
3r  r  J 

EA'T 


MiT< IF > 

c  iwl=i,k'hl 

t~l  .lt.tm^ut  (N't* )  an  to  e 

rn^t 

=  .*  L  r  £  T  'I.  TWL  ,NM* 

n  nr, 


NT  17  (  1C,  TH  ,  N*  Tt  N'D  ,  TfllF  (1  ,NM)  ,  3  LF4T(  1  ,H|_  ,NM*  ,t 

lc>  =  o a ^ f I s  >  ♦  Q^an’ORs ( iwl> *atPa*FL » Tr ) 

>  =  O(IF)  f 

T  NUr 

e*l 


LFa* 


XUF  A  T 


CKi  U  Vi  wi  Vi  vi  ViViO 


PRO  SPAM  APLOT  f OUT PUT, T A PE6= OUTPUT,  TAPES  =  5 IV,  TAPE40=514,  PLOT1-  J 
THT  S  f'Ror.oan  PLOTS  thp  DATA  GE  NEP  AT  E  0  9Y  T9AP-M'_  and  st  op  eo  on 
FILES  tapes.  Twe  OUTPUt  plot  TAPE  is  TAPE  39. 

DP3  SPAM  APLOT  was  WRITTEN  <3Y  <  AMAN  AVIDYNE,  SE°T.,  1973. 

COMMON  /'■TOPE/  0AT(200*S),  100(20),  00,  TO,  III,  TI«(f,00) 
COm-!ON/<ACOM/NCPW,  iriN(21,ICOUT(l>,INP<3>  ,  NOP  APS,  NOGPAf.hcpt 
DIM-  NOION  NPO  (5)  ,  XMT  (?)  ,  Ymt(2)  t  xrm  ,  YT(4)  ,IT2(7»  ,  IT3(5>  ,  IPS'*  <6 ) 
DIMENSION  IPMT(l) 

DATA  NY / 1 C /,  NY/13/ 

DATA  XT/4HTIMF,4H  (SE,4WC>  / 

DATA  YT  /  AH  TEMP ,  4  H  ,  (0,4HFG  p,4H>  / 

DATA  NTl/40/,  NT2/25/,  NTV22/ 

DATA  TT?/4HPcA<, 4H  rL  U , 4  HY ,  5 , 4 HTU/ 1 , 4HN 2 /S , 4HE0  =,4H  / 

DATA  TT3/4H°ULS,  4ME  OtJ ,  4H=  A  T  1 , 4H0N ,  ,  4HSEC  ,4H=  / 

DATA  IFHT/4HP5.3/.  IC/4/,  IP/5/ 

DATA  I°SS/1,3»5,6,4,2/ 

DATA  AL/P.C/ 

DATA  NDIMEN/200/ 

CALL  PLOT  (0.,0.,-3) 

CALL  FIRST  C 1 0  0 » 

20  IDIDT  =  f 

p-AI  (6)  (IDD(I)  ,  1  =  1,  20  )  ,  00 , TC, III, NO°,NSP 
T"  ( EOF( ®) )  900, 40 
40  wpitt  (e,f7> 

il  -  in  -  i 

50  I p.  D  T  =  IPLOT  1 

I  =  tplot 

PEAD  (6)  TIM(I),  (OAT  II,  J)  ,  J=l,Ill> 

Ic  ( EOF ( S ) )  123,70 

70  W=ITE  (6,8)  riN(T) ,  (OAMI,  J)  ,J  =  1,  TL1T 

GO  To  50 

120  NC  =  IL 1 

IPi  dt  -  IPLDT  -  1 

ELIMINATE  MELTED  REGION. 

Do  ?C0  J=1,NC 
DO  1 5  C  1  =  1, IPLDT 
N  =  T 

IF  < DATfT, J) .ED. 0. 0)  GO  TD  180 
153  CONTINUF 

np:(j)  =  tplot 
GD  TO  200 
130  N°C<J)  =  N-l 
20 0  CONTINUE 

scale  the  TA9LE. 

N=T  A  nS (NPC ( 1 ) > 

CALL  KALE  (TIN,AL,N,1,XM,OY,OELX) 

CALL  KALE  (OAT, AL , N,1 ,YM,PY,OELY) 

IF('jr.ED.l)  GO  TO  3  50 
X  MI  N  =  XM 

YMA  Y  r  YM*DX*AL 
YMT  N  =  YM 
YMAY  =  Y  M-f  D  Y*  AL 
"IOC  =  1 
DO  3  C  C  T  SC=2* NC 
N=TAPS(NPC(ISC) ) 

ML DC  =  NLOC+NOIMFN 

call  kale  (TIM,AL,N,l,Xw,DX,OPLX) 

call  KALE(DAT(NLOC) ,AL,N,1,vm,DY,DELYI 

TEST  =  X«+nx*AL 

IFtYM.LT.XMIN)  X  MI  N  =  X1' 

IFt TEST. GT.XMAXI  XMAX  =  TEST 
TEST  =  YM*9Y*AL 
IFt YM.LT .YMIN)  Y  MI N  =  YM 
IF  < TEST .GT.YMAX)  YMAX  =  test 
300  CONTTNUE 

XMT(i)  =  XMTn 
XMT<2)  r  XMAX 
YMT(i)  =  YMIN 

APLOT 
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L 


CVIJ-I.Cur 


v>T  1  ?  1  =  Ymi\  x 

Ci l -  KflLrryf'T,/\L,?,i,XM,r)y,r)eLX) 

CALL  KALF  <  YMT,  A ,  2  ,  1,  VM,  OY  ,3  EL  Y) 

THC  *.  XFS. 

350  XAX  =  O.c 

IF  (  v»*#L  t  .  o.  0.  AND. -vy /OY  .1  T.  Al  )  X®X=-YM/1Y 
CALL  <C.C«  *ax,XT,-.NX  ,Al  ,C.  0,  xm,dx,dflx,o» 

y  '  X  =  o .  c 

If  <  V'l.LT.O.O.axn.-yx/OX.LT.flL)  YAX=-XM/OX 
CALL  KAXIr (VAX,  0 . o ,vr ,nY ,AL,O0.0, YM,OY, CELY, -It 

’lr’4W  TITLFS. 

H  =  0.1 

f*UL  K4  °  Y"  (0.3,5.5,H,r.)O,0.C,\'Tl,3t 

CALL  x4CY‘M0.o,f.*,-.,tt?,C.G,NT?,3I 

call  XANUm  (9°9 . 0, 5. 3,H ,On,C . o , TFMT, T C, IP> 

IF  (VOP.F0.3)  33  T1  4  CO 

CAL.  x  A  a  Y  m  (9.9,5.i,H,  TT7,0.0,NT3,3t 

CALL  XANMM  C9q9.0,5.itH,TC,0.0»IFMTtIC,I?> 

PICT  rfiLiH  "J’Vf, 

400  ’JFJ  =  I^LOT/5 

no  sou  t=i, nc 
<  =  NHI^ENM  1-1 )  +  1 

►IP  :  V'PC(T) 

IP'  =  I?SP<T» 

CAL.  KAlTNF  (TIY  (  1  )  ,OAT  (K)  ,‘IP,  XM,OX,YM,DY  ,NF?,TPS) 
SCO  ‘v1  p T  t  r  <5,?t  T 

TfP T fa  TF  mFl T  CONDITION. 

m  =  P  .  1 

CO  -Cj  J=1,NC 

If  (  MFC  (J)  .FT.  I  ’LOT*  GO  Tn  P.CC 
T  =  A'P?(J) 

K  -  NDIMFMMJ-lt  ♦  T 
X  =  »TI  I  )-XY)  /OX 
Y  =  f  CAT  (  K) -YM)  /  PY  -  .5  *  H 

r  A  L  L  X  A  C  Yw  <X, Y. H, IOU",  ]  .  C  ,  -  7  ,  3) 

SPG  "CNriKiUC 

C'l  £CY  FOP  last  plot. 

Tr  (  f'OCPAF.PO.l)  GO  to  701 
N-’I  T  F  <5,4t  VOGFAp 
Go  TO  900 
TOO  'jPIXF  (*,,5) 

V0  NO*,?/if=>jogoaf*1 

AOIflMop  TO  MFXT  FP  A  *'F  . 


CAL.  KAYAUS  (AL*-t.5,C  .0* 

GO  rO  20 

~L G  F  THF  PLOTTING  SYSTF-l. 


►IOG  ■'  A  F  =  flOG^AP  -  1 
'l-'T'F  (5,5)  ►IDG^Ap 
CA_.  PLOTF  CNN* 

OT3r> 


7 

9 


FO?mAT  .'TATSMt'NTS. 

F  O  ? 1  A  T  |ph  CU'Y'T3,1X,7MPL0TTFQ/I 

FOVfT  IFX,7H**'‘****IT,1X,?2HC-PAPHS  PLOTTFO  .****#*•  ) 

F  O  ?  '•*  A  T  {  Fv,34H*»  ****•■  FT  :>ST  G"A  PH  PLO  TTFR.  ***•*.*»  J 

fo?v<at  (  fx,  7M»  »»*-►♦*  I  t,  iy,  3  7HGPAPHS  PIOTTFD,  SYSTEM  CLOSEn 

I 

FOP-tAT  (13H10ATA  3L  °T  Tr  0  /’ 

FORMAT  (1H  ,  F15. 6, 5  y,  r,F!  s.  6) 


APLOT 
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:X>wK)OUO^CKK)CKAKiOOOUO.#OOtKKKJUUUw)  lK>  AK  A  K)UUl  )C)Ci  i;  OU 


subroutine  rro3TfMi 

OPEN  plot  filts  FOP  OFF -LTNt  PLOTS. 

M  =  MUMPpp  OF  GRAPHS  TO  nc  DLOTTEO. 

comm  on/K  acOM/NCPW,  KIN<?)  ,  TCOUT(l)  ,  JNP<3) ,  NGRAF>  ,N0GPAP,*1CRT 
DT  M  r  MF I  ON  I(2>,<  I  *) 

DIMENSION  PLOTIMJ) 

0 A  T  3  I  /  4  H  (  ,4444 1)/,J/4U(A4) /,K/4H(  ,4H18A4,4H)  7.L/4/ 

DATA  PL0Tin/lC47?09  HILL  ,  1  OHL  EE  RH0NF  ,10H  ? 7? 1990 
N:=I4aP(M) 

XMS<-1C*M 

VMA*=11.4 

CALL  PLTI03f PLOTin, XMAX,Ymax»1 . 0) 


M3?r  =  FLAG  F0=  TTD.F  or  DLOTTTNG. 

REPUI^EC  4  T  SOME  INSTALLATIONS  FOR  PROPER  FILM  A  D  V  6  NOP 
AND  OPTS  IN  PFFET. 

0  -  ORIGIN  CONTROLLFC  oy  ijspp. 

1  =  ORIGIN  MOVE 0  UP  AND  0VF3  2.0  INCHES  FOP  USE  with 
QUICK  PLOT  ROUTINES. 

NCRT=G 

IF  (M.GT.G)  NCRT  =  1 


IF  (M.GT.C)  CALL  < AVANS  <3.C,1.0> 

I3M  7<S0  -  A  VCO. 

CALCOHp  YOOEL  890  CATMOOF  PAY  TU°F  pioTted. 

MAXtMUH  PLOT  AREA  WITHOUT  0IST0otfON|  =  7. 5X7. 5  INCHES. 

CALL  IDFRM8 ( *  AVI DYNF  ' ,  ' A5 ' , 0 , T , 39) 

N=IA °S(H) 

IP  M  xs  LESS  than  7EPO.  UPEP  is  E  X° EC  T17  D  TO  ADVANCE 
FTL*  SET  OHO  AVCO  ID  FRAME  kith 

CALL  CALFHP(C. 0,0. 0,0,2) 

IP  J SrR  CALLS  <A  V  A  Nc  ORIGIN  WILL  *5  MOVED  UP  1.05  tmchES, 
ANT  °TGHT  1.25  INCHES  -  INTENSITY  WILL  PE  CHANGED  3Y  -5. 

CHECK  FOP  OUTC<  °LOT  OP  USER  PLOTTING. 

MOV-  PEN  TO  RIGHT  OF  I.D.  FPAMF  AND  SET  INTFNSTTY 
t n  A  L0WEp  LEVEL  fop  pastpr  *r H*  MOVFMENT. 

TF  (m.gt.ci  CALL  KAVANS (1.05,1.25) 


EFT  COUNTERS. 

NGRAFS  r  TOTAL  M'JMPFp  CP  G 9 5  75  H  S  TO  OP  GENPOATFO. 
NOG? AF  =  v U H R E R  OP  NEXT  G3  A  c  h  to  of  CLD  TTFD. 

NGRA  FSsN 
NOG  ? AF=1 

NC»w=L 

icrjri)  =  ni» 

IC I  N  (  2)  =  I  (  2) 

TCD JT(i) =J 
IN'*  (  1 )  =K  ( 1  • 

!»»(?)*< (?) 

V°<  ?)=<(JI 

no I r  p  (5,1) 

1  POR-AT  (  3 1 UI 3L 1 T  HONTINE  Aei.0T  INITTALT7r0) 


ptpST 
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V  )  ■  O' i  f  I  <.ALINr<'O'fYv,N,yM,'">Y*YMt:iYfY,<0r'PS* 

'  -LI  T  A  r  *-515  - 

v'i  ;  v  r~s  of  tat?  patmtS  to  of-  plotti-o.  <  =» f aci 

v>.  -  y  "'i''f:!T>!S(rc  ''r  BAT'  n  0  T  N  T  S  T0  ''c  PLOTTFO.  fP^ALI 

V  r  |  4nrt  or  3AT6  cOINTr  TC  OF  “LOTTES.  (INTEGER) 

•4  ‘if-iTTjc  m  i ‘<o ^ o t _ ^  a  point  oy  point  plot  wtm 
mo  c. tm:>-  ;  tins  ltnee. 


Y“  ~ 

■-X  JT-II'I  "ft'  lr 

pp  U  A  P  T  a  0  L  ” 

or 

ABSCISSA 

.  <  9r  A  L  ) 

1Y  = 

GOAL-  r  a r  t 5 p 

(UNTTO/INSM) 

FOP 

ABSCISS 

a.  <  R -  A  L » 

Y s<  = 

vIlT*!|J4  ’/I  Ur 

or  v  a p’ a  b lc 

TO 

°F  PLOT! 

r0  ON  OFOINaTr  .  (OP'D 

V 

p  0,  4  L  r  F  »  0  T  0  p 

(  I  N  I  T  S  /T  NS  M ) 

OF 

VA PI  4 OLE 

T0  BE  PLOTTZO  ON 

BP  ->  I M  4  T  ~  ,  (5!-Ul 

<  =  AOt-T'rIVir  O3  'It  GA TT  Vc  COD'-.  (INTEGER! 

tf  cosrrrv^,  a  social  c^a^arTFo  f^oi  t-« r  symbol  Tj<BLr  is 

TO  ar  =>\_  x>  T  T  ~  r- .  IMF-  Nil'"'F’  OF  THS  CsA^ftCTE»  IS  GI\*FN  !♦' 

y  ricr  Cl). 

rr  t'rr.ar'-tfE,  T  ~>t.  ir£5  SUPPLTfrS  TH£  PLOT  CODE  c09  MTS 

owv  tD-ojaL  rwa=>ar.TFt  in  <oscs. 

crrsej  -  LM'YO-o  rp  w  j  O.JS  BFSC9T0ING  SYMBOL  . 

^nr-|C£(T)  r  5TAF>T  PP  t  V  M  ^ , 

tt  <  =  ( i  >  vo  social  symbol  ts  rn  ,3r  usfo,  aNO 

(?)  FV”*  f»OI'|T  IN  THE  Ta«Ltr  IS  T0  “F  3L0TT-r). 
'MGNITinr  op  <  I 'PI  SATES  POINTS  IN  TABLE  to  BE  3LOTTFn, 

<  =  1,  z.4 fpy  point . 

<  =  ?,  f1/p-~y  pgrTP  on’>|T, 

*  -  u  ^\lrR  Y  TMIC.Q '  POJMT, 

P  T  P  . 

01')  -  NSTO*1  YW  C 1)  ,  Y1'  ( «j  ,<o  O^S  ( 1 1 
Ir  C  N  .  £  "4  .  n  I  EETJPV 

oufr-iF  ptp  CONTINUOUS  LT  *,r. 

TPt'J.LT.Cl)  GO  TO  ?0 

T1!  4  OOMTTVJOIJS  0IRWF-. 

y-  (  <  v/  »  1 1  -  ym»  /p  < 

Yr(TIMll-YM)  /  0  V 

' a l l  PLnT(y.Y,T) 

00  )P  1=1 ,M 

XrrC  YU  (J  )  -VMt  /P< 

Yr(  /  '*  (X)  -Y'O  /  0/ 

I  I  gall  ri  it  ( v,  v ,  ?» 

~H*Cy  F  IF  SPlOTPL  SYMBOL". 

TF  ( Y  .FN . C  )  GO  TO  5  : 

uo,/  r  P(->J  T0  p  t  5  ^  t  POINT, 

n  Y;(VII(<|  - Vf)  /s  Y 
V  -(  T  )'  (  1  I  -V“)  /OY 
Pi'Ll  01  0 T  (  X ,  V  ,  1 1 

"tT  l~0(M) 

sur-y  p-ii.  (IS vt  oi!PPLTFO  "O0c. 

”r  (y.LT.G)  GO  TO  <=>  0 
L=< 

MO')"  r-YOPFP(i) 

•>  ■>  TL  r  L 

0"  a r  I  =  IL,M,  . 
v-(  Y )/  (  J  I  -y  •'»  /n< 

Vs|YW(fi -vMt /OY 

un  0- L.  YAO yu (X , Y  ,  x . 1 , yoor G (T< ) , r . 0  A <00F , T( 

L  t  p  t  or  \> 


SfLTNF 
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<^3  C.U-  PLOT(X,Y,«» 

or|  jDfj 

USE  R  SIJPPLTFD  'OOF  . 
e.n  KOD- =-20-KOOFS<l» 

IXr? 

Lr  *  < 

n  r o 
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IKK  >0  V  >V  J  (KK)  OCIO  UU.O  vKKKKKK)()(KKHKKK)  J(KKKKHK)0vKK)O'K)C)O 


sui r  out  r me  <al^ cn,  at  ,no,  j,flMiM,ofl,nEL  i 

SU“  3  OUT  T  ME  KALE  PROVIDES  SCALING  INFORMATION  FOR 
THr  AVIS  °  OU  T  T  N  "  KAVTS. 

CALLING  ARGUMENT  5  - 

TN®  J  T  - 

0  =  APOAY  Oc  data  POINTS  to  OF  EXAMINEO.  (RE4L ) 

41  =  AVIS  LTNGTH,  TNCHFS.  (REAL) 

NP  =  NUMBER  0r  PO  TNTS  TO  ME  SC4NNE0  I K  ARRAY.  ( INTEGER ) 

J  =  JUMP  CODE  WHOSE  MAGNITUDE  DETERMINES 

T Mc  D4T4  TO  OF  SCALED. 

J  =  1.  USE  EVERY  POINT. 

J  -  2,  USE  EVFpy  SECOND  POINT. 

U  =  3,  JSF  E',CRY  THRID  POINT. 

ETC. 

(INTEGFR) 

OUTPUT  - 

AMT-J  =  STARTING  value  for  VAPIAOLE  on  AXIS.  (REAL) 

DA  =  SCALE  factor,  UNITS  PER  INCH,  TO  OE  USED  ON  AXIS.  (REAL) 
DEL  -  UNITS  PER  TIC  MARK.  ( Rc  AL) 

ALLOWAOLE  OP'TSIONS  -  1.0, 2.0,5. 0  *  POWER  OF  TEN. 

C  A  J  ”  I  (in  -  WHEN  “LOTTING  ON  GRAPH  PAPE®,  PLEASE  NOTE  THAT 

THE  NUMBER  OE  MAJO=>  TIC  MARKS  IS  DETERMINED  RY  THE  OATA. 
UNLIKE  OTHER  PLOT  PACKAGES  (SUCH  AS  CALCOMD,  ftC.) 

THIS  ROUTINE  wtli  NOT  ALWAYS  PRODUCE  ONE  MAJOR 
TIC  MARK  PER  INCH. 

DIMENSION  0(1) 

on  r  5  T  n  TOTAL  NJHRFR  OF  “OINTS  IN  ARRAY,  N. 

N=N“* J 

A mj *!  =  MINIMUM  VALUE  OF  OATA  IN  ARRAY. 

AMA  <  =  mA  MIML/m  VALUE  OF  OATA  IN  AROAY. 

AHT'jrO(l) 

AMAV=AMIN 
DO  2 P  1=1, N,J 

T  =  TEST  VALUE. 

T=D( I) 

IF  (AMIN.LT.T)  go  TO  10 
A  HI  N  =  T 
GO  Tp  20 

10  IF  (AMAX.GT.tj  so  TO  20 
A'*AX  =T 
?C  CONI  TMUr 

IF  ( AmAX ,EO. AMIN)  AHAV=1.0C1*AMIN 
JF  < AMA X.LT. AMIN)  AMAX  =  0  .99* AMIN 
IF  ( AMAX.NE. AMIN)  GO  TO  25 
NRI T  f  (6,1)  AL,NP 

1  format  (64H1***“0SSI“LE  E®ROR  OFTECTED  3Y  KALE.  OATA  TO  ®E  SCALED 
1  ALL  7ER3.  /15H  AXIS  LENGTH  =  E13.5,1H,T6,8H  POINTS.) 

SPECIAL  OF  RUG  OUTPUT  FOR  CDC  6600  . 

I=L1PF(D) 

WRITE  (6,2)  I 

2  FORMAT  (2flH  STARTING  AD0RFSS  OF  DATA  =  08/) 


AMA V=l. 0 

MINIMUM  AMD  MAXIMUM  SHOULO  FALL  ON  TIC  MARKS. 

ALL3WADLF  TIC  MARKS  =  1,  2,  0“  5. 

R  =  PANGF. 

kale 
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o- i  4  f.  y-a  mtn 
-  P/iL 

TM  -  atOGJ  p  (POL> 

T-I'JT  (0UM> 

W_  -  ’PL/  (1C.  0**T  > 

T5  CPOL.GT.1.5)  GO  to  yo 
I>rLO.  0.:**T 
go  in 

11  TP  (  POL.  r,T  .  2.  (1)  GO  TO  4  0 
orL=?.c»ic.o**r 
go  Tn  sc 

UQ  IP  (OOU.GT.S.C*  GO  T0  45 
r)rLT^.0*lC.3**I 
GO  rO  50 

43  0"L-1C.O**(I+1> 

50  DIM  =  t  *IN/DEL 

I-I',IT  (QI|Mt 

tp  ( cloit(ii .r,r. oun  i=i-i 

A'4I'<  =  rL0flT  (I)  *0"L 

•p=A4A.  Y-4MIN 

O'M  =  l,i  OC1*OE_**L  ♦  aMI‘’ 

I F  C  iMAX.LT  •  0  U  M  *  446  Y  =  OU-1 

MT  r  MU'lSPR  05  T  TO  Mfip<0  ’EY  ON  O  GTA-tTIM3  VALUE. 

N»=  T  *'T(  3/DPL-O  .  3  C  0  1>  ♦  1 
TP  (  MT,  3  T  .  JMT  <  A'_  1  )  GO  TO  *>5 
Ofl  =  0  PL*  P  LOAT (MT> /AL 

ocT  )rN 
EN0 

<HE 
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SU-A-  nUT  I  Mc  <AN'J«  (X,V,H,F,ANGLE,TFMT,NC,N3) 

RUT  3  CUT  T  NF  KAMUM  OLnTS  NU^gE3?  ACCORDING  TO  FORMAT  IFvr. 


r^U-  T  “G  ARGUMENT  2  - 


u 

F 

A  M',  '_  c 
IF«(*’ 


►r 

n3 


rCCRnTNATES  OP  LOWER  LECT  CORNO°  OF  NUMBER 

to  re  dotted,  in  inches,  from  the  origin.  <°pa_) 

HEIGHT,  TM  INCHES,  OF  PLOTTFO  NUMTE3. 

NlJMt3F°  TO  qr  cLoTTF0. 

ran  IF  EITHER  INTEGER  03  OFC  TMA  L. 

anojlar  ortentatton,  in  degrees,  of  plotted  numofr 
■l^ftS'JCEO  COUNTERCL  JCITWISE  PROM  HOP  I7  ONT  4  L  AXTS.(RFAL» 
Pt)»MAT  JNDE3  WHICH  NUMqcp,  F,  IS  TO  OF  WRITTEN. 

AMV  ALLOWS  3  L  F  F03TRAN  FORMAT  THAT  CORRESPONDS 
TN  TVP?  TO  F.  (HOLLERITH) 

NtJMOF3  OF  CHARACTERS  TN  F 03 MAT,  IFMT.  (INTEGER) 

N!)M OE 3  OP  CHARACTERS  RESULTING  from  format.  (E.G.  El<=>.<* 
3ESJLT5  IN  IF  CHARACTERS.)  (INTEGER) 


C  nM  M  ON/<  AC  OM/MC  P W .  ICIN<  ?) , I C OUT ( 1 ) , I NP ( 0) , NO3 Ac 5 , NOG RAF, N CRT 
OTM-M<;nM  MFHT(IF)  ,  TF*n  (  L)  , M0UT(1»1 
o  A  T  A  NFMT(i)/AH(  /  ,  NF  MT  ( 1  M  )  /  <*H )  / 

IF  fMC.EO.O)  GO  to  20 

N  =  MUMTE  OP  ALPHAMERIC  WORDS  IN  IFMT. 


N=N' /NC HI 

IF  f M*NCPW.LT.N0)  N  =  N  ♦  1 
J=’ 


CILL  FORMAT  ARRAY. 

on  l f  T  =  l ,  N 
NFMT  r  J> -IFMT (T) 

J-  J  +  i 

IF  (J.LT.1F)  MPMT( J»=NFMT<13) 

CONCEPT  TO  ALPHAMERIC. 

WRITE  (!*0,NPMT|  E 
R-WTMO  U  o 
NrN? / NC  PW 

IF  < M«N0PW .LT .N3 ) N  =  N+ 1 

pc-at  (<*0 ,  TNP)  (NOUT(T),T  =  l,N) 

oc  jfj  Mn  ur 

PLOT  CODE. 

CALL  <43 YM (X,v,H,N OUT, ANGLE, NR, 3) 

RCT  )ON 

TNO 


<ANUM 


V)(KK)(KK)A  »C1.)OOJOUUU.)OJ'.A)OUl.KK)l)^KnOJUU*  >OC?.KK>OCKKK  )U.  )U.>OC>UsKHKKK) 


SM3RoUTTNF  HEIGHT,  IJCD, ANGLE.NCHAR, IPfNC) 

SUBROUTTNP  K  A  SY  M  DRAWS  Al  OHNUMrftlC  INEO°mATION 
ANT  SPECIAL  SYMBOL5;  at  ant  angle  ANT  ANY  SI7P 
ON  the  PLOTTING  A=>tA. 

TN  ITS  STANDARD  EO.RM  SYMBOL  IG  USED  TO  PRINT  T-XT 

mat-ptal  along  axes,  and  ’■itle  graphs,  the  letters 
\  TH=OUGH  7,  DIGITS  0  through  9  AND  ALL  OTHER  STANOAPD 
EORTPAN  CHARACTERS  ARE  AVAILABLE. 

IN  AN  ALTERNATE  FORM,  IT  MILL  DLOT  A  CENTERED  SYMBOL 
FROM  The  S°ECTAL  SYMBOL  TABLE,  OR  A  USER  SUPPLIED  SYMBOL. 

THE  A°GUMFNT  NCHAp  EBNTROLS  THE  FORM  TO  RE  USFO. 

ARGUMENT  L  TSt  - 


XT,  v  T  r 


ANG.F 


NC  H  A  R 


TPENC 


:  COOROINATES.TN  TNCHE  S , OF  FI THCR 
(1*  THE  LOWER  LEET  CORNER  OF  THE 

e T R S T  CHARACTER  OF  STANEIARU  TEXT  TO 
r»E  PRINTED  -  OR  - 

(?)  THE  CENTER  Oc  THE  SPECIAL  SYMBOL 
TO  RF  PLOTTED. 

:  HEIGHT,  IN  PITHEP  INCHFS  OR  CENTIMETER;. 

OE  THE  TEXT  CHARACTERS,  DR  SPECIAL  synROLS 
TO  IF  PLOTTE0. 

IF  HEIGHT  T3  GREATER  THAN  0.0,  the  UNITS 
ARE  ASSURER  TO  Dr  TN  INCHES. 

TF  HPISHT  IS  LcSS  than  o.O,  T Hc  UNITS 
ARE  ASSUMED  TO  RE  IN  CENTIMETERS. 

TEXT  MATERIAL  TO  RE  P°  TNT  FD ,  OR 
SPECIAL  SXMTOL  COOC ,  DEPENDING  U°ON 
THE  VA'.UF  OF  NCHAR. 

STANDARD  ALPHANUMERIC  TEXT  (BOO  OR 
A-TYPP  FORMAT)  SHOULD  be  L EF T- JU S TI F IEO 
AND  CONTIGUOUS  IN 
(1)  A  SINGLP  VA<?IA9Lr, 

(?)  AN  ARRAY,  OR 

(?)  A  HOLLERITH  LITERAL  (IF  THE  COMPILER  PEPMITS). 

SPECIAL  SYMBOL  CODES  SHOULD  RF  INTEGER  FORMAT. 

:  ANGULAR  ORIENTATION.  IN  OFGREES,  AT  WHICH  CHARACTER 
LINE  IS  to  DF  PRINTED.  POSITIVE  VALUES  APE  MEASURED 
COUNTERCLOCKWISE  FROM  THE  HORIZONTAL  X-AXIS,  AND 
NEGATIVE  VALUES  ARE  MEASURED  CLOCKWISE. 

:  CONTROL  CONSTANT  FOP.  TYPE  OF  MATERIAL  TO  RF 
PRINT'D  ON  PLOT. 

A  POSITIVE  NCHAR  INDICATES  STANOARO 
A  LDH  A  NUMERIC  T^XT  IS  TO  9F  PRINTED.  NCMAR 
TS  THE  NUMBER  OF  ChARACTFpS  TO  BE  PRINTED. 

A  NEGATIVE  NCHAR  INDICATES  A  SPECIAL 
CENTERED  SYMBOL  IS  TO  BE  PLOTTED.  IE  THE 
ABSOLJTE  VALUE  OF  NCHAR  IS  BETWEEN  1  AND  14, 

A  SPECIAL  SYMBOL  FROM  THE  TABLE  IS  orinteo. 

IE  THE  ABSOLUTE  VALUE  OE  NCHAR  IS  GREATER 
THAN  14,  THE  USER  SUPPLIES  HIS  OWN 
SPECIAL  SYMBOL  IN  I8C0. 

IE  MORE  THAN  TWO  WORDS  ARE  REQUIRED  FDR  the 
DEFINITION  OF  THC  NEW  SYMROL,  EXPAND 
THE  DIMENSION  OF  NOD  . 

THE  DIMENSION  OF  THE  NEW  SYMROL  IN  THE 
CALLING  PROGRAM  IN  NOW  ASSUMED  TO  3E  TWO. 

PEN  CODE. 

?,  PEN  DOWN  DURING  MOVE  TO  (X3.Y3). 

3,  °PN  'JP  DURING  MOVc  TO  (X3,Y3>. 
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►imn  concs  for  standaro  pqrtran  chapactezs  a?? 

OnNTAINED  IN  NORMAL,  ANO  cPF0  T  AL  CYM90LS  IN  NOOES. 

NSYM  TS  the  N!.mE°  OF  SPECIAL  c Y  MHOL  CODEC  IN  NOOFS. 

COMNON/K ACOM/NC3w,  ICIN( 2 ) ,  ICOUTCl » , INP( 3) ,  Nr,=  AFC,  NOGRA*.  NCRT 
DIMr  NS I  ON  ICHARS (10' 

DIM-  ms.T  ON  T3C0(  1  >  ,  IN(  li.  ),  tnoph  (SI)  ,  NODES  (90)  ,  NORMAL  ( 205  1 
OT H '  MSION  NORN (51 ) 

0  AM  vsa  VE/3  .  0/.  YSAWE/0  .  0/ 

0  4  T  a  NOaMAL/5713, 1757. 5744, 62 03,0757, 

11744, 6253,132?, 2110, 6667, a7?6,  440-, 6263,1722,  2110, 

2501n, 2526, 1757, 6644, 6380, 2040,1717, 4400  ,2387,8310, 

3 25? 8,  175  7, 88  44, 2  21  3,  5  36  2,8657,  1728,  2510, 0  0  55, 6644  , 

46253, 2301,07 44, 5 086,665 7, 1726, 25 10, 5061, 625 7. 172 2, 

52110  , 44 0  0,2110,5  061,6253, 1322,  2617,  5766,4400 , 5310, 

66044, 44  CO, 0607, 5786, 0  64  4, 2 763, 6101,  6167, 44 00, 5713, 

70717,57 17, 44 0  0, 2213, 5762, 6657, 1726, 20  GO, 05 44, 6766, 

05717,  76 2 t,  4400, 6765, 1722, 2110, 6000, 27 44, 23 67, 60  CO, 

86 057,  27  44, 0206, 0  050, 1044, 170 3, 575o,  0717, 4400, 6303, 
11218,0757,4400,5121,4026,6544,6267,2322,4003,0744, 

26 76 3,  0727, 2744, 2217, 5 7 62, 6150,  1025,2617  ,5766,4400, 

76 761, 0321, 2725, 6544, 6 723,4400, 2 61 7, 5 766 , 6 25 7 , 1 37 2 , 

44400 , 6387, 2726,4400,6767,8857, 1726, 2213, 5340, 2357, 

5440.  ,  67 63, 2723,4400, 6 32  7, 402  3, 6744, 6300, 2745, 0007, 

6 440r, 63 23, 67 27, 4 05U, 1044, 6G1C, 212 2, 13 67, 6717, 26 25, 10 44, 6763, 
7 132 7, ?81 7, 6744, 636 7, 6 020, a7?3, 4 400, 66 67, 0707, 6773, 4400, 6767, 
66677 ,60? 7, 4400,6 010, 212 2, 1363, 6744, 17 26, 2? 17, 5 35 2, 6657, 1740, 
9162  7 ,44  00,63 07,  2  344, 6 T57, 0317,  2344,  66  23,4  003 , 074  0,63  27,4400  , 
1  1 60 5  ,  3615, 1607, 4400, 2 21  3, 57 62, 615 0,1025, 2617, 5766, 400 3,  07 44/ 

Tara  morn  /iho, ihi , ih?, ih7, ih4 ,ih5, ih6, 

11H7,  1  H6, 1HQ, 1H-, 1H  ,iw. , 1HF,  1HI,1HG,  1HU,1MP,  1HF, 

21HV ,  m(  , 1H) , 1H= , 1HT  ,  1HM ,  IMS, 1HA, 1H/,  1HC* 1HL,  1H3, 

3  1 H  N  •  1  HX,  1HY,  1 H  Z  ,  lHr,,lHD,  1HH,  1HJ,  1HK,  1H°, IHO,  1HV, 

41HUI,  1M*,  1H,,1M7,  1HX,1H<,  lHi,lHX/,NN0RM/51/ 

DATA  TNO-?M/l,4.7,l3,2C,24,  29,76,3  9, 

14?  ,  5  6,5  7,58,61 , 5  5, 6  9,  75 ,7'9,  a 4,  83,91 ,95,99, 

21 3  2.  106, 109, 116, 120, 122,  17 7, 130, 137, 140, 14 3,  147, 

7151, 157, 181, 185, 189, 173, 177, 134,1B8, 169, 194, 193, 91, Q5, 3 8, 99/ 
OAT  a  NODES/4033,  7  7  77,  77  73,4044,4053, 

1  7175, 57 1  7,  35 31, 1  35 3, 40 44, 400  3,  773 7,  0  7 40, 4400, 00 0  3, 

70 70  0,  70  30, 0040, 4400, 0073, 77 GO,  3377,  0040, 44 00, 43 03, 

37^0  7  ,  70  0  7 ,4044,  C  O'”  7,0  77  0  ,  7  0  0  3 ,  4  044,  OC77,3373,  770  0  , 

44-44, 4077, 3377, 77 00,6020,0040, 44 00, 0007,00 77, 0033, 

6  004 C ,  44  0  0 . 40  7  3,6  266 ,776  6,2637,  262  2,  3322, 6? 40, 44  00, 

6 n 007,  0  700,  37 77, 0073, 3700, 4044, 003 3,  73 37, 7700, 40 44, 

70037, G3G0,4044/,NMOncS/6?7 
nAT"  IN 71, 6, 11, 15, 20, 25, 2°, 33, 37,43,48,  66, 

16?,6f/,NSvM/14/ 

DATA  LIFT/ 40 7, NEXT/ 44/ 

O A T A  PimaO/1,  746329252F-Q2/ 

CHECK  FOR  EMPTY  RIJN. 

IF  (MCHAR.E3.C)  RETURN 
Xl  =  *  7 

if  t y 7.F0.999. 0)  xi=ysAve 

Y 1  =  Y  3 

IF  I Y 7. E0.999. 0)  Y1=YCAVE 

MOVr  PFN  TO  (XI, Yl»  . 

CALL  pLOT  (XI, Yl.TPFNCJ 

CHANGE  ANGLE  TO  RADIANS. 

T  =  A  N  GLE*  P 1018  C 

initialize, 

EINT=SI N(T) 

Cr»5 1  =COE  ( t  ) 
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H  =  HC  Tf.MT 

change  centimeters  to  inches  tf  necessapy. 

IF  (H.LT.O.O)  H=-H*?.54 

OF'.  =  OISTANCf  IN  INCHES  BFTWFFN  NOOES. 

0FL  =  h/7.  C 
DEL  C  =OEL  *  COST 
OrL  S=OEL*STNT 

TST4s>t  =  1 

OFTFPMINF  TY«>r  f  CHARACTERS  TO  BE  PLOTTEO. 

tf  (  MCHAO.r.T.  0)  GO  TO  7  0 

SPECIAL  SYMBOL. 

N  =  - N  C  HA  0 

X2  =  X 1 ♦H*0 PST 

T’:r 1 ♦H*  SINT 

IF  tH.GT.NSYH)  30  TO  50 

SYMBOL  IS  CONT  AI Nr  0  IN  T  AnL  F  . 

IM3r  Y  =  IN<N) 

10  ICH'°-NOnES(INO“X) 

TOS  3  T  r  1 

°L  0  T  NOOES. 

20  NO3F=ICHAR/10C 

ic  <  TPflPT.  E3 . 2)  NOOF=TCHAF-10C*NOOF 

IF  \IOOE  =  1,0,  LIFT  FFN  4N0  GO  BSCK  FOR  ANOTHER  »OINT. 

IF  (MOOF.EO.LIFT)  GO  TO  SO 

IF  NOOF  =  44,  END  OF  THIS  SYMBOL. 

IF  < NODE .EQ. NEXT)  GO  TO  170 

SEPARATE  MODE  INTO  (X ,  Y )  COOROINATES. 

IXrNCBE/10 
Tm  =  mOOF-J  3*IX 
IF  (IY.GT.4)  I X  =  4  -  I  X 
TF  (TY.GT.4)  IT  -  4  -  I v 

OBTAIN  COOROINATFS  OF  NOOF, 

X=X1 ♦OELC*FLOAT(  IX)  -DELS ’FLOAT  (IY) 

Y  =  Y  1  ♦ DEL C*  FLOS  T ( I Y ) +OFLS*cLOAT  ( IX ) 

PLOT  °0I NT . 

TF  (  TSTAPT.EQ.l)  CALL  PL0T(X,Y,3> 

IFTAPT=TSTA»T  fl 
T°E  M~? 

TO  CftLL  FLOT  (X,y,ioen> 

OBTAIN  NEXT  POINT. 

IPA  ?t  =  ipaptH 
TF  (TPART.LE.2)  GO  TO  20 
INOF  X  =  l NDrXF 1 
GO  TO  10 
50  Iocn=T 
ISTAPTrl 
X-Xl 

y=ti 

GO  TO  30 

<A*NM 
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(sc  tn-»-  x=nnooeg*i 
J  =  'l-  ?c 

■in  S?  T=  l,  J 
K-M'iOnFS  +  I 
6E  MOD  r  ^  <  < »  =inC0«r» 

r,n  m  n 

stamrard  fortran  io^f. 

70  *■}  w  =  1 

vr  r  ; 

0 *-•  I  r  T  TO  OENTPR  pr  GRIP. 

X2=X1+T.  C*  nEL3-PELS> 
v7-=/i+3.r*(9-L'«-Pt'Lr> 

007  =  7 ,0*PcLC 
7.:*n^LS 
»1  t-TTEC-l 

Tn4  3  T  =  1 

0OT4TM  SOOF  OF  CHARACTER. 

WRITE  (  40,  ICOUTI  Ir'Cr'<NW) 
orjjm  4  0 

RF4T  <40,TCIN>  ( TTHAPSl It ,1  =  1, NFPWI 

REWIND  40 

Rp  Nf  O'*  r  =  TFHflPS  (Lr  TTr  F) 

no  ar  jtfn  =  1*  NNORO 
IF  I  NOOE  ,f  0.  NOR'UI  TEN)  *  OR  TO  06 
06  CON1  T  I  Wt 

WRITE  (6,1)  M0Tc  ,TR(-n  (\'W) 

1  F0R-AT  (  6e.HP**»W  A  IN'INf,**  •.  N  ON  -  ST  AN  0  A  RD  FORTRAN  RVMOOL  rMCO'J*‘>TFRED 
1  n»  k-fiGYM.  A1,2X,A10/) 

RCT J°N 

OBTAIN  StARTING  POINT  TN  T  A  OLE  . 

36  IMO-X  =  INC»M<ITFM> 

ORTiTN  CHARACTER  COOP  FRO"  T  A ’LE « 

ITT  IPHA  »  =  NOR''AL  t  INDEX  I 

IOfl'JT  =  i 

PLTf  NO  OF. 

no  not f  =  r phar/io: 

IF  (  TPART.FO.  ?>  ''ODF=IPHAP-1CO*NOGF 

IF  MOPE  =  40,  L  r  FT  FFN  A  N">  GO  RACK  FOR  ANOTHER  °0INT. 

Tc  < NODE. EO. LIFT >  GO  TO  140 

IP  NOPt  =  44,  END  OF  TM I S  CHAFAOTEP. 

IF  ( MOOE.FQ.NFXT )  GO  =  0  tog 
SPtJAr'ATR  NOOE  INTO  IX, x)  EOOFOI  NATES  • 

IV  *  MOpr/lC 

TV  :  AGOF-10*!* 

IF(j'y.GT,4)  tx  =  4  ~  TX 
TF  (IY,GT.4»  TT  r  4-TY 

rLTT  NOOE , 

XpX?40ElC* FLOAT!  I X  >  -0 F L S *FLri AT  ( I v  I 
Y=f  ?  ♦nFLC*  FLOAT!  I*  >  *P  EL  6*  r  L  R  A  T  ( I X  ) 
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IF  < TSTAPT  .E0.1>  CALL  PLOTtx 

IST5DT=TSTART*1 

IP- 'J  =  ? 

120  CflL.  PLOT  (X.Y.tPFN* 


y  ,3> 


OBTAIN  N  F  X  T  PniNT. 


IPA5^  =  IP  ART* 1 
IF  ( TPAPT.LE. ?)  00  TO  110 
I  NO :  X  - 1  NnFX«- 1 
50  TO  IOC 
1«.0  TPEN  =  3 
X=T  > 
y=y  ? 

IST1PT=1 

50  to  120 


nfxt  CHapflCTEo. 


150  LFTTFP=LEtTE<?M 

CALL  PL0T(X2, V2, 3) 

Nn=vjr*i 
IFT1PT=1 
X?  =  <  2  +  0C7 
Y2=T  ?»DS7 

Tc  (  NO  •  G  T  .  H  A  R )  f,0  TO  15C 
TF  « LETTFP.LE.NCPH>  GO  TO  90 

INCREASE  WORD  COUNT. 

NW=N W*1 
50  TO  flO 

PFS'T  PEN  AT  MEN  <X?,Y2). 

150  XSA«JF  =  X2  -3.GM0ELC  -PELS) 
YSA  i/ r  =Y?  -3.0M0FLC  +OFLS> 
155  CALL  PLOT  (XI,  Yl,  3t 
PFT  JPN 

1*0  XSA;F=X2 
YEA VF=Y2 
50  TO  155 
FMT 
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SUlROUTINr  KAVA'JS(X,Y> 

SUVOUTINc  KAVft  MS  ACVANCFp  T  H r  FILM  TO  THC  NFXT  FPAMF,  AMD 
opS-Tc.  Tnr  L00T3AL  0»TGTN  PREVIOUSLY  DFFINF9. 

4  GALL  WITH  EITHER  X  OR  Y  =  9^9.0  WILL  RESULT  IN  thf 
J  T  OIJFL  y  OEFIWFO  VALUES  OF  X  ANO  Y  9EIUG  USFO. 

data  xsave/o. ox, ysave/o . ox 

Xl  =  x 

IF  ( XI. FO. 999.0)  X 1 =X  S A  V  F 

XOAVFrXl 

Ylrx 

IF  ( Y1.F0.999. 0)  Y 1  =  Y  SA VE 
YRA Vr=wl 

FOR  CALC  0MF  PEN  PL  0TTFRS  .  USE 

CALL  PLOT<  XI  ,  Yi, -■>) 

CALL  PLOT  (X1.Y1.-3) 

FOR  CATMOOF  RAY  T(jqr  PLOTTED,  USE 

CALL  CAlCMP(0.0,C.0,C000,?» 

CALL  CALCMPt  xi, Y 1, -6, 3) 

TO  AnvAMCF  thf  film,  ano 
TO  :’F9ET  thf  LOGICAL  ORIGIN. 

RET  Jon 
FN0 

<6 VANS 
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subroutine  kaxis (x ,*, la  sel.nc,  al,  angle,  amtn.oa.ofl,  itu>n» 
subroutine  <axis  draws  a  labeled  arrs. 
cal.ing  argjmpnts  . 

x  =  X-ODORDINATF  IN  INCHES  OF  THF  START  OF  THE  ORIGIN.  (REAL) 

v  =  Y  -COOROI N  a  T  F  IN  INOFS  OF  the  START  OF  THE  ORIGIN.  (PfJLI 

LA  B  :  L  =  ALPHAMERIC  tpXT  Tp  gr  drintfd  along  axis.  (HOLLERITH) 

N~  =  NUMB PR  gp  CHARACTPB3  IN  LABEL  (  INOL JnIN3  ALL  r»LaN<S». 

ft  POSITIVE  vaujp  WILL  PL  ACE  THF  TITLE  ON  thF 
CPUNTERCLOCKWIS”  RIDE  OP  Thf  AXIS.  USUAL  Y-AXTS. 
a  NF  Ga  T I  Vr  VAL'Jr  NTLL  °LaCr  THP  LaorL  ON  the  CLOrKWISE 
STOE  OF  The  avi3,  USUAL  X-AXTS. 

(INTEGER) 

at.  =  length  of  axis  in  inches.  (read 

ANGlF  =  AM GUL AP  OPTENTATTON  OP  AXTR,  OpGREFS. 

r  DEGREES  =  HOR IZONT  AL  TO  RIGHT,  POSITIVE  ANGLES 
APE  MEASURED  COUNTER CLOCK WISP. 

DEAL) 

AHTM  =  MINI  HUH  VALUP  0C  VARIABLE  ON  AXIS. 

VALUE  RETURNED  nv  <alp.  (RPAL) 

HA  =  PCALr  PASTOR,  UNTTS/INCH,  OP  VARIABLE  ALONG  AXIS.  DEAD 

OrL  -  UNITSZTir:  HftPK,  RETURNED  BY  K^LE.  (REAL  ) 

ITJRN  r  NUHRPR  0C  00  DFG?EE  POTATION  OF  NUMBERS  INOTCATING 
VAL'Jr  OF  MAJOR  TJC  MARKS.  (INTEGER) 

COMMON/KACOH/MSPH, ICIN( ?> , ICOUT (1) , INO(I) ,NGRAPS,NOGPAE,NCPT 
ill  ME  MSI  ON  LABFL  (  1 )  .  JFMT  (  1?)  ,  IPMT  (  3) 

RATA  JFMT/t*HF6.1,l»HF6.2t4HFf,.3,2HI?  ,  UH,  1HX  ,  4H  2H1Q  ,  ?HI  3  , 

12HI5,2HI6, 2HI4, 2HI 3,2HT1/ 

HGTNUm  =  HEIGHT  OF  NUMBERS  ALONG  AXIS.  INCHES. 

OSTC  =  DISTANCE  BFTWEFM  AXIS  ANO  TOP  OF  NUMBERS,  AND  3FTWEPN 
NUMBERS  ANO  AXIS  LABEL. 

H  G  T  ~  v  M  =  MAXIMUM  ALLO'JARLP  H  FI  GHT  OF  SYMBOLS  IN  AXIS  TITLE. 
TNCHES. 

OATA  HGTNU«/t).10  0/,DSTO/0.0-J/,HGTSVM/o.  130/ 


CFNTER  -  NUMBER  CENTERING  VARIABLE. 

°p  0  1  L  FM  - 

S01r  SYS  TP  MS  RIGHT  JUSTTPY  NUMBERS  WRITTEN  UNDER  FORMAT 
STATEMENTS.  OTHERS  LEFT  JUSTIFY  THFSE  NUMBERS.  THUS,  ON  SOME 
MACHINES  T HE  MAJOR  TIC  Mfto<  VALUES  MAY  BE  SHIFTED  ETTH-’p 
LEPP  OR  RIGHT  Op  WHERE  Tnry  SHOULD  IE. 

SOLJTTON  - 

C  E  N  t  p  R  PERMITS  CENTERING  "UMBERS  FOR  PARTICULAR  SYSTEM  IN  USE. 

rrNTPP  =  C.o,  LEFT  SIOF  Oc  NUMBER  WILL  LINE  UP  WITH  MAJOR  TIC  MARX 
R  rN  T  PR  =  1.0,  RIGHT  SIDE  OF  NUMBER  WILL  LINE  U=»  WITH  MAJOR  TTC  MAR 
C  PN r  p  R  =  0.5,  MTODLE  OF  NUMBER  (AS  FORMATTED)  WILL  LINE  UP  WTTH 
MAJOR  TIC  Mao<.  THUS,  WITH  AN  15  FORMAT, 

T H p  THIRO  SPACE  WILL  LINP  UR  WITH  Thf  tic  MARK. 

OATA  CENTFp/0.5/ 

SIDE =1. 0 

IF  (NC.LT.O)  STOF=-i.O 
OSTz  OSTC 

IF  (  SIDE .EO. -1.0 .AMD. ITURN.FO. 0)  DS T  =  DST+HGTNUM 
CTT :  r  ,  0 

if  (ITURN.NE.0)  CITa-1.0 
XOB  = 1,0 

if  ( SIDE. EO. FLOAT ( ITURN) )  XOR=0.0 
Turr  A  =  A  NGLE*  0. 0174532 °5 2 
SINT sSIN(THFTA) 

POST  =COS (THETA) 

T-:  ANGLE*  EL  OAT  (TTJRN)*90.  C 

KAXIS 
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-  C.  'lru**??1*? 

“  YM  -  r  J  ►)  (  T  M  ) 

r;  ipifntation  !’n-?^rf'Tio'Jc'i 

''■'Jr~I''-»Y0->*'-i‘O*CTT-;EN^E7?*C0ST*frTT«-t.rM 

'"ni  =  ;t  )"*xn=*;orT*pirfr."“jTt^»srNT»(ciT  +  i.; 

plot  rir  harks. 

O  z  C’A  JC,~  IF  OTJ  n|_F, 

M*  r  'rj'll'O  OF  FIT  UJ=<F. 

A”lY=Li«Fk|4.j^»rij 

-  -  A  ‘.1F-1MT  *4 

M7= ' M  T  < P/OFL  +  3 . 1 > 

y>:< 

v?-y 

TT~  mkT'ut  -  3.1  INCHES. 

xa  =  f  *■  .  i  *sn«F*si  i>  t 

Yf  =v  -3.  i*snF‘I1ST 

<F  r  CO ‘1ST  ANT  INC  R  F  MF  NT  ALONG  X  AXTS. 

Y  C  -  CONSTANT  INC^rNFNT  Ai_OMG  Y  AXIS. 

Y  C  :  ■>  r  I  »OOST/OA 

Y  " .  r  '  F  l  *  S  T  V  T  /  0  A 

r  V_  PI  IT  (  XA  ,  VA  ,  X  l 

■'ALL  cLOT  (XA.VA,  ?) 

J-l 

OP  i r  1=5 , NT  M 
CAUL  PLOT ( Y3, vs, ?) 

y-l-Y  4  y  ' 

v  T  :  /  n  *■  V  C 

sal:  clottyb, YT, ?i 

x  ".  =  v  a  4-x  c 

V.'irr  v  A  fY"' 

SAUL  CL',T(XA,YA,  ?) 

PLOT  ►'TM!'3  Ftr  MAP<%  1I/INTF5VAL. 

N  r  ‘  -  1C* 'IT 

Tc  <  J.FQ . ?)  GO  TO 

J  =  ? 

XA=yR+G.0G»STOE»STMT 
.  0F»FinF»','v;T 

xr=-r.i»vr 

YCr-C.l*yC 
G  •'  "  n  5 
Yfs-  13.  .J*xr 

VF;- 1 ' ,F*YC 

MFI’F  \J  ALU-  AT  TIC  ''AF<T. 

NT=NT*1 
A  -  A  F  ( 0  c  L  1 

tf  (a.lt.i.q.ob.a.gf .tOGor.CT  go  to  so 

I'T-TT  format. 

I  =  r  « T (AMT") 

J=T'T  (,)FL) 

rr.  i 

T  F  =  ' 

TFMf  OrjFXTMl 

IF  <ff  •TN.NF.-10ai0.rj>  GO  TO  ?L 

TFM4(J|=JFut(0) 

GO  r  o  ->s 

o  =  A  ’  ~  1  .i  •'  A  X  ) 


j-e  ro 


AD-A074  887 


UNCLASSIFIED 


KAMAN  AVIDYNE  BURLINGTON  MA 

trap-ml-a  two  dimensional  thermal  response  code 

NOV  78  W  N  LEE 

KA-TR-154  DNA-4770F 


F/6  20/13 

TAILORED  FOR  TH~ETC(U) 
DNA001-78-C-0057 
NL 


nno  otinci  ono 


IF  IO.LT.A0S  (AMIN))  3  =  A  3  S  (  A  M  T  N  ) 

3=3»0.1 

IF  ( O.Gr.  1  POO. )  GO  TO  25 
I’s'* 

ifmt  id  =  jfmt no) 

TF  ll.GE.  100.0)  50  TO  25 

T“=T 

IFMT  I  1)  =  jfmt  (11) 

Ir  (AMIM.LT. 0.0  .OR.  DEL. LT. 0.0)  50  TO  25 

T»:  p 

TFMT (1)=JCMT(U) 

IF  (  ".GE.1G.0)  3  0  TO  ?5 
IFMI ( 1) - jfmt (12) 

I=>=1 

?5  FIP=FL0AT(IR)*H3TMUW 

Xr’=X-FI0E*DST*SINT*.FIo«00M 
Y*=<  *STnE*OST*OOST  -FI^’CO'ii 
no  TC  11  =  1, NT 

CAL  _  KANUMIX8»Y3,HGTNU“, I,T, IFMT,IC,IR> 

Y0=)(O4.y~ 

Yq  =  Y04YC 
TC  1=10 

<♦0  IF  (NC.fq.oj  ro  TO  <,5 

LAV L  AXIS. 

S  =  cLnaT(Nr)*snr«.i.n 
M=4L/S 

IF  (  H.GT .HGTSYM)  H  =  HG  TS  Y M 
T0=H*S 

X=»=Y*0.5*  (  AL-TO)  *  COST-SI  OF*  SI  NT*  < HGTNUM*2 . 0*OST0 *  0.  5*  <1 . 0  -SI^E)  *  H) 
IF  ( ITU^N.NE. 0)  XO=XO-( FIR+OSTQ) *CTN 

YB=v«-0.5*  (AL-TQ)  * S IHT *S IOE*COS T*  (MGTNU*«-2.  0*  OSTO*  C.  5*  <  t .  0  ~SinEI  *Hi 

IF  ( TTUPN.MF.Q)  Y=A=Ya-(FID  +  osT0)*STN 

T=IAOS(NCI 

call  <ASYM(X3,Y3,H,LAOEL,AMGLF, I, J) 

<♦5  ’’F  T  J  o  N 


or-iviflL  NOTATION. 

53  0=4.0010 (A) 

IF  (  p  .GT  .0.0)  3  =  01-0.02 
I=INT  (0-0.01) 

IF  (I.GT.C.0P.I.LT.-2)  GO  TO  80 
IFMr (1)=JFMT(TI 

IF  (*».LE.-1.0>  I FMT (1 )  =  JFMT ( 2) 

TO  =  *4 

T3s=> 

A  =  A  M  T  n 

FI>=5.0*HGTNUM 
X0=<-SIOFTOST*SINT+FIO*CON 
Y3=v»5IOE*PST*COST-FIP*CON1 
70  OO  7?  1=1, MT 

call  <ANUM(X0,Y3,HGTNUM, a,t,ifmt,ic,ir» 
X3  =  X °  4XC 
Y^  =  Y  R4YC 
72  a=440FL 
GO  TO  40 

EX»3NFNTIAL  F03MAT. 


80 


Q0 


4=10 .P**T 
4M=4MIN/4 
404=OFL/A 
NT  =  N  T  -1 

IF  ( AOA.GF.0.01) 
1  =  1-1 
GO  TO  80 
IFMT ( 1) =  JFMT ( 2) 
IF  (AOa.GF.i.rj) 
IC=  '* 


GO  TO  °0 

IFmt( d  s  JFMT (1) 
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KAXIS 


ot*i>  ;k)u 


CI’  =  F  .  cr  HGTNUM 
X  ■»  i  x  -  GI  0  F  *  05  T  *  SI  N  T  ♦  F  T  »*  0  no 
YO  =  y»0T0F*QST*3O5T-FTP»OP‘ii 
lOOnnijr  j=i,MT 

CAL-  <fl  JU*<  <X°,  *3  ,  ur,  tmijm  JM  r  ,IFMT,rc,TP» 
X3=X"+X5  ’ 

Y=)=YO«-YO 
110  A 

IF'IT  <11  =  JFMT  <?) 

IF'IT  <  ?>  r  JTMT  (<51 
IFir  <  7>  =  JFMT  <P,> 

Tr-1  ?. 

IP=1 7 

X3=x  "-0.4YCTN 
YO=YO-O.U*STO 

115  CAL'.  K'AMUM<X^,Y3,HGTNU~(,A«,TtIFMTfrC,  I’l 

TN3-PT  r YPCMEMr. 

0*HGTMJM»rTr'-HGTNtri*ST^ 

Y3=  YG»Q.  0*HGTND'1*STm-Hr,  TNU*'*CTN 

3.0  =  F  0HARACr~95  TF  f j  > M -i £ f?  «•  }  FO*  X1G. 

TF'ir  <l)  =  JFMT<t») 

IFr  ? 

I=>=? 

IF  <T.LF.°.AND.T.GF.-o>  GO  TO  117 

pr  7 

IFMT  <i) = JFMT (7) 

117  CAL.  YAOU*  <XO,  y ° , HGTNU^.I , T, IFMT ,TG , IP) 
GO  fO  ltd 

fmo 


KAXIS 
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